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INTEGRAL EQUATION FORMULATIONS OF SCATTERING 
FROM TWO-DIMENSIONAL INHOMOGENEITIES IN A CONDUCTIVE EARTH 

John Robert Parry 
ABSTRACT 

The electromagnetic fields scattered by finitely conducting cylinders 

of arbitrary cross section in a conductive half-space with an arbitrary 

earth-air profile are calculated. The following conclusions of importance 

to geophysical exploration are reached: 1) The ratio of real (H) to 

imaginary (H) is a function of traverse position x and ground conductivity 

as well as the cylinder conductivity*^ ; 2) Topography can give rise 

to a tilt angle of about 5° at an operating frequency of lOOOhz., ground 
-3 

conductivity of 10 mhos/m and normal incidence of an Ey-polarized plane 
wave; and 3) In no case was a zero phase observed, even for perfectly 
conducting scatterers. 

i 

The problem is formulated by choosing an integral representation for 

the electromagnetic fields in each homogeneous region present. By enforc- 

! 

mg the boundary conditions on tangential E and H, a set of coupled inte- 
gral equations results which can be solved numerically for the unknown 
equivalent surface current densities on the interface bounding each homo- 
geneous region. Once these current densities have been estimated, the 
fields can be calculated at any point from the general integral represen- 
tations. 

The validity and accuracy of the integral equations are demonstrated 
by comparing numerical results with analytical results for scattering 
from circular cylinders in a conductive whole-space assuming plane wave 

iv 



and line source incident fields. It is shown that three-dimensional source 
configurations can be considered by expanding the primary current distribu- 
tion and the field it radiates into a Fourier integral over a continuous 
mode distribution. 
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CHAPTER 1- 


INTRODUCTION 


The determination of the electromagnetic fields scattered by conducting 
inhomogeneities reduces to the solution of a vector Helmholtz equation. 
Unfortunately, the classical method of solution yields analytical results 
only for those inhomogeneities whose surfaces coincide with coordinate 
surfaces of orthogonal coordinate systems in which the vector Helmholtz 
equation is solvable by the method of separation of variables. In these 
special cases, which include wedges, spheres and circular, elliptic and 
parabolic cylinders, the fields in each homogeneous region are expressed 
as a sum of mode functions, and the unknown constants are determined by 
matching boundary conditions. This theory has led to useful solutions for 
scattering from objects of simple geometry in a conductiv e who le-space when 
the characteristic dimensions of the inhomogeneity do not exceed' about 25 
wavelerigths. In the past, the important geophysical problems of scattering 
from multiple conductors, conductors of arbitrary shape, and conductors in 
a conductive half-space have received little attention because of the 
intractability of their formal solution. 

The possibility of a general numerical solution of antenna and 
scattering problems was suggested first about ten years ago by Sinclair 
(1959), and was based upon an integral equation approach investigated by 
Albert and Synge (1948) and Synge (1948). In this short paper, Sinclair 
pointed out that with the availability of high speed digital computers 
it is possible to obtain numerical solutions of electromagnetic boundary 
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value problems where suitable analytical solutions are often difficult to 
obtain. He notes that while . . analytical solutions may be aesthetically 

t 

more pleasing than numerical solutions, . . . this is of little interest to 
the engineer concerned with the design' of a particular electromagnetic 
system. " 

Of perhaps more importance is phat he notes also that . . even 
when analytical solutions can be obtained for electromagnetic problems, it 
is usually found that a great amount of calculation is required to obtain 
numerical values, as, for example, when .wave 1 functions in a given coordinate 
system have to be computed. In many cases, direct numerical solution of a 
suitable differential or integral equation can yield the same accuracy of 
result with little, if any, increase in the amount of computing required.'-' 

In particular, we might note that D 'Yakonov (1959 a,b) has published a 
solution to a circular cylinder and a sphere in a conductive half-space, 
and we have yet to witness numerical results for either of these particular 
problems . 

Some of the first numerical results of electromagnetic scattering 
problems were presented by Mei and Van Bladel (1963 a,b) and Andreasen 
(1964, 1965 a,b) . Later published work includes significant papers by 
Mitzner (1967, 1968), In these papers, the electromagnetic fields exterior 
to the mhomogeneities are written m terras of integrals around the contour 
bounding each mhomogeneity and in terms of a surface impedance boundary 
condition. These integrals involve the two-dimensional Green's function 
and equivalent electric and magnetic surface current densities at the 
boundary of each mhomogeneity. The unknown equivalent current densities 
are estimated by enforcing the electromagnetic boundary conditions at a 
finite number of points on the boundary of each inhomogeneity. In so 
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doing, a number of sampled values of the unknown current distribution 
is determined, from which the scattered fields can be calculated at any 
point using the general integral representations. 

In the present study, integral representations are chosen inside the 
scatterer as well as exterior to the scatterer. This approach removes 
restrictions upon the conductivity and curvature of the scatterer inherent 
in the above formulations, although it leads to more complicated integral 
equations. However, it leads also to a consideration of the problem of 
scattering from cylinders in a conductive half- space. 

Since the coutour can be deformed into any shape desired, it is 
possible to formulate and solve any scattering problem within the storage 
and time limitations of a computer solution. Although this study will be 
restricted to scattering from two-dimensional inhomogeneities, it is shown 
that the incident field is arbitrary. 

The purpose of this dissertation is to consider the theory of integral 
representations as applied to solving two-dimensional geophysical scattering 
problems. In particular, examples will be given to demonstrate the validity 
and generality of this approach. However, due to the length of this subject, 
a general numerical analysis of geophysical scattering problems must be left 


to a later study. 



CHAPTER 2 


DERIVATION OF INTEGRAL REPRESENTATIONS FOR E AND H 

Sliver (1965; p. 201) has shown that if a two-dimensional waveguide 
is homogeneous in structure along the axial direction; only two independent 
field components; E y and Hy; exist; and all other quantities may be derived 
from these. Since similar conclusions can be drawn for two-dimensional 
scattering problems; the first task in treatment of scattering problems 
is to obtain general representations of the axial field components. 

Integral expressions for E y and were obtained by Papas (1950) 
and Borgnis and Papas (1955) in investigating scattering from circular 
cylinders of infinite conductivity. Andreasen (1965 b) showed that their 
work could be modified to general integral representations for E y and Hy. 

In this paper; Andreasen assumed that the surface impedance of the mhomo- 
geneity was very small with respect to the impedance of the exterior region 
so that his final integral representations are not sufficiently general 
for earth scattering problems. 

We will follow Andreasen's development to obtain general integral 
representations for E and H interior and exterior to the cylinder. In 
section (2-5); it will be shown that the general integral representations 
reduce to those given by Andreasen if it is assumed that the scatterer 
has a low surface impedance and tint the skin depth is much smaller than 
the radius of curvature. Although the integral representations m the 
exterior region are modified only slightly by making these approximations; 
the advantage obtained is that the solution of the resulting integral 
equations is much easier. However; if the problem of a cylinder m a con- 
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duct iv e half-space is to be examined, the general integral representations 
must be considered. 


2-1 Derivation of Integral Expressions for E y and Hy 


We begin with Green's theorem (Stratton, 1941, p. 165) for any two 
scalar functions U. and at" , 


J 


Yol 


i U. VV - AT V\ | diCool) - J> j 


- AT 2 Z <4 (greet) 


(2-1) 


In equation (2-1), i- s the derivative in the direction of the outward 

normal to the volume of interest, which is defined in Fig. 1 as the volume 
between the surfaces and C^. 

If we add and subtract k M.V on. the left-hand side of equation (2-1), 
we obtain 

ju^+k^AT- j - axXl Z i (dre a). (2-2) 

fell ^ J 

Assuming that the scatterer is two-dimensional, then an integral 
expression for r W> the total electric field intensity in the source 
free volume between and C„, can be obtained by setting u. = E (j5) and 


at - G(p, j»') where E^jj) and G(p,|5') satisfy the equations 


-S(jr-p') 


(2-3) 


The vector p' refers to a source point on the surface and p refers to a point 
in the volume between and on the surfaces and C£. Applying equation (2-2) 
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to (2-3), we obtain 


J 


V.i 


ar«.d ^ 


However, equation (2-4) gives Ey(p) at a source point in terms of the 
total fxeld around the boundary. To obtain Ey(p), we invoke the reciprocity 
relation of the observation point and the source point. As a result, we may 
interchange p and p 1 in equation (2-4) and replace G(j3 ’,p) by G(p,p ') pro- 
vided the Green's function is symmetric. 

Morse and Feshbach (1953, p. 808) show that G(p,p') is symmetric if 
it satisfies some homogeneous boundary conditions on the boundary surfaces 
in both the p and p' coordinate system. Since the choice of which function 
we adopt for G(p,p') is optional provided equation(2-3) is satisfied (see 
Jackson, 1962, p. 18), we will ensure that G(p,p') is symmetrical by choosing 
the whole space value of G(p, j>'). Thus, a homogeneous Dirichlet boundary 
condition exists at infinity since G(p,p') satisfies the radiation condition. 

Equation (2-4) now becomes 


- - 


j ji') - & Cjs, p') $£ 3 . tpO 1 

c)n' J 


_ Gtf.j*') g> E g ( p ') 1 

'' M V' J 


(2-5) 


The total electric field intensity at any point is represented by the 


sum of an incident electric field ®y(f) and a scattered electric field 

E (p). 
y 
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£jif) * ejip + e“cp> 


(2-6) 


If we let the outer surface G ^ approach infinity, then the second integral 
of equation (2-5) becomes • 




f l is ' 


Ejn P ') p t^ F ') . s( fiP o^T(r) ? is'. 
^ ' • a a' J 


(2-7) 


Stratton (1941, p. 360) states that every two-dimensional electromag- 
netic field at great distances from the source can be represented by linear 
combinations of elementary wave functions of the form - x h (df^ (fep)jt * 

Assuming that the source of this wave is on a line, then the elementary 


wave functions are independant of 8 and only (jj^ is required. Thus, as 

<.k o 

p — *• QOj these (circular) cylindrical waves will behave like JL ' 

. .<!).. . 1 /7 

when H c Up] is expanded for large p. ' 

As C 2 approaches infinity, the scattered field E y C Ce’> and G(p,p‘) will 
behave as cylindrical waves for each point p ' on C 2 . Consequently, both 
functions will behave as except for constant factors. As a 

/f 

result, the second terra of (2-7) becomes 

r r ^ , iWp"' i ikp' 


A B 


k x l ff ^ 


ix r > i^)i 
'[7f> vf T,'{SjT J J 




'and the difference of the two terms under the second integral tends to zero 
for each p', and the contribution of the second integral is zero. 
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This is not true for the first integral since by assumption the 
incoming plane waves never vanish at infinity. In fact, the first integral 
is equal to - E y<P)> which can be seen by applying equation (2-2) to 

( V%- ) Ej (f )' =o 

(V 1 -*- \f)G (^p') ~ - S(p-j5') 

Equation (2 5) noW becomes 



] = {eiiptgtf) -F*( f ')|i( p , f -')jol s '. < 2 - 9 > 

By comparing equations (2-6) and (2-9), it is evident that the 
scattered electric field intensity is given by 



(2-10) 


• t* SC 

Similar expressions for (?) and Ky ( 
ft- = H^(j>) and AT = G(p,g') in equation (2-2' 


can be obtained by setting 
where 


(vN-kM tf) = d 
( V z t-^) G( pvpO - - ?') 


Following a development which parallels the above, we find that 



(2-12) 




G ( p , f -)^( p -0 . Us 



(2-13) 
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Equations (2-9) and (2-10) and equations (2-12) and (2-13) are 

integral expressions for and f ( P> respectively. However, since a 

t t / 

normal derivative of E^(0') and H^(£ ) is required rather than a tangential 
derivative, these expressions are not in a particularly convenient form. 
Andreasen (1965 b) has shown that by manipulating' Maxwell 's equations, the 
normal derivative of E and Hy can be transformed into more suitable 
forms. Using these representations, an integral equation solution for 
the unknown boundary values of E and H can be formulated easily. 


2-2 Suitable Expressions for the Normal Derivatives of E and H 

y y 


Starting with Maxwell 's first two equations and assuming a time 
dependence of , then 


VxE = - 5>S - 1 uu> H ' 

V x H <) P + J == - i tof \ £ 

& 5 


(2-14) 

(2-15) 


from which we may write 


iyttW - = C L 


(2-16a) 


ZyUu) H, 


= dE 


•oc 




- 


(2-16b) 


lyxu) H e = 


1 fy ^ > 


(2-16c) 
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( T- l \A> £ ) -rr C — £ Hy 


(c r-lu>i) 


(2-17b) 


L<r-lw0 Eg = £1^ - ih y H# . 

3/k 


(2-17c) 


In equations (2-16) and (2-17), we have assumed that the o 

<ttu 

dependence of all field components is of the form Jt 3J . Such an axial 
dependence would arise when a primary three-dimensional current distrib- 
ution and the field it radiates has "been expanded over a continuous 
mode distribution according to 


air J _ 


ilj. . 


(2-18) 


However, this problem will be discussed in more detail in Chapter 6, 
when three-dimentional source distributions are considered. It is suf- 
ficient to n'ote here that in the special case of a plane wave obliquely 
incident upon cylinders in free space, ky reduces to the axial component 
of the free space propagation constant. Whenever a plane wave is incident 
normal to cylinders in a whole-space, is zero. 

Substituting (2-17c) into (2-16a), then 


ljx\0 H. 


e>Wn _ c lu H , 


a) E f, j 

if 


(2-19) 
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Rearranging m terms of H , we find that 


« lL ^ _ Cr-twg) ^ „ 






(2-20) 


where 'y x ■=. V- 4 — Vy“ 


( 2 - 21 ) 


Similarly,^ substituting (2-17a) .into ‘ (2- 16c) 


y*.u>U € - ^ > 


(2-22) 


or 


^ * (.«■- iwe) + t-^y jjj 






(2-23) 


Defining the positive tangental direction on the contour by 


A A. / A 

5 ' n - *J > 


(2-24) 


and adding equations (2-20) and (2-23), -we find that the transverse 
magnetic field intensity, H , is given by 


= ^7 u + Cr-tiofO(v ft Ey)x^ 


(2-25) 


Following an analogous argument, we can obtain an expression for 

the transverse electric field intensity E. . Substituting (2-16c) into 

tr 

(2-17a) 


C <T- tw£)f 


A. 


tkj - t-W H. 




' CylXU) 




- Hx > 
^2 


(2-26) 



and substituting (2-16a) into (2- 17c) 


(<T- lut) E - 


Xf 



l / w 



E z - 

if 



(2-28) 


or 


^ t 


* yc 


u> 








>E V 

Ji 


(2-29) 


Adding equations (2-27) and (2-29) 



i/uo (? t Hj ) X j + i^V«, Ey - 


(2-30) 


Taking the dot product of § and equations (2-25) and (2-30), 


& N s = i-ty ^ 4- itoe) Ml. > 

<) n/ 


(2-31) 


and 


Eg - y itJ + tky 

3b. ' 5s 


If we multiply (2-31) by yub and rearrange, then 


(2-32) 



13 


or 


<3>e. 


34. = 


5 a' 


l + .gjgy. 0>Uj > 

Jt X ^ 


(2-34) 


where Z is the intrinsic impedance of the medium and is given by im 

k 

Similarly, equation (2-32) becomes 


Mm = - i* 1 F - ks M 


^ a' 


ek 


?F 17 


(2-35) 


If we substitute equation (2-34) and (2-35) into equations (2-9) 

and (2-12), we would obtain general integrodifferential representations 
t t 

for E y (e> and H y (e) in terms of the axial and tangential components of 
E and H on the contour. By applying the appropriate boundary conditions, 
i^e could determine E^, E g , and H g on the boundary and, with these 
quantities known, be able to calculate the scattered fields at any point 
m space. However, to conform with some of the terminology in the liter- 
ature, and to obtain a more physical interpretation in certain instances, 
it is desirable to introduce the concept of equivalent electric and magnetic 
surface current densities (see Harrington, 1961, p. 106) . 

In brief, the equivalence principle states that many source distrib- 
utions inside (outside) a given region can produce the same field 
outside (inside) that region*. As a result, equivalent electric and 
magnetic surface currents can be introduced to generate the same field 
external (internal) to a region as do the original sources internal 
(external) to that region. For our problem, these surface currents are 


* Note that this informs us that the inverse problem is not unique. 
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chosen such that the boundary conditions, continuity of tangential E and 
tangential II, are satisfied. Thus, the equivalent electric surface current 
density is set equal to the tangential component of the magnetic field 
intensity and the equivalent magnetic surface current density is set equal 
to the tangential component of the electric field intensity at the boundary. 

From the continuity of tangential H, we obtain 


A 

/rt 


or 


* UU- ) = o , 


(2-36) 


m. x H«*t - ^ * K > (2-37) 

where rl is taken to be the outward normal from the scatterer, We conclude 
then, that the equivalent electric curface current density components are 
given by, 

=■ Ky i (2-38a) 

A xH^=> * - K s * (2-38b) 

Similarly, the equivalent magnetic surface current density components are 
g iven by, 

(2-39a) 
(2- 3 9b) 

Upon introducing equations (2-38) and (2-39) into equations 
(2-34) and (2-35), our final expressions for the normal derivative of 

V* 

Ey and on the boundary become 


£" s - ' My j 

* - M s '* 


-$4;l 

Z a 


k, 

X J 


r 


(2-40) 


l) Uy 

bn' 


_ l y" M», + Jv C>M S 


(2-41) 
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2-3 General Field Representations for the Exterior Region 


Coupled integrodif ferential representations for E^(p) and H^(p) 
are obtained when equations (2-40) and (2-41) are introduced into (2-9) 
and (2-12). 




and 


4- f M-ip') J 

-'Ci chi/ 

Jtp)* if) - [ 6(p,p') l f'> I ^ 

Jc. ( ?k 2k J 


(2-42) 


+ f MjpO 2£ ( f,f'> • 


(2-43) 


It is shown in Appendix A that the two-dimensional Green's function 


is 


«u,r> - i »o 't-xif-f'i) , 

4 


(2-44) 


where H 0 is the Hankel function of the first kind and order zero, 
representing an outgoing wave, and X is the transverse propagation 
constant, and is related to the propagation constant A of the medium 
by equation (2-21). 

The derivatives of the current components in equations (2-42) and 


(2-43) can be eliminated analytically with integration by parts. On 
carrying out the partial integration, and introducing equation (2-44), we 
find 
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Ej If) = ejlp - st C kj(f') 

4 k 4 :. 


+ kj , j^jjo on^pt^ 3 ' 4-.L f M s y)^(y)|?-fj) ^ } 
4 k v. r ^S' ^ 4» e^a* 


(jO * (jf^ 4 X- f H CpO W c Wp-p' 0 ^ s 
-lit* f M s (p') rf (yi ? _ ? '<)Js / +± ( K/pO ^(Vlp-p'i) is'. 

lit Jc W 4 I W r r 

The normal and tangential derivatives of H p l^-p')) are 
given by 

3§y^ = • 


Since 


aZ a (-.*0 - - -i -h ,(3) j 

<* 3 3 

where (j^ is any Bessel function, 


and 


then 


2 t 3 ) * C-PS.ij) , 


azu - *) » & Z»0o . aoo » y[-AZ-^) ♦ Z 

7^ TT^) a* l* 3 - J«u 

Thus, the transverse derivative of is given by 




at Jk + a t 

dx 1 eU' . 


(2-45) 


(2-46) 


(2-47) 


(2-48) 


(2-49) 


(? /U-*') 1 * c?-r) r ) 



17 


= Y r U-V)-S * tf?(YI ? -jr'i) 


If-f ' 1 


Y C Ctrl j*> £ + s'^w ji 2 ) 


(2-50) 


where is defined in Fig. 2 as the angle formed by the $ axis and the 
line between the point of observation and a point on the contour. 

From Fig. 2, we have 

ft/ a e* & +• SUvvoC Z (2-51) 

i * 

/n = Get (cc+ f6o w ) A 4- 5 Lw (_<st_ +■ 1 8 o" ) g - - ft ^ (2-52) 

5 * (.“£.+ fo°) a. + SLW C«t+1»*')2 4 60 S* 2. (2-53) 


Using equations (2-50) , (2-51) and (2-53), the normal and tangential 
derivatives of H^(y)^-y’l) become 


i Ho'LVIj-p'l) 

1 - y [ftrf fl 4- S Wv |S s£v\ «c ) C# S F- /( ) 




» Y 1 a-*) rfHVIp-^0 i 

(2-54) 

and 

d (.Yl/i-S'l, 
> 1 

) s Y" 5 i/H (jk-oi.) (, Y ) * 

(2-55) 


hf 

When equations (2-54) 'and (2-55) are introduced into equations 


(2-45) and (2-46), we find 
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E j(f) " ifU/p-p' i)d*' 

+ Ljjvjf ( Ks(?0 svKCp— t) H?Ulp'p'l)^' + f M s (f') 0rtlp->x) l-|f/&ip-f'l)c£ s J 

"Ik 4 ^ 4, (2-56) 


aLiu 

h‘ ( F ) = Hj ifU _ll P M if) «*> Ul f-f-'l)^' 

42k J c 

i 

- M f Mi if') s~<.p-0 Hf?lVlf-f'l)ol«' +il f Mj'Wf--' 


42 k 


(2-57) 


Equations (2-56) and (2-57) are the desired coupled integral rep- 
t ^ t * 

resentations of E /?> and Hy(j>) for the exterior problem. The point of 
observation, 5 , is exterior to the inhomogeneities present, and 
represents a contour integration around each scatterer present. Once the 
equivalent surface current densities are known, the field components 
are determined from equations (2-14), (2-15), (2-56) and (2-57). 

From equations (2-14) and (2-56), the scattered transverse magnetic 
field intensity is given by 


LixlO 


( e j 3) * _L(25* - 2 ) 

yXLJ V ** )sK J 


(2-58) 


It follows from equation (2-49) that 
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Similarly, 
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and 







(2-63) 


(2-64) 


When equations (2-56), (2-59), (2-60), (2-61), (2-62), (2-63), and 
(2-64) are introduced into equation (2-58), we find that the scattered 
transverse magnetic field intensity is given by 
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(2-66) 


Similarly, from equations (2-15) and (2-57), the scattered trans- 
verse electric field intensity is given by 


•sc 


tr* 


J .Vx(h;5)= -ii ixg'*. *u/ c i) 
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(2-67) 


It follows from above that 


E~ = it? J Hyp llftflf.f'l) (,^(SA - cos j3 z 
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+ W? Ctflp-?'0 ( z ) I ds' 

V i p^V' i 

+ ^ ^ r K s (j5') y ct-^c H 5 o i - cos p 2 ^ 


- H^C^lp- p'O ( svwLa,^.-c<)^ - c°s (v.ji> -<*)£)(. <ds / - 
^ I j5 - jf' ! J 


2-4 General Field Representations for the Interior Region 


We wish to obtain general integral representations of E, yW and 
H*” (p) for the interior problem. As in the exterior problem we begin with 

y 

the modified form of Green's theorem, equation (2-2), and set the two 
scalar functions U. and /vr equal to E^Q) and G(p,^') respectively. Since 
our volume of interest is bounded only by an outer surface C , Green's 
theorem yields the result that 


£-; ip - - 


— 'C, L <)a J 


(2-69) 


In equation (2-69) } rv is the outward normal to the volume of interest 
(Fig. 3) and is equal to -a* of the exterior problem. G*" (|>, g') is the whole- 
space Green's function of the interior region, and is given by 




(2-70) 
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Thus, our integral expression for Ey(^) is 


-j *- "J 7 {^'F^ c ? /)_ b j ( P^ l N 


£?(?)« - 


p. fl') ( ots*' . 


(2-71) 


By comparing equations (2-71) and (2-10), it is evident that our 

t _ 

integral expression for for the interior problem is equal to the 

negative of our integral .expression for E^({5) for the exterior problem. 
Thus, we may write the integral representation of (§> from equation 
((2-56) as 


’’ 4* t Jc, 

(2-72) 


- iV; 




Since similar conclusions can be drawn for H (g), we have from 
equation (2-57), 

' 4*t\4r ' 


- th. ' c ° 3 v t ~ u) ! f • - 


(2-73) 


From equations (2-66) and (2-68), the transverse field components 


are given by 
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(2-74) 


(2-75) 
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2-5 Representations for E and H When the Scatterer Has a Low Surface 
Impedance 


The integral representations for E and H simplify considerably 
when the scatterer is assumed to be highly conducting . The reason for this 
is that E and H take a particularly simple form just inside the scatterer, 
and may be assumed to be zero further from the boundary as a result of 

attenuation. 

I 

The angle of refraction at any local site on the contour is given by 
Snell's Law as 


CoS 





(2-76) 


In equation (2-76) , 0-j^ is defined from Fig. 4 as the angle of incidence 
and 0 2 as the angle of refraction. It is evident that if , 0^ 

is approximately zero even far grazing exterior angles of incidence. As 
a result, the total field at any point just inside the scatterer may be 
represented by a wave propagating normally away from the coutour. Following 
Stratton (1941, p. 354), we may write that the transverse electric portion 
of the field must satisfy the relationship 


i-s ~ ~ U t , 


3 


(2-77a) 


(E^ - E n = = 0) and the transverse magnetic portion of the field 

must satisfy the relationship 


e j =2 ‘j Hs > 


(2-77b) 
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(E =E =H =H =0). Z. and Z. are the cylindrical surface impedances 
s » y n iy 

which Stratton (p.360) shows are given hy 

^ h (■ 5) 

^ y (2- 78a) 

AMT $,<*> 


,/r* i>' ^(s) 
A ~3j- 


(2-78b) 


represents a radially outward travelling wave and is given by 

{/ki-kt f ) . (2 79) 

Since the wave is normal to the interface^ A. is zero in equations (2-78) and 


(2-79) and we find that 


*- 5 * f-x - C M 

LA(. h<? tfef) 


(2- 80a) 
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When 1 a is Large, we can expand H"( ) as 
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(2-80b) 


(2-81) 
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so that 


(0 


jClk.* 

«r (Kf> 


-i Jt 
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- - £. 


(2-82) 


where j kip / >>/ ^ j^pl» (\ and ~~ ~ Af j ) £ "A. . 

Equation (2-82) is a good approximation if the curvature is small compared 
to the attenuation per unit distance. In this case, equations (2-80) 
reduce to plane wave surface impedances given by 


and 



(2-83a) 



Ay 


i 


(2-83b) 


where tbe subscripts on Z are taken to agree with the subscripts on E 
to indicate that the plane wave impedance is calculated using the propa- 
gation constant of that component of the electric field. 

As well as being a very simple representation for the internal field, 
equations (2-77) are important in that they relate the equivalent magnetic 
surface current densities to the equivalent electric surface current 
densities. When equations (2-38) and (2-39) are introduced into equations 
(2-77), we find that 




(2-84a) 


M 


- z ijK r 


(2~84b) 


As a result, our field representations for the exterior region, equations 
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(2-56), (2-57), (2-66) and (2-68), reduce to 

Ejl /0 - 2** f Ky Ip') CVJ^^l)is y 

4 k Jc 

4k J c, 4 -4 ^ ’ ' 

' (2-85) 
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ip yp )oU' + lV f K s (p') C6 5 ( f--0 Yl p-p'i) L'. 

1 Jc, 


(2-86) 


The scattered transverse field quantities become 

Vr - U s f k„tf) H®Cns.5'i)6i»^ U>' 
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(2-87) 
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CHAPTER 3 


SCATTERING FROM PERFECTLY CONDUCTING CYLINDERS 


We wish to evaluate scattering from cylinders for the limiting case, 

Z = Z. =0. Further, we will assume that the field is constant in the 

i„ i 3 

y 

y direction (ky - 0) as is the case, for example, when a plane wave is 
incident normal to the cylinder axis. By making these approximations, 
equations (2-85) and .(2-86) uncouple and reduce to very simple integral 
representations. As a result, the solution of the integral equations 
encountered in solving for scattering from cylinders of arbitrary impedance 
is more easily understood, and in addition, many of the numerical problems 
which are met may be studied individually. 

Mei and Van Bladel (196.3b) and Andreasen (1964) have examined the problem 
of scattering from perfectly conducting cylinders. Andreasen (1965a) also 
has examined the problem of scattering from perfectly conducting bodies of 
revolution. In this chapter,, the principle difference from their studies 
will be in the physical situation to which the solution is applied. 


3-1 Basic Integral Equations 


If we assume that Z,- = Z. =0 and that k. = 0 m equations (2-85) 

^-s iy . .y 

and (2-86), then we find that the most basic Integral representations 
for E and H are 

y y 



(3-1) 
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and 








(3-2) 


or 


and 



(3-3) 



A f ks(f') U? LMf-f'i) , 

4 Jc i 


(3-4) 


where we have referred to the earth parameters by the subscript 1. 

The electric field intensity radiated by a line source may be 
determined by setting the surface current density in equation (3-3) equal 
to a delta function source. With 

Kj tjr'J - (3 ' 5) 

where I is the current in amps on the line source, 
equation (3-3) becomes 


e D i f ) = ' I (3-6) 

i 1 

From equation (3-6) it is evident that the integral representation for the 
scattered field given by equation (3-3) may be interpreted as the summation 
of a continuous distribution of weighted line sources around the cylinder 
contour, radiating into the exterior region. In fact, it is shown m 
Appendix B that equation (3-3) may be derived from just such a physical 
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argument, rather than the mathematical approach used in Chapter 2, A similar 
use of physical reasoning could have been used to derive equation (3-4) . 

Before calculating the scattered fields from equations (3-3) and (3-4), 
we must determine first the current distribution about the cylinder. To 
do this, we enforce the boundary conditions on tangential E and H. For 
perfect conductors, these boundary conditions are 

% (3-7a) 

= F . (3- 7b) 

If we let the point of observation approach -a boundary point given 
_ // 

by the vector p and apply equation (3- 7a), then equation (3-1) reduces to 


Bj (yJ'O *. jj, t/> f Kyts')U 0 ' (3-8) 

1 JC, 

Equation (3-8) is a singular Fredholm integral equation of the first kind 

and can be solved numerically for the unknown function K^QO once the 

contour has been specified. If we assume that the incident plane wave 

possesses only an axial component of electric field intensity (hereafter 

referred to as E -polarization), then we may write that E^p*) in equation 
y y * 

(3-8) is given by 


t^") ~ x. 


(3-9) 


where H q is the magnitude of the magnetic field intensity in air, and ^ 
is the direction of propagation of the incident wave in Fig. 2. 

Similarly, when the boundary condition given by equation (3- 7b) is 
applied to equation (3-2), we find that 
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- - K*( f") - jA f (3-10) 

A \ 

Equation (3-10) xs a singular Fredholm integral equation of the second kind. 
If we assume this time that the incident field possesses only an axial 
component of magnetic field intensity (hereafter referred to as H^- 
polarization), then we may write that Hy(^") in- equation (3-10) is given 

by • , „x ,, tk + Z" ) 

If ' ° H 0 x- * (3-11) 

Once equations (3-8) and (3^10) have been solved for Ky(|T) and 

K s (p0> we can calculate the scattered axial components of the field 

from equations (3-3) and (3-4) . The scattered transverse components 

become , from equations (2-87) and (2-88) 



1^-jr'lK - ariLip-*) 2 

Jjf-f'l 



312 Numerical Solution of Singular Fredholm Integral Equations 


(3-13) 


Kopal (1955, Chapter VIII) discusses several numerical methods used 
to solve Fredholm integral equations . However, the common approaches 
used in solving integral equations which occur m scattering problems 
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may be divided into two general philosophies. These are referred to as an 
expansion of the unknown 

1) In functions of full range support, and 

2) In functions of subrange support. 

In full range support, one function is chosen to represent the unknown 
around the contour, whereas in subrange support the unknown is expanded 
in functions over small intervals. In investigating geophysical scattering 
problems, an expansion of the unknown m functions of subrange support is 
preferred to an expansion of the unknown in functions of full range support 
since 

1) The matrix is better conditioned 

2) The integration time is faster, and 

3) The edges of rectangular inhomogeneities can be handled 
better. 

Having decided on expansion of the unknown in functions of subrange 
support, there are a variety of methods available for interpolation in the 
use of these expansions. A suitable method is to expand the unknown in a 
set of N algebraic functions and require the integral equation to be satis- 
fied at N points. We have not investigated another popular method, that 
of sinusoidal interpolation, since we do not have a good estimate for the 
wave number in sinusoidal interpolation. 

We will discuss two sets of algebraic functions which can be used. 

The first set of functions has the advantage that a solution can be obtained 
clearly and easily, whereas the second set has the advantage of being a 
more accurate representation of the unknown. 



3-2-1 Constant Current Density Approximation 


We assume m all discussions that the contour of the cylinder has been 
divided into N straight intervals as is shown in Fig. 5. Then we can choose 
to represent our unknown current density by a set of algebraic functions 
consisting of a constant current density over each interval. This is 
represented graphically in Fig. 6, and can be written as 


Kl 

K if") = 2! K-UjCs') 


(3-14) 


where 




a 


is 


the 


t Isevtf Le.r<c 
contour coordinate. 


'-,W 


s is the midpoint of the jth interval, 

w^ is the half-width of the jth interval 

✓ 

and K is the amplitude of the electric or magnetic surface 
current density in the jth interval 


When equation (3-14) is introduced into (3-8), we find that 


Ho 


i 1?, Cfc* fa + Z " c o! - fa ) 


M 


p 1 K, au'l tftMp-f'l)" 12 ' 

j-l Jvi J 1 * 


tO 


ah 1 Kj K ' 
4 i-w- 


(3-15) 

The resulting integral can be computed numerically over each interval 
of the contour. Thus, equation (3-15) 'reduces to a linear equation with 


N unknowns . 
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If we enforce the' boundary condition at the midpoint of each interval, 
we will obtain N linearly independent equations with N unknowns, and the 
problem has reduced to solving a matrix of the form 

(^XKj) = (/}) . (3-17) 

Using standard matrix inversion techniques, N sampled values of the current 
distribution are obtained, 

(fcy) * (B) (fl) (3-18) 

Equation (3-18) is an approximate solution of equation (3-8) . 

3-2-2 Quadratic Current Density Approximation 


Unfortunately, the constant current density approximation given by 
equation (3-14) is not always accurate, especially over those intervals 
on which K(s ') varies rapidly. Rather than take a smaller sampling interval 
in these regions, a better approximation is to choose a set of algebraic 
functions which consists of a quadratic function over each interval. Thus, 
we may write 


K C$0 =r Ur Cs') 
J 


(3-19) 


where 


Uj is') = j rtj + ts'-Sj) + CjU'-Sj)*- > U r v d )<s'^ ( ywj) i 

Equation (3-19) contains three constants per interval which are determined 


by invoking three boundary conditions per interval: 
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1) Enforce the usual E.M. boundary conditions at the mid- 
point of each interval ^ 

t 

2) Enforce continuity of the quadratic function at the edges 
of each interval, and 

3) Enforce continuity of the first derivative of the quad- 
ratic function at the edges of each interval. 

However , this would increase the matrix size by a factor of three, and 
hence the solution time by a factor of 27(3 ), although superior accuracy 
would be achieved. 

An approximate solution which yields an N x H matrix is that which 
assumes that the only unknowns are the set K. where 

j 

K. ~ K.(Sj) — amplitude of the current density at the mid- 
point of each interval. 

With this set of unknowns, we interpolate over each interval by fitting 
a quadratic function to the amplitude of the current density at the mid- 
point of the jth interval and the midpoint of the two adjacent intervals, 
as is shown m Fig. 7. While the accuracy of , this set of quadratic functions 
is not as great as that for the original set chosen, it is still greater 
than the accuracy of the constant current density representation, and yet 
it requires essentially the same computer time as the constant current 
density representation. 

" 1 

If we fit equation (3-19) to the amplitude, of the current density 
at the midpoint of the jth interval and the midpoint of the two adjacent 
intervals, we find that 
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K o . A. 


k- = A - Bj (w t- w ) t C ( w. +w, V 
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K i*. * A J + YW' * YVV' . 


(3-20) 


By allowing the half width of each interval to be distinct in equations 
(3-20), we have allowed ourselves the freedom to choose smaller intervals 
for those regions where the current density varies most rapidly. 

Solving equations (3-20), we find that 

"i- k j 

B o ■ K.UI - Kj-Aj K,„ 


(n-A J Xw j + 'V 1 ) 
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where 
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Introducing these values for the coefficients into equation (3-19), equation 
(3-8) can be rewritten as 
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As before, the remaining integrals in equation (3-22) can be performed 
numerically, and thus the integral over each interval contributes to three 
coefficients in each row of the coefficient matrix. It should be noted, 
however, that care must be taken in programming equation (3-22) to ensure 
that the Nth and first intervals are coupled whenever the contour C-^ is 
closed. 

It is possible to rewrite equation (3-22) in terms of alone ^ 
and we find that 

$>; rt + 
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(3-23) 


However, equation (3-22) is preferred for its greater ease in programming. 

Equation (3-22) reduces to a linear equation with N unknowns, similar 
to equation (3-16). By enforcing the boundary condition at the midpoint of 
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each interval, we obtain a modified coefficient matrix which yields an 
improved estimate of the N sampled values of the current distribution. 



f B T PJ 


(3-24) 


A similar interpolation procedure can be applied to equation (3-10), 
assuming H - polarization, and we would find that 
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(3-25) 


where K g ^ is the amplitude of the current density at the point ^ where the 
boundary condition is being enforced. 

Once equations (3-22) and (3-25) have been solved for the N sampled 
values of Ky and Kg, the scattered field quantities can be calculated from 
equations (3-3), (3-4), (3-12), and (3-13). 


3-2-3 Integration Through the Point of Singularity 

When the integration is carried out over the interval m which the 
boundary condition has been applied (the ith interval), the Hankel function 
possesses a singularity and the contribution of this interval must be 
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evaluated analytically. Fortunately, the argument of the Hankel function 
is small in this interval, so that we can make the following approximations: 
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(3-26) 


and 
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(3-27) 


where )f / = In )(” «= .57722 = Euler's constant, 

s 1 is the contour coordinate measured from the center of the 
ith interval (see Fig. 8), 

and $ is the distance of the point of observation above the 
center of the ith interval. 

If tklfp.p'i 6 - 5 ) an error of less than 0.1% in the value of ri^(k|^-pM) 
and lij\k|^-£'l ) is obtained when using equations (3-26) and (3-27). 

Approximating the Hankel functions by equations (3-26) and (3-27), the 
contribution of the ith interval to HT and H is evaluated with the point of 
observation a distance S above the contour. Having obtained these 
integrals, then the contribution of the singularity to the boundary con- 
dition is taken as the limit as £ approaches zero. It should be pointed 
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out that since the original boundary conditions on E and H. (equations 
(3-7)} are derived by a similar limiting process, this procedure is completely 
rigorous . 

If the point of observation is taken to be a distance § above the 
center of the ith interval, then the contribution of the singular interval 
to equation (3-22) can be written as 
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(3-29b) 


(3- 29c) 


Equations (3-29) have been evaluated. in Appendix D assuming that 
H‘:’U m s' 1 " ) is given by equation (3-26) . After the integrals have 

been obtained, we find that 
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(3- 30c) 


When equations (3-30) are introduced into equation (3-28), we find 
that the contribution of the singular interval to the coefficient of 
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and the contribution to the coefficient of K - is 

Yi-1 


= K^Ul * (3-33) 

Equations (3-31), (3-32) and (3-33) are to be substituted for the 
right-hand side of equation (3-22) whenever j = i. 

A similar limiting process is required to evaluate the contribution 
of the singular interval in equation (3-25), We have from Fig. 8 that 
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(3-34) 
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It follows from above that the contribution of the ith interval to equation 
(3-25) can be written as 




(3-35) 


(3 -36a) 


(3-36b) 



Equations (3-36) have been evaluated m Appendix D assuming that 

(kj/s"*+ s' x ) is 8^ veTl fey equation (3-27). After the integrals 

have been obtained, we find that 
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When, equations (3-37) have been introduced into equation (3-35) , 
then the contribution of the singular interval reduces to 


Kr 


which must be added to -K in equation (3-35). 

s i 

integral equation for Hy polarization as 


(3-38) 

Thus, we may rewrite the 
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3-3 Symmetry Considerations 


Andreasen (1964) has pointed out that a considerable amount of 
computer time can be saved whenever the cross section of the scatterer 
is symmetric with respect to the z axis. In’ this case, it is possible 
to expand the incident field into a sum of an even mode of x and an odd 
mode of x and solve separately the matrix equations for the even and odd 
modes of the equivalent current density. 
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For example^ assume Ey-polarization and expand the incident field 


as 


L — * ^ 

e j ■ v f 6 * . 

and the propagation constant as 

A = ^ ** 4 . ^1 i 

Then it is evident from equation (3-9) that 
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(3 -42 a) 


(3 -42b) 


When the even and odd mode of the equivalent current density are separa- 
ted, the coefficient matrix for each mode is given by 
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(Bi -) = 2 (B- + B/ . ) 


L = /. N/x 
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(3 -43b) 


where we have assumed that the first interval is the mirror image of the 
Nth interval with respect to the z axis . 

The even and odd modes of the equivalent current density are deter- 
mined by solving separately the matrix equation obtained when equation 
(3 -43a) is equated to equation (3-42a) and equation (3 -43b) is equated to 
equation (3-42b) . Thus, the total equivalent current distribution is given 


by 
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Since matrix inversion routines are proportional to K , it has been 
four times faster to invert two matrices of orderN/2 than it would have 
been to invert just one matrix of order N, 

Whenever the scatterer is symmetric with respect to the direction 
of propagation, only an even mode of the equivalent current exists, and 
assuming vertical incidence, it is given by 
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(3-45) 


The even coefficient matrix is given still by equation (3-43a) and, 
having solved the matrix equation obtained by equating equation (3-43a) 
and (3-45), the total equivalent current distribution is given by 
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3-4 Numerical Examples 


If the integrands of equations (3-3), (3-4), (3-8) and (3-10) are 
smoothly varying functions along each interval, then it is sufficient to 
approximate the integrand m this interval by a parabola. Thus, we may 
use Simpson's rule with n = 2 to integrate numerically across the interval. 

For example, the approximate value of the^integral across the jth interval 
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of equation (3-22) is 
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where the positron vectors pj_ and have been defined in Fig. 9. 

To demonstrate the validity of the integral representations and their 
numerical solution, we will compare the numerical results with the analytical 
results obtained for the case of scattering from circular cylinders. It 
is shown in Appendix C that the transverse magnetic field intensity which 
is scattered by a perfectly conducting circular cylinder in the presence 
of an Ey-polarized plane wave is 
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where 8^ * 1 | } ^ - <? 

( 2 - m i\ 

and R is the radius of the cylinder. 

The induced surface current density around ,the cylinder is 


I I 

>,n). an. % 1 Z r*g. 

irfe.R hSU-JO 

and 

M S C§) = o. 


(3 -49a) 


(3 -49b) 


The distance h m equatxons (3-48) and (3-49) is the depth from the earth- 
air interface to the center of the cylinder. Since some point must be 
adopted as origin for a phase and attenuation reference, the earth-air 
interface has been chosen to facilitate comparison of these results with 
those of the half-space problem considered in Chapter 7. In addition, 
we will chose our incident field to be.- the transmitted field of a plane 
wave incident upon a conductive half-space. This choice of incident field 
also will facilitate our 1 ":dis cuss ion in -Chapter 7 when we wish to estimate 
the significance of coupling between the conductor and the earth-air 
interface. It should be noted, however, that this still is a "ficticious” ' 
model in that reflections from the earth-air interface are ignored when 
computing the currents induced on the surface of the conductor and 
also, the earth-air interface is ignored when computing the scattered fields 
at this boundary. 

We have from Stratton (1941, p.,493) that the transmitted electric 
field intensity of an Ey-polarized plane wave incident upon a conductive 
half-space is 
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where ^ is given from Snell's Law, equation (2-76) as 
Gtr$ * / 1 - * * 


(3-50) 


We will consider only the special case of normal incidence since in 

» 4 

most problems of geophysical interest the transmitted field propagates 
normally away from the interface even for grazing angles of incidence. 

In addition, this approximation yields the' greatest saving in computer 
time as a result of the symmetry. Thus, equation (3-50) reduces to 
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(3-51) 


/ 

From equation (3-51), it is evident that the desired incident field for 
Ey-polarization is obtained by setting equal to .180° and replacing 
z , H 0 by 


az, h* 


(3-52) 


in the pertinent equations. 

The unknown equivalent electric surface current density around the 

cylinder was determined from equations (3-22) after Z,H had been replaced 

i o 

by (3-52) . The cylinder coutour was approximated by inscribing 10, 20 and 
40 intervals as is shown in Fig. 10 for N = 20. By increasing the number 
N of sampled values, we estimate more accurately the cylinder contour. 
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the integration around the contour, and the current density distribution. 

Figs. 11 through 16 demonstrate that the integral representations 
derived for Ey scattering from perfect conductors are valid, and that 
the numerical solution rapidly converges as N increases. We have assumed 
normal incidence m this example and that: the depth to the top of the 

cylinder is 20 m, the cylinder radius is 100 m, the incident field frequency 

' -3 

is 1000 hz, and the conductivity of the whole space is 10 mhos/m. 

Since the induced current density is symmetric with respect to the 
z-axis, only part of the results have been plotted m Figs. 11 and 12. 

It is seen that the maximum current amplitude occurs on the shadow side 
of the cylinder, as discussed by King and Wu (1959) for the case of 
scattering from circular cylinders in air. However, since the cylinder 
is situated in a conducting medium, the minimum current amplitude is 
due also to attenuation trf the incident field as it propagates through 
the ground. The maximum errors occur at the points of maximum and 
minimum illumination, although the errors rapidly decrease to less than 
VL as N is increased. 

In Figs. 13 through 16, the scattered magnetic field intensity 
obtained from these N sampled values of the current density has been 
plotted for points of observation near the cylinder on the plane 20 m 
above it. It is seen that a slight instability in the field occurs 
close to the cylinder, especially for N = 10, but disappears as N increases. 
This problem, which will be discussed below, is caused by the inaccuracy 
of a parabolic approximation to the integration along that part of the 
contour in close proximity to the observer. However, it is sufficient 
to observe here that ,the instability obviously decreases as the sampling 
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width decreases. We see that for N = 40, the maximum error is less 
than 0.5%. 

It can be shown by following a development similar to that of 
Harrington (1961, p. 235) or Appencix C. that the axial magnetic field 
intensity which is scattered by a perfectly conducting circular cylinder 
m the presence of an Hy-polarized plane wave is 
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As m equations (3-48) and (3-49), we have chosen our origin to be at the 
earth-air interface. 

It follows from equation (3-51) that the transmitted magnetic field 
intensity of an Hy-polarized plane wave incident normal to a conductive 
half- space is 

ikz" 

(3-54) 
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Thus, the desired incident field for Hy-polarization is obtained by 
setting equal to zero and replacing H Q by 

£ H« , 

, t + 

in the pertinent equations . 


(3-55) 
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10, 20 and 40 sampled values of the equivalent electric surface 

current density were determined from equation (3-39) after H had been 

o 

replaced by (3-55), considering the same cylinder studied for Ey-polar- 
ization. Figs. 17 and 18 show the axial (y) magnetic field intensity 
obtained from these sampled values for points near the cylinder. Again, 
an instability is observed, and although the amplitude of the oscillations 
is larger than for Ey- polarization, it also obviously decreases as the 
sampling width decreases. Since these oscillations are meaningless with 
respect to the actual problem, no attempt has been made to define them 
properly. We see, however, that good convergence of the Hy-polarization 
integral representations is obtained for N = 40, with the maximum error 
being less than 1.0%. 

We have already pointed out that the instabilities arise from the 
inaccuracy of a parabolic approximation to the integration along that 
part of the coutour in close proximity to the- observer . This can be 
seen in the follbwing way: Assume that the ’equivalent surface current 

density, across the jth interval is approximately constant. Then the 
contribuiton of the integral along the jth interval -to the scattered 
field in equations (3-3) and (3-4) can be written as 
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Since the Hankel function is singular when both the real and imaginary 


parts of the argument, are zero, it is evident that the value- of the 
argument of equations (3-56) and (3-57) will pass through a pseudo- 
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singularity for points of observation close to the contour. Fig. 19 
represents a hypothetical plot of the integrand of equation (3-56) or 
(3-57) for a point close to the contour, and the "best" parabola which 
would be interpreted for this curve. It is seen that the area under the 
parabola is not a good approximation to the true area. 

If the point of observation is to one side of the interval, as in 
Fig. 20 , then a parabolic approximation to the integration of equations 
(3-56) or (3-57) is much better. Although the parabola still does not 
represent a good approximation to the integrand, it is seen that the 
area under the parabola could constitute a fair estimate of the area under 
the hypothetical curve since the positive and negative errors tend to 
compensate one another. 

Inasmuch as this problem arises when the argument of the Hankel 
function is small, it can be overcome by making a small argument expan- 
sion and integrating analytically across these intervals as was done in 
section (3-2-3). For the best accuracy, this should be done for points 
of observation just to one side of the interval, as in Fig. 20, as well 
as for points of observation over the interval. Thus, small argument 
expansions will be used for intervals near where the boundary condition 
is being enforced, as well as when calculating the scattered fields. 

Before the small argument expansion is made, it is convenient to 
rewrite the integral representation in a more suitable form involving as 
few analytical integrations as possible. It is shown in Appendix D 
thatj in general, the contribution of the jth interval to the boundary 
electric field intensity (equation (3-22)) can be written as 
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The geometrical symbols in equation (3-58) are defined in Fig. 21. 

<2 - 1 2 1 

Using the small argument approximation for H + s' ) given by 

equation (3-26) , equations (3-60)' have been evaluated m Appendix D. 

Similarly, it is shown in Appendix D that in general the contribution 
of the jth interval to the scattered magnetic field intensity (equation 
(3-12) can be written as 
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(3-63c) 
(3 -63d) 


Equations (3-63) have been evaluated in Appendix D‘ after approximating 
Hf (kj7 £ + s' 2 ) by equation (3-27). If equations (3-58) through (3-63) 
are used to calculate the contribution of the' jth interval when both 
l k l(l^ - and IkjJlp - •are less than 0y3, then the observed in- 
stabilities disappear and an accuracy comparable to that obtained on the 
flanks of the scattered field results for points of observation near the 
scatterer contour. 


A similar analysis must be carried out for H^-polarization. Again, 
it is shown in Appendix D that in general- the contribution of the jth 
interval to the boundary axial magnetic field intensity (equation 3-39) 
can be written as 


H . ik.S -B 


(3-64) 


where B is given by equation (3-62b) . The approximate contribution of the 
jth interval to the scattered axial magnetic field intensity is given 
by equation (3-64), but with opposite sign. 

To demonstrate that equations (3-58) through (3-64) are necessary 
to predict accurately ’the scattered field for points of observation close 
to the scatterer without taking an excessive number of sampled values of 
the surface current density, we have chosen a particularly good example. 
The real and the imaginary parts of the axial magnetic field intensity 
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have been plotted in Figs. 22 and 23 assuming an Hy-polanzed plane wave 

incident normal to a conductive half-space. The depth z-^ to the top of 

the cylinder is 3 jn, the cylinder radius is 100 m, the incident field 

-3 

frequency is 1000 hz, the conductivity of the whole-space is 10 mhos/m, 
and the cylinder contour has been approximated by 20 intervals . It is 
seen that the undesirable oscillations have been eliminated completely 
by using the small argument approximation given by equations (3-64), with 
a maximum error of about 1% being obtained through 20 sampled values of 
the current density. 

The important conclusion to be drawn from Figs. 11 through 18 is 

that if the numerical integration is inaccurate, a significant change m 

all the sampling widths brings about a significant change in the predicted 

scattered field. We are led to believe that if the predicted fields then 

converge, the numerical integration and sampling have been performed 

correctly. In future discussions, a small argument approximation will 

be used whenever (|k^|jp - g'|) — .3 for any interval. 

To illustrate the generality of the program, we will examine the 

field scattered by a vertical slab which is 30 m wide and 300 m deep, 

assuming an Ey-polarized incident plane wave. The general cross section has 

been drawn in Fig. 24 (a) and a detail of the contour at one corner has 

been drawn in Fig. 24 b. A sharper corner for the slab could have been 

chosen, but this would require a tighter sampling interval to define 

properly both the contour and the current density. We have assumed that 

the depth z to the top of the slab is 20 m, the incident field frequency 

-3 

is 1000 hz and the conductivity of the whole-space is 10 mhos/m'.. 

The variation in the magnitude of the equivalent electric surface 
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current density from the top (s 1 = 0 m) to the bottom ( s ' = 312 m) of the 
slab is shown in Fig. 25. The variation of the phase of the equivalent 
electric surface current density around this part of the contour is shown 
in Fig. 26. Since only slight differences are obtained in the current 
density distribution for 30 and 42 sampled values, we have assumed that 
a convergent solution has been found. 

A very interesting variation m the magnitude of the surface current 
density is seen to occur around the bottom limit of the slab. We see 
that the current density increases around the bend m the contour and 
then is constant on the bottom flat portion. (Actually, there is a 
very slight decrease in the amplitude of the current density at s 1 - 312 m)„ 
Van Bladel (1964, p. 388) points out that although the total current in the 
neighborhood of a comer remains finite, the' current density becomes 
infinite at the comers. Thus, we may expect an increase m the magnitude 
of the current density wherever the radius of curvature becomes small. 

For this reason, a small sampling interval is required m these regions 
so that the current density is defined properly. 

If we wish to study scattering from a rectangular slab, we should not 
attempt to place a sampling point at a corner. Instead, Andrea sen 
(1964) recommends that the current density should be sampled evenly to 
each side of the comer.- 

The scattered magnetic field intensity calculated from these 
sampled values of the current density has been plotted m Figs. 27 and 28. 
It is important to note that the phase of H is dependent upon the position 
of the observer in space. This result will be true for any cylinder 
when the electrical parameters of the earth are similar to those of this 
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example. 

In Figs. 29 through 32, analytical results are given for the 
horizontal magnetic field intensity scattered by a perfectly conducting 
circular cylinder in a whole r space having various conductivities. We 
have assumed normal incidence of an Ey-polarized plane wave in this example 
and that: the depth to the top of the cylinder is 20 m, the cylinder 
radius is 100 m, and the incident field frequency is 1000 hz. 

It is seen that for extremely low conductivities of the whole-space, 
the magnetic field intensity has a large amplitude and is predominantly 
out-of-phase with the incident field. As the conductivity of the whole- 
space increases, we observe that the amplitude of the magnetic field inten- 
sity decreases, the in-phase component becomes larger than the out-of-phase 
component and then becomes smaller again, and the peaks of the vertical 
field shift closer to the crossover. 

The decrease in amplitude and change in phase of the scattered field 
are caused by both the decrease m wavelength with respect to the radius 
of the cylinder and by an increase in the attenuation of the fields as 
the conductivity of 'the whole-space increases. This can be seen intuitively 
from equations (3-49a) . When the wavelength in the whole-space is very 
much greater than the radius of the cylinder, jk^Rj « j and the order 
zero of the Hankel function becomes dominant. Thus, using the small 
argument expansion for rf^(kjR), the surface current density is given by 
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(3-65) 


and we see that it is independent of and 90° out-of-phase with the 
incident field if k^ is real (which is true when the whole-space has a 
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low conductivity) and h is small with respect to the wavelength. As 
a result, the magnetic field intensity near the cylinder will be large-, 

90° out of phase with the incident field, and the field will be that of a 
line of current. 

In Figs. 29 and 30, the wavelength of the incident field is always 
very much greater than the radius of the cylinder and, as a result, the 
magnitude of the scattered field is large for the three cases. However, 
a decrease in the amplitude of the out-of-phase component and an increase 
m the rn-phase component is observed since higher order terms become 
important in equation (3-49a) as the conductivity of the whole-space 
increases and since becomes larger m equation (3-65) . In addition, 
there is some phase shift as both the incident and scattered fields 
propagate through the earth. 

In Figs. 31 and 32, we see. that the amplitude of the scattered fields 
becomes significantly smaller as the wavelength of the incident field 
approaches the radius of the cylinder and as the skin depth (S') decreases. 
We see also that the peak of the vertical field component in Fig. 32 
shifts towards the crossover as the conductivity of the whole space 
increases. The reason for this is that the surface current density shifts 
to being predominantly on the illuminated side of the cylinder as the 
wavelength decreases. When the wavelength is very much greater than the 
cylinder radius, the current density is constant around the cylinder and 
this current radiates like a line source located at the center of the 
cylinder. However, as the wavelength decreases, there is an increased 
tendency for the current density to be zero on the shadow side of the 
cylinder, as is seen in Fig. 11, so that the location of an effective 
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line source is above the center of the cylinder (i.e. nearer the observer). 
Thus, the peak in the vertical field component will be closer to the cross- 
over. In addition, we observe rn Figs. 31 and 32 that the spaci'al wave- 
length of the fields becomes apparent in the data for high conductivities 
of the whole-space. 

Since the wavelength (A) and the skin depth ( £') in the whole-space 
determine the behavior of scattering from perfectly conducting cylinders, 
a similar trend will occur as the frequency of the incident field increases. 
In Figs. 33 and 34, we have assumed normal incidence of an Ey-polarized 
plane wave and that: the depth z^ to the top of the cylinder is 20 m, 

the cylinder radius is 100 m, and the conductivity of the whole-space 

-3 

is 10 mhos/m. It is seen that as the frequency increases, the amplitude 
of the magnetic field intensity decreases and the peaks of the vertical 
field shift closer to the cross over, as was observed in Figs. 29 through 
32. Note, however, that the m-phase component is larger than the out-of- 
phase component at low frequencies since the ratio wave length/ skin depth 
is as large as in Fig. 31 for <57 = 10 ^ mhos/m. Nonetheless, a large 
amplitude is observed at low frequencies since the wavelength is very much 
greater than the radius of the cylinder. 

In Fig. 35, the conductivity of the whole-space has been set equal 
to zero so that the effect of -the wavelength alone in determining the 
magnitude and phase of the scattered magnetic field can be studied. 

It is seen that as the wavelength decreases, the amplitude of the scat- 
tered magnetic field decreases. However, unlike scattering from perfectly 
conducting cylinders in a conductive whole-space, the out-of -phase 
component is always larger than the in- phase component. 
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The results of Figs. 29 through 35 haye been summarized in Figs. 

36 and 37 by plotting the peak value of |H ) against log (A / R) and 

X 

the phase of the peak value of JhJ against log (A/R) for two ratios of 
wavelength to skin depth. Note, however, that these results assume that 
the conductivity of the cylinder is infinite. If finitely conducting 

I 

circular cylinders are to be examined, the wavelength inside the cylinder 

must be considered also. Thus, a third axis of log ( A ^/R C y^.) could 

be plotted out of the page and the area between A/& =0 and A /£ = 6.28 

in Figs. 36 and 37 would become a volume distribution. 

Fig. 36 illustrates that the general effect of wavelength upon the 

magnitude of the scattered field is to increase )h 1 as A/R increases 

x 

for a fixed point of observation. Consequently, as the frequency of the 

incident field is increased, the magnitude of the scattered field will 

decrease. However, Fig. 36 points out also that if Ai 5R, the magnitude 

of the scattered magnetic field intensity is larger in a conductive earth 

than if the same wavelength had been used in free space. 

• This interesting result is due to the fact that the magnetic field 

intensity transmitted at the fictitious earth-air interface increases' in 

magnitude as the conductivity of the earth increases. Inasmuch as the 

transmitted electromagnetic field is considered to be the incident field 

m these examples, this indicates that the equivalent electric surface 

current density will increase with conductivity, whenever attenuation 

t 

is not important since = H g at the surface of a perfectly conducting 
cylirlder . 

This effect can be understood best by studying the transmitted 
magnetic field intensity. If we assume that k is equal to k Q in air, 
then it is evident from equation (3-51) that the magnitude of the trans- 
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mitted horizontal magnetic field intensity will be given by 



2 

l + Z, M. 




(3-66) 


However, m a conductive medium where A/ S = 6.28 and A = A of 
1 air 

equation (3-66), k will be given by k Q + ik^ and the magnitude of the 
transmitted magnetic field intensity will become 
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if attenuation is not significant. Thus, it is „evident from equation 
(3-67) that for a given wavelength, | | is larger in a conductive 

half-space than in air whenever attenuation is not too important. 

Intuitively, it is expected that attenuation will not be important 
for R and 1 kR| 1 . Consequently, it follows that the curve of 
! | versus log (A/R) will be displaced downwards with increasing A / <$ 

for ■ A ^ 5R and displaced upwards with increasing A/ £ for X 5R. 

Fig, 37 illustrates the variation of the phase of H with increasing 


x 

A/r. it is seen that there is a minimum in the plase of H x for 

A~ 30R and that the effect of the conductive half-space is to reduce 

the phase of H . 

x 

Figs. 38 and 39 compare the analytical solution with an approximate 
solution given by Meyer (1963) which sometimes has been used to predict 
plane wave scattering by conducting cylinders in a whole-space of low 
conductivity. For comparison, we have assumed normal incidence of an Ey- 
polarized plane wave and that: the depth to the top of the cylinder 

is 20 m, the cylinder radius is 100 m, and the frequency of the incident 
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field is 10 tnhos/m. It is evident that for perfectly conducting 
cylinders, Meyer's solution fails to predict the amplitude and phase of 
the scattered field even for free space exterior to the cylinder. The 
reason for this discrepancy is that his solution ihas not accounted for 
the contribution of the incident electric field intensity, which is of 
paramount importance in a whole space of low conductivity. 

3-5 Accuracy and Limitations 

The validity of the integral representations 'for scattering from 

perfectly conducting cylinders has been demonstrated by comparing the 

numerical results with the analytical results for the case of scattering 

i 

from circular cylinders. In general, however, the accuracy of the 
numerical results can be demonstrated only .by observing the solution 
convergence as the number of sampled values is increased. Mei and Van 
B’ladel (1963b) and Andreasen (1964) have suggested that about 10 sub- 
divisions per wavelength are sufficient to give a convergent result. As 
a general rule, however, this seems to be applicable only to scatter ers 
m free space. 

In Figs. 40 and 41, numerical results are compared with analytical 
results for the horizontal magnetic field intensity predicted when a 
cylinder has been sampled 14 times and 36 times per wavelength. It is 
seen that only when the cylinder has been sampled 36 times per wavelength 
do we obtain a maximum error of less than 1%. Thus, we stress that only 
by checking the convergence can we ensure accuracy of the numerical 
results. A general rule on the sampling interval could lead to an unrec- 
ognized error since the skin depth and the wavelength interact to influence 
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what sampling interval is required. 

In Fig. 42, numerical results are compared with analytical results 

v 

when the cylinder of Figs. 40 and 41 has been sampled 14 times and the 
order n of Simpson's rule across each interval is 2 and 4. It is seen 
that increasing the order n from 2 to 4 increases the maximum error in 

i 

the predicted results. However, since no significant change is observed 
when n is increased to 6, this remaining error must represent the inac- 
curacy of the current density representation. 

It should be noted that increasing the number of sampled values of 
the current density is not always the most economical way to increase 
solution accuracy. The reason for this is that it is more desirable to 
obtain solution convergence by increasing the order n m Simpson's rule 
rather than to increase the number of sampling points whenever the 
parabolic approximation to the current density is accurate and the 
scatterer contour is described accurately. This procedure is more 
economical since the computer time required to solve a set of N linear 

equations is proportional to N , whereas the computer time required to 

2 

set up the coefficient matrix is proportionsl to the product of N and 
the order of Simpson's rule. 

In Figs 40 and 41, it is seen that a small error is obtained in the 
predicted result when the cylinder is sampled 36 times per wavelength 
and that an n of 2 is sufficient. . Although this sampling interval improved 
the current density representation, it also increased the integration 
accuracy to greater than that of Fig. 42 with n equal to 4. Thus, it is 
evident in this example that once the current density is sampled 
accurately, the order 2 in Simpson's rule is sufficient for integration 
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accuracy. We note thereby that increasing the number of sampled values 
of the current density is the best general test since it checks both the 
integration accuracy and the accuracy of the current density represen- 


tation. 



CHAPTER 4 


SCATTERING FROM CYLINDERS WITH LOW SURFACE IMPEDANCE 

In Chapter 3, the problem of scattering from perfectly conducting 
cylinders was examined assuming both Ey and Hy polarizations . In 
addition, the representation of the unknown surface current densities 
and the solution of the resulting singular Fredholm integral equations 
was discussed. In this chapter, we wish to consider the more general 
case of scattering from cylinders which have both a low surface impedance 
and a small curvature compared to the attenuation per unit distance. 

Andreasen (1965b)' has considered this problem, but it should be noted 
that although he states that the field must not penetrate deeply into 
the scattering body, his solution also requires that the contour cur- 
vature is small compared to the attenuation per unit distance. Mitzner 
(1967, 1968) has extended this work somewhat by introducing curvature 
dependent boundary conditions. However, as in Chapter 2, the analysis 
of 'this section closely follows that of Andreasen, with the principle 
difference being the physical situation to which the solution, is applied. 

4-1 Derivation of the Integral Equations 

It was shown in Section 2-5 that if the surface impedance of the 
scatterer is low with respect to -the surrounding whole space, and if the 
curvature of the scatterer is small compared to the attenuation per unit 
distance, then the axial field components in the conductive whole-space 
can be written as 
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/ 

where the earth parameters have been referred to by the subscript 1. and 
the cylinder parameters have been referred to by the subscript 2. 

As m Chapter 3, the unknown current densities, Ky and Ks, are 
utsuermined by enforcing the boundary conditions on tangential E and 
H and solving the resulting integral equations. For finitely conducting 
bodies, these boundary conditions are 
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and 


x ( H* - ~ o- 


(4-3b) 


Since we assume that the surface impedance of the scatterer is low 
with respect to the surrounding whole-space and that the curvature of 

l 

the scatterer is small compared to the attenuation per unit distance, 
the representation of the fields inside the conductor is particularly 
simple. Thus, we can write from section (2-5) that just inside the 
contour boundary 


< - - ** S , 

H t - - K " 


(4-4a) 


(4-4b) 

If we let the point of observation approach a boundary point given 
by the position vector and enforce (4-3), then equations (4-1) and 
(4-2) reduce to the desired coupled integral equations 
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(4-5) 
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The integral equations can be solved by dividing the contour into 
N intervals and expressing .each integral as a weighted sum of N sampled 
values of K ( jj ) and K g ( jj ) . A system of 2N linear equations will result 
and these can be solved for the unknown sampled values of K (p ) and 

y v 

K g (^f ). However, the solution of coupled integral equations will -be 
discussed in Chapter 5, since in this chapter we will discuss only E - 

y 

polarized fields which are incident normal to the 'conductor . In this 
case, equations (4-5) and (4-6) uncouple. 

If we assume that k^ = 0 in equations (4-5) and (4-6), then the 
integral equations uncouple and can be solved separately for each axial 
component. We find that 


EjV) * Zjiptyf'') +Ai >f vp 
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For simplicity, we will confine our study on this chapter to - 
polarized incident fields (equation (4-7)). It should be noted, however, 
that the numerical solution of equation (4-8) is the same as the numerical 
solution of equation (4-7). Its solution would not present any numerical 
complications . 

It is interesting to note that the first integral of equation 
(4-7) is the field scattered by a perfectly conducting cylinder and that 
the additional terms account for the finite conductivity of the scatterer. 
It was seen in Chapter 3 that on the basis of physical reasoning, the first 
integral could be thought of as the electric field due to axial currents 
at the surface of the conductor. These axial currents are equal to the 
transverse tangential magnetic field intensity and are present even at 
the surface of perfectly conducting scatterers. 

In a similar manner, the additional terms can be thought of as the 
electric field due to transverse magnetic currents at the surface of the 
conductor. These transverse currents are equal to the axial electric 
field intensity and are present only at the surface of finitely conducting 
scatterers. As a result of the approximations made in deriving equation 
(4-7), we have been able to set the transverse magnetic currents equal to 

-y?> vt>- 



Thus, we could have derived equation (4-7) on the basis of physical 
reasoning alone. The portion of the scattered electric field due to an 


axial electric current would yield a term given by equation (3-8). 

Similarly, the portion of the scattered electric field due to a transverse 

i 

magnetic current would yield terms given by equation (3-10) with H (»*) 

. y 

i „ 

replaced by E (p ) and K replaced by M . 

y * s s 

The numerical solution of equation (4-7) is the same as that of 
equation (3-8) or equation (3-10). Thus, the contour is divided into 
N sampled values of K (p ) and each integral is expressed as a weighted 

y p 

sum of N sampled values of K^(|'). By forcing the boundary condition to 
hold at the midpoint of each interval, a system* of N linear equations will 
result and these can be solved for the unknown sampled values of K (7 ). 

y r 

Integration through the point of singularity in equation (4-7) 
also presents no problem since the contribution of each integral in the 
singular interval has been examined in section (3-2-3). Thus, we can 
rewrite equation (4-7) as 

1 

Ejtp'l - ZijCriVP + AiL f VP Hjtkilp-f't)'**' 

1 A -'/l 


f y?'> & if-*) 

4 Jf. 


(4-9) 


where it is assumed that the contribution of the first integral in the 
singular interval is given by equations (3-31), (3-32) and (3-33) and 
that of the second integral is zero. 
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Once equation (4-9) has been solved for N sampled values of Ky(j3), 
the scattered magnetic field intensity is given by equation (2-87) with 
= 0. We find that 
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K_ tj?') I ^ s( ^f ^ i) 
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- - coU2p,-«)z) Z Js'. (4-10) 

As in Chapter 3, it is necessary to make a small argument expansion 

to the Hankel function and integrate analytically whenever the point of 

> 

observation is close (|k H p - p'j - -3) to the interval over which the 

integration is being performed. Both integrals of equation (4-9) have been 

i 

examined under these 'conditions in a similar analysis m section (3-4), 

Thus, the first integral is replaced by equation (3-58) and the second 
integral is replaced by (-Z ) times equation (3-64) . 

The first integral of equation (4-10) was examined also in section 
(3-4) and should be replaced by equation (3-61) whenever the point of 
observation is close to the contour. However, the second integral has not 
been examined and it is shown m Appendix D that in general the contribution 
of its jth interval to the scattered magnetic field intensity can be 
written as 



J 



S ( C CoS oL,+ -_i [ ( & S'p) ivwit + 2.S (yc^-Sot) 

k. 





75 


^ S(Csvv^- <Ccsoij) + S^E-f) C^S^j- 2S6 


si'V'C’O 2 


where 


C- “ f* ^ -• b ) (X7 “ c j - 1 ) - ( ^ -2c- X 7 +cjj h )l 
L ?cr - J 


+ *3*. jzz5.il* * iliS 

t i Aj L A j 1 


^ A j k ls-» \ X i C3 7-' c i 

h • J ^ L Aj 


- Xg - 2c- t T 7 + c. 7 


(4- 12a) 


3 \*S - iL±ltr t-C- 1 *) - ll 7 -zc 
A i 


lIl±&eJ 1 * 3 
% J 


+ k jj<. + ii 

T^i 


niaiiis] - a j 

tj J TTXfl 


(4- 12b) 


^ * j" ^ ' 1>V ^ T to ~ %c \ ?/o + c j -£? *} 

LA- t J 


V 3 ^ J 

l±r~i Aj 4j J 

(4-12c) 



(4-12e) 



76 


and 



(4- 13a) 


(4- 13b) 


(4- 13c) ' 


(4- 13d) 



(4- 14a) 


(4- 14b) 


(4-14c) 


(4-14d) 


(4-14e) 
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The undefined parameters are given by equations (3-59) and are defined in 

fig . 21. 


Using the small argument approximation for 
H ( h J + s' ) given by equations (3-26) and (3-27) ^ 
(4-14) have been evaluated in Appendix D. 


/<5*+ s ) and 
equations (4-13) and 


4-2 Numerical Examples: Cylinders in a Conductive Whole-Space 


Equations (4-9) and (4-10) were programmed assuming that the incident 
field is the transmitted field of a plane wave incident normal to a con- 
ductive half-space (equation (3-51)). Their validity can be demonstrated 
by comparing the numerical results with the analytical results obtained 
for the case of scattering from circular cylinders. 

It is shown m Appendix C that the transverse magnetic field intensity 
which is scattered by a finitely conducting circular cylinder m the presence 
of an E -polarized plane wave is 


and 



(4- 15a) 
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R xs the radius of the cylinder 

and the derivatives in (4- 16a) and (4- 16b) are with respect to the 
argument kR. 

The equivalent surface current densities around the surface of the 
cylinder are 


(4 ' 17a) 


and 




*L V oo 

‘ i rN.J^U.lO+aX U,R» , 


(4-17b) 


where a^ is given by equations (4-16)^ and the derivatives in equation 
(4- 17a) are with respect to the argument k^R. 

In the case of scattering from conductors which have a low surface 
impedance and which have a small curvature compared to the attenuation 
per unit distance, it is expected that 

, (4-18) 

The distance h in equations (4-15) and (4-17) is the depth from the earth- 

air interface to the center of the cylinder. 

* 

Rigs, 43 and 44 demonstrate that the integral representations derived 
for E scattering from highly conducting scatterers are valid. We have 

y 

assumed normal incidence in this example and that: the depth to the 

top of the cylinder is 20 m ; the cylinder radius is 100 m, the incident 
field frequency is 1000 hz, the conductivity of the whole space is 10"3 mhos/m 
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and the conductivity of the cylinder is 10 mhos/m. The equivalent electric 
surface current density has been sampled 40 times and, for this example, 
a maximum error of less than .3% has been obtained. 

In Figs. 45 "and 46, it is shown that the curvature must be small 
compared to the attenuation per unit distance if equations (4-9) and (4-10) 
are to predict the correct results. In this example, we have assumed 
normal incidence of an E -polarized plane wave and that: the depth z 

y 1 

to the top of the cylinder is 20 m, the cylinder radius is 100 m, the 

incident field frequency is 1000 hz, the conductivity of the whole space 

-3 

is 0. mhos/m, and that the conductivity of the cylinder is 10 mhos/m. 
Since the impedance contrast is large, it is valid to assume that the total 
field at any point just inside the scatterer may be represented by a wave 
propagating normally away from the contour. However, since the curvature 
is 0.01/m and the attenuation per unit distance is 0.002/m, the small 
curvature approximation has been violated. As a result, the convergent 
numerical results predict a field whose in-phase components are m error 
by 20% and whose out-of-phase components are in error by 1000%. The 
numerical results are like those of a scatterer which is more highly 
conducting than the actual inhomogeneity. 

4-3 Numerical Examples: Topographic Scattering 

The integral representations derived in Chapter 2 can be used also 
to investigate topographic scattering problems by considering the contour 
to be very large with respect to the wavelength of the incident field 
and the topographic region of interest. As becomes large, the fields 
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will be determined only by the interface on the upper part of the cylinder 
and, as a result, they will become independent of the contour chosen out- 
side of this central region of interest. Thus, the contour integral 
can be replaced by an integral over (-oo,°o) and, assuming k^ « 0, 
equation (4-1) becomes 

00 ky(|fO JU' 

oO 


VfJ- J r 



(4-19) 


where parameters of the air are referred to by the subscript o, and 
parameters of the earth by the subscript 1. 

The equivalent surface current density in equation (4-19) can be 
thought of as being made up of two-partsJ 


1) a portion due to fields scattered by the topography, plus 

2) a portion due to fields reflected by the half-space. 


Since the topography will be confined to a central section of the contour 
in all examples that will be discussed, it follows that outside a region 
bounded by (-a, a), the equivalent electric surface current density arises 
from fields reflected' by the half-space alone. Thus, outside of the region 
(-a, a), the equivalent surface current density is known and is obtained 
from the analytical solution to the incident field impinging upon a flat 
half-space. Consequently, equation (4-19) can be rewritten as 
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where s is the x coordinate of the position vector , 
i 

i 

and K is the equivalent current density which arises from 

y 

the incident field alone. 

In this way, an infinite integral equation has been transformed 
into a finite integral equation plus several infinite integrals. It 
should be noted, however, that no approximations were made in going from 
(4-19) to (4-20) since the interval (-a,a) always is taken large enough 
to obtain a convergent solution in the region of interest. It does 
assume, nonetheless, that the half-space contour is flat outside of 
(-a, a) . 

It is possible also to transform the infinite integrals to integrals 

that can be evaluated numerically over a finite interval. By rewriting 

the infinite integrals as 
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it follows from (4-19) that for p above the half-space, the integrals 
• over (-oQ, otj ) must equal the electric field intensity reflected by a flat 
half-space. Thus, equation (4-20) reduces to 
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It is important to remember, however, that m equation (4-22) the contour 

i 

integrals involving are along a flat half-space while the contour 

integrals involving are along the topographic, profile. In addition, 

(4-22) assumes that |T is above the half -space. 

By following an analysis similar to section (4-1), it is found that 
the desired integral equation is 
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As in earlier examples,. we will confine our study to normal incidence 
of an E -polarized plane wave. Thus, it follows from equation (2-38a) 

y 

and Stratton (1941, p. 493) that 
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It is interesting to note that equations (4-24) and (4-25) can be 
obtained analytically from equations (4-19) and (4-23) if we assume that 
an E -polarized plane wave is incident normal to a flat half-space. In 

y. 

this case' the equivalent electric surface current density is constant 
and' equation (4-19) becomes 
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It is shown in Appendix E that 



U' 


, c 1 

I- a 





(4-27) 



84 


provided Re(k), Im (k) 0, Re(S)>0, and Re(|j.)> -1. 

Setting 5 = 0 and |j. = -1/2 in (4-27), we find that 
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and setting D= 1 and p. = -1/2, 
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Thus, equation (4-26) becomes 
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Taking the limit of (4-30) as § approaches zero and introducing the 
result into (4-23), it is found that with (-a^a) equal to (0^0), ihen 


ejtf) = 4±z, kj . 


(4-31) 


Note, however that to obtain (4-31), ^*3 — l. in 

equation (4-23) has been replaced by Z-^y(^‘ i ) ,K (^0 • This has been done 
since the singular contribution of the last term of (4-23) had been added 
to (4-23) in deriving the integral equation and now must be subtracted 
before the last term can be ignored. 

Thus, it follows from (4-31) that 
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(4-32) 
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which agrees with (4-24). Introducing (4-3i2) into (4-30), we find that 
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(4-33) 


which agrees with (4-25) . 

In deriving (4-22) and the subsequent results, it was assumed that 

the point of observation was above or level with the flat half-space. 

However, if the topographic profile includes a valley which extends below 

the half-space, it is not valid to replace the integrals over (-co,oo) 

rsf I gc _ 

in equation (4-21) by E whenever A is below the half-space. 

y ' 

lo establish what value the integrals do yield, it is necessary to 
remember that m applying the equivalence principle in section (2-2), 
the same fields were obtained only within the volume of interest. For a 
point of observation outside the volume of interest, the predicted fields 
will be zero since the boundary conditions have been satisfied by the 
equivalent surface current densities (see Harrington, 1961, p. 106). 

Thus,- the total field will be zero when the point of observation is below 
the half-space, and it follows from (4-19) that the integrals over (-. 10 , 00 ) 
must be equal to the negative of the incident electric field intensity. 

This result can be obtained analytically if it is assumed that the 
incident field is an E^-polarized plane wave incident normal to a flat 
half-space. In this case, the equivalent electric surface current density 
is constant and for a point of observation below the half-space, the infinite 
integrals of (4-21) become 
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Using the results of (4-28) and (4-29) , equation (4-34) yields 
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Introducing the analytical values for from equation (4-24), the 
contribution of the infinite integrals is found to be 


. ‘ Z " H " ■*" > (4-36) 

which is the negative of the incident field. 

We can summarize these results by rewriting the integral equation 

(4-23) as 
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where -E (^" *) is the electric field intensity reflected by a 

flat conductive half-space and is to be used when the point 

of observation is above or on the flat half-space, 
i 

E (p" - ) is the incident electric field intensity and is to be 

y r 

used when the point of observation is below the flat half- 


space. 



87 


It should be remembered that the contour integrals in (4-38) are along a 
flat half-space while the contour integrals in (4-37) are along the topo 
graphic profile. 

If we assume tint an E -polarized plane wave is incident normal to 

y 

a flat half-space, then K 3 ' is constant and equation (4-38) can be re- 

y 

written as 
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(4-39) 


Special care must be taken at the edges of the interval in which 
K is assumed *to be unknown to ensure that the unknown parabolic current 

y 

distribution is continuous with the known current distribution. This 
is accomplished by assuming that the Oth and Nth + 1 intervals lie out- 
side (-a, a), that and + ^should be approximately equal to 
and to ensure that the parabolic fit to the current density is not 
violated. 


If, however, the incident field is an E^-polarized plane wave 
incident normal to a half-space, then it is valid to assume that W 


o 

and W, , are infinite since the parabolic fit to the current density 
N + 1 

is not violated. In this case, consider the general integral 


M > (4 ' 40) 

W- 
J J 

where A is a constant and F ( ^ ) is the kernal of the integral, 

to represent each integral along the topographic profile in equation 
(4-37). Then the contribution of the first interval to (4-40) is 
(compare with equation (3-22) as ^ — ►*«>) 
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Similarly, the contribution of the last interval ^ to (4-40) is 
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In particular, when the boundary condition is applied in the first 
interval, (4-41) reduces to 

Ak 3.rV»r ,A, '> (4-43) 

s,-w, 

and when the boundary condition is applied in the last interval, (4-42) 


reduces to 






(4-44) 
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An expression for the scattered magnetic field intensity is obtained 
by following a similar analysis for equation (4-10). Assuming that the 
contour is flat outside the interval (-a, a) and that the equivalent surface 
current density arises from the incident field, then it would be found 
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where s^ is the x coordinate of the point of observation, 

H tr (p ' ) is the magnetic field intensity reflected 
by a flat conductive half-space and is to be used when the 
point of observation is above or on the flat half-space, 

-H (p'Xs the incident magnetic field intensity and is to be used 
when the point of observation is below the flat half-space. 


and 
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(4-46) 

1 

As in equations (4-37) and (4-38), the contour integrals involving 
are along a flat half-space while the contour integrals involving K 

y 

are-along the topographic profile. 

It is shown in Appendix D that the trigonometric expressions can be 


written as 
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where the parameters ’ are defined in Fig. 21. 

Thus, expanding the trigonometric expressions according to equations (4-47), 

I 

noting that the half-space is parallel to the x-axis (sc = -90°), and 
assuming that the incident field is an Ey-polanzed plane wave incident 
normal to the half-space, then equation (4-47) can be rewritten as 
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where I 1 , I , I J , I , I _ and I are given by equations (3-63a), 
•i- 2 5^ 6 9 i-U 11 

(3-63b), (4- 13a), (4- 13b), (4- 14a), (4-14b) and (4- 14c) 
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* respectively, with a = -a-s < and b" = a-s t • 
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It should be noted that equations (4-40) through (4-44) must be 
considered at the edges of the interval (-a, a) to ensure continuity of 
the known and unknown current distributions in equation (4-45) . 

Equations (4-37) and (4-45) were programmed assuming that the 
incident field is an Ey-polarized plane wave incident normal to a conduc- 
tive half- space. Their validity can be demonstrated by comparing the 
numerical results with the analytical results obtained for the case' of 
reflection from a conductive half-space, equation (4-24) . 

It should be noted, however, that it is not necessary to sample the 
contour at least ten times per wavelength to describe adequately the 
equivalent surface current density while testing the validity of the 
integral representations. The" reason for this is that with normal inci- 
dence the equivalent surface current density is constant. As a result, 
the contour can be described by any number of intervals if integration 
accuracy is maintained. 

To take advantage of this fact, the integrals over the contour 

section (-a, a) m equation (4-37), were replaced by equations (3-58) and 

(-Z ) times equation (3-64), and those in equation (4-45) were replaced 

2 a 

by equation (4-11) , In this form, all independent integrals in equations 
(3-58), (3-64) and (4-11) were integrated numerically with the integration 
tolerance being specified by an integration sampling rate per wavelength. 

The concept of an integration sampling rate per wavelength was 
chosen to correspond to the integration accuracy that would be obtained if 
the contour had been sampled evenly at this same rate. For example, if 
a contour section one wavelength long had been sampled previously by 10 
intervals of equal width, then the same accuracy can be achieved by one 
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interval with an integration sampling rate of 10. This is equal to an 
n of 20 in Simpson's rule for this interval width, but in general n will 
be different for each interval, although the integration sampling rate 
will be the same. 

Table 1 demonstrates that the integral representations are valid 

for topographic scattering problems and summarizes the accuracy and 

convergence of the solution as the numerical variables are changed. It 

should be remembered that since the incident field is normal to the 

interface at all points and that since the curvature is zero, equations 

(4-37) and (4-45) represent an exact formulation of this particular problem. 

Thus, Table 1 is a good indication of solution accuracy and convergence 

under various numerical conditions. 

In this example, we have assumed that an E^-polarized plane wave is 

normally incident to a horizontal interface and that: the height z^ of 

the observation point is 100 m above the interface, the incident field 

frequency is 1000 hz, the conductivity of the upper region is 0, and the 

-3 ' 

conductivity of the lower region is 10 mhos/m. In particular. Table 1 

assumes that the point of observation is at x - -400 m. This point was 

chosen over x — 0 (the center of the ha If -space contour) since H always 

-14 

is given by a computer round-off error of about 10 at the center of any 

contour symmetric with respect to the z-axis. Thus, by choosing x = -400 m, 

H is proportional to the accuracy of the solution rather than the computer 
z 

round-off. 

In Table 1, "a" is the half-width of the contour section (-a, a) 
within which the surface current density is assumed to be unknown, N is 
the number of sampling points of the equivalent surface current density 
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± a 


Real (H ) 


Imag. (H x ) 


Real (Hz) 


Imag (Hz) 


N 


Analytical Results 


,9894 


.01251 


Int . rate 


Re 


Numerical Results 


600,000 


1.0127 

.9944 

.9897 


.00511 

.01104 

.01266 


8x10 


-3 


2x10 


-3 


7x10 


-5 


2x10* 


5x10 


-4 


2x10 


-5 


14 

14 

14 


10 

20 

20 


.3 

.3 

.5 


300,000 


1.0124 

.9944 


.00528 

.01105 


9x10 


-3 


-.-3 


1x10 


-3 


.9942 

.9917 

.9897 

.9896 

.9895 


.01103 

.01280 

.01268 

.01267 

.01267 


1x10 


-3 


3x10 
5xl0~ 4 
4 


12 

12 


5x10 


5x10* 


1x10 


-5 


4x10 


-7 


3x10 


-6 


7x10 


-5 


1x10 


-5 


8x10 


2x10 


-5 


12 

22 

22 

12 

12 


10 

20 


30 

10 

20 

20 


30 


.3 

.3 


.3 

.3 

.3 


.5 

.5 


50,000 , 


1.0129 

.9943 

.9908 

.9896 


.00520 

.01093 

.01275 

.01264 


8x10 


-3 


3x10* 


2x10 


5x10* 


1x10 


-4 


3x10 


8x10 


1x10 


8 

8 

8 

8 


10 

20 

10 

20 


.3 

.3 

.5 

.5 


20,000 


1.0131 

.9946 

.9897 

.9896 


.00546 

.01103 

.01269 

.01268 


9x10 


3x10* 


2x10 


5x10* 


9x10* 


8x10 


-6 


5x10 


-6 


3x10 


-6 


6 

6 

6 

6 


10 

20 

10 

20 


.3 

.3 

.5 

.5 


5,000 


.9894 

.9894 

.9894 


.01251- 

.01251 

.01251 


1x10 


-14 


1x10 


-14 


1x10 


• 14 


1x10 


-14 


1x10 


-14 


1x10 


-14 


4 

4 

4 


10 

20 

20 


.3 

.3 

.5 


_ j ^ 

Table 1. Magnetic field intensity reflected by a horizontal interface between 


air and a conductive earth: Ey-polarization 3) x = -400 m, Z Q = 100 m, 

f = 1000 hz, =0.mhos/m, and — 10 mhos/m. 
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within (-a, a), int. rate is the integration sampling rate per wavelength 
and Re is that value for which the small argument solution to equations 
(4-37) and (4-45) is used if (k^rj ^ Re over the entire width of a sampling 
interval. 

It is evident from the data of Table 1 for a = 5000 that equations 
(4-37) and (4-45) are an accurate formulation of the problem whenever 
jk Q a|<<: Re/2. This is not surprising since the small argument approxi- 
mations are very accurate for these small arguments and they are used to 

t 

evaluate all integrals when the contour section (-a, a) is so short. 
However, when jk aj > Re/2, both Simpson's rule and the small argument 
approximations are used to evaluate equations (4-37) and (4-45), in which 
case the usual solution accuracy of 1% or better is achieved. 

It is evident from the data of Table 1 for a = 300,000 that the 
solution accuracy can be increased by raising the number of sampling 
points within the interval (-a, a) or by increasing Re, as well as by 
increasing the integration sampling rate. This result is actually more 
a function of the manner in which the program has been written than 
the numerical nature of the solution. The reason for this is that the 
small argument approximation is used only when jkrj £ Re for all points on 
an interval. 

As an example of why this is a programming effect, it should be 
noted that the small argument approximation would not be used to calculate 
the contribution of an interval to the fields iif|kr a ( ^ .01 and|kr^|yo.4 

(where r and r are the radial distances -to each edge of the interval) 

a b 

and Re were .3. Instead, Simpson's rule would be used and significant 
integration errors would be incurred unless a large' integration sampling 
rate had been specified. ' However, by choosing an Re of .5, the more 
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accurate small argument approximation would be used for this interval 

and solution accuracy would be improved substantially. Alternatively, 

the contour could have been sampled more frequently and the interval of 

this example might now be sampled twice such that jkr^j a/ .01, |kr | A j .09 

and |kr ( . 4. Thus, even with an Re of .3, a more accurate solution 

b 

would be obtained since the small argument approximation would be used 

for the interval (r , r ) whereas numerical" integration would be used for 

a. c 

the interval (r c , r ) . 

It is apparent that the above situation can arise in equation 
(4-37) whenever small interval widths are adjacent to large interval 
widths and/or in equation (4-45) whenever the point of observation is 
close to the contour (|kr|«A). As a result, it is necessary to choose 
small contour sampling widths near the points of observation and gradually 
increase the contour sampling widths on each side away from the points 


of observation. In this way it is possible to sample accurately the 
equivalent surface current density with the smallest possible matrix 

-5 

size 'and still be assured of obtaining a solution accuracy of about 10 
To illustrate the application of the method to topographic 
scattering problems which cannot be handled analytically, the field 
scattered by the hill of Rig. 47 has been considered in Figs. 48 through 
51. In this problem, the contour has been computed in meters according 
to the Gaussian distribution 


He, 3 kt 

In Figs . 



5 IOC • (4-49) 

48 through 51, it is asaimed that an Ey-polarized plane 


wave is incident normal to the topography and that: the point of 
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observation Z is 150 m above the half-space, the incident field frequency 
o 

1 

- is 1000 hz, the conductivity of the upper medium is 0, and the conductivity 

_3 

of the lower medium is 10 mhos/m. 

The vertical magnetic field intensity has been examined in 

Figs 1 , 48 through 50 to establish solution convergence. This field was 

chosen over the horizontal magnetic field intensity since it is more 

sensitive to both the inhomogeneity and round-off errors. 

It is evident from Fig. 48 that for this example, a contour 

section bounded by (-10, 000^4- 10,000) £ie (-X air /30, + A a j_ r /30^, is 

sufficient to describe the unknown surface current density. In Figs. 

49 and 50, the sampling density required within the horizontal region 

(-10,000, + 10,000) has been tested for solution convergence. On the 

basis of Figs. 48 through 50, it was decided that an accuracy to better 

than 1 % (at the peak value of H ) could be obtained if the contour was 

z 

sampled 44 times within the contour section bounded by (-10,000, + 10,000). 

In Fig. 51, the horizontal magnetic field intensity scattered 
by the hill has been plotted assuming a constant flight level of 150 m 
above the half-space and also^ assuming a contour flight level which is 
150 m above the interface. It is seen that by contour flying, the peak 
electromagnetic response of the hill has been reduced by a factor of 
two. In addition, the electromagnetic response has been held to a 
constant over the central portion of the hill. However, even by contour 
flying, there has been an approximate increase in the peak value of 
real (H x ) of 127 as a result of topography. 

It should be noted from Fig. 47, however, that since the curvature 
approximations have been violated in this example, real (H) will be about 
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20% too large whereas imaginary (H) may be meaningless. Nonetheless, 
it has been feasible to estimate these errors since the curvature 
violation and electrical parameters are approximately the same in each 
problem. As a result, it is possible to estimate that for the contour 
flight example of Fig. 51, there should be an actual increase in the 
peak value of real (H ) of about 10% as a result of topography. 

* X 



CHAPTER 5 


SCATTERING FROM CYLINDERS WITH ARBITRARY IMPEDANCE 


In Chapters 3 and 4, an integral equation solution to scattering 
from perfectly conducting and highly conducting cylinders was discussed. 
In this analysis, the general integral representations were investigated 
under simplifying assumptions so that each numerical difficulty could be 
studied individually. 

Having examined some particular problems encountered in formulating 
an integral equation solution to geophysical scattering problems, we are 
prepared now to examine electromagnetic scattering from cylinders with 
arbitrary impedance. Once this analysis has been completed, any two- 
dimensional problem can be formulated and solved (within the time and 
storage limitations of the computer), including the very important 
problem of scattering from cylinders in a conductive half-space. 


5-1 Derivation of the Integral Equations 

*• 


t t 

The general coupled integral representations of p ) and H^(^) 
m the exterior medium were given in Chapter 2 (equations (2-56) and 
(2-57) as 


Ej(j>) =• e*tp) - z, y, 1 f KjCpouf Cv,t^-p'0 d. s ' 
T V ^ 

4k X 

+ l A ( As ^ 

^ c l 98 
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H*(fU Hjip)+ ^ f 

I I c, 

- it,/, f H 5 (f') s~ (f-) H? LVf-t'l) y 

Jc, 

+ iX f Mf'J Mlji-a) , 

^ -Jr. 


(5-2) 


where the earth parameters have been referred to by the sub- 
script 1. 

The transverse field intensities are obtained by adding the incident 
transverse magnetic field intensity and transverse electric field intensity 
to the scattered transverse field quantities (equations (2-66) and 
(2-68)), and we find that 


Hfctp) • t£i P > { yfOF,U l>Pjf ')<h' 

4% Jc, 

- -M, 1 f y " Kif>V<;,j»,f'><h', 

■i ** 1, 

and 

thf) »-E r £ tpUi)^.f m 

-C J C| 

- klkf N,r>h(V, lf , r )i ! ' + a 1 ^f Mf'JFsU,,*,*')!.', 

i, ivi, 


(5-3) 


where 


lp y ij') = Hf Of lp-f '|)( ^ - cosp£ ) ^ 


(5-4) 


(5-5a) 


(5-5b) 
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z ^ ^1 f i j> ) - e-t’sCp* ~ •*) H 0 - cos p,z. ) 


(5-6) 


_ - cosC2f>-*#.)z) . 

V\f-f' ! (5-5c; 

Similarly, the integral representations of the fields inside the 

scatterer are obtained from equations (2-72), (2-73), (2-74) and (2-75). 

It is found that 


- Ki ip 9 ^ W ( ^114 \f-f' i)i s 

«h 4, 

- ±L [ MlWp-.htfpOtlj.f'l) Is', 

■ ^ c l 

Hhjr)* _ ^ f M (s') 

J Hh j c , 

+ i f Elitp') S'"* 

Jc, 

- i_L f K s(f 0 (7-4 Xf U 2 ljr-jI'l)J>' ) 

4 JC , 

tf ) = -iX f 

<£ Jc, 

+ X4 f Mpo ^CX t ?ip')I>' + f fjhijf ,7')^ , 

4r, - 
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Ejlf) « - iX f My(fQlf,LY, l f,pA/ i 

4^2 J C 

+ kjX f . 5 X f 

4 k£ -'c. 4 kj. 4 


(5-9) 


where the cylinder parameters have been referred to by the subscript 2, 
and ) and ) are given by equations 

(5-5) . 

The unknown current components K^(^), K g (|-), M (jp) and M ( ^ ) 
are obtained by enforcing the boundary conditions on tangential E and”H 
and solving the resulting integral equations. We have from equation (2-53) 
that the tangential unit vector is given by 


^ = - SCvr. U l 


/k *■ e°£ <* 


* * > 


(5-10) 


where is the normal into the cylinder at the boundary 
point. Thus, the tangential components of E and H in the transverse plane 
are found by taking the dot product of (5-10) with equations (5-3 ) , (5-4), 
(5-8) and (5-9). We find that exterior to the cylinder, 

( j>) = ^ (p) 4. i X ? Kj (p') 

«F J, 1 


- It,, y f K 5 lf) - jM H s (f')C,lY,jr l p')J«' } 

4 k, Jq ^,W \ 

d 

£^(p) = £/*') &^ x , f )f y > ^ 

- ^ y j 1 4 . ^ j * 


(5-11) 


(5-12) 
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Interior to the cylinder, 



(5-14) 

where S*FCK t p,^') , 

G.ty.j r) - - i^Cyip-j'O cos-ip-^O 

G^^p ,^0 = - jVhC^-^cosC^-^l) W 0 ( Yl p-p'l) 

+ s5*tap— ) 

Xlf-f' 1 

G^C^Pijs-') - - tos (yi-^VosCp-^O H^tyjp-p'i ) 


(5 -15a) 
(5-15b) 

(5- 15c) 


(5-15d) 


If we let the point of observation approach a boundary point given 

// * **^ ^ 
by the position vector p and equate tengential components of E and H, 


we obtain the desired coupled integral equations: 
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A J c 

I 


1. •'"[ 


'«, hi tfftv, i?. P 'i)+ p'Ol^s 

“IT ”” J. J 


- L 

4 


M s (pO css^-sOj"^ + Hi (Xjp - p'OjdU ^ 


Hj tpi- -_L [ M^p') 
AX 


€&< \f - p ' 1 ) + lL $^f-f •) 




I f M,(s') Sinlji-^r k*l H^C^lp-A'O + tfWf-fOW 
4-1, ITET J 

Kc. ip') iff 0 + X. H/lXlp*- p'^)j <^ 3 j 


--M 

4 Jc 


s»H. (p) *» -j- 


I f K,lp[ &,Of„^p + 6, Wilis' 

•+ i ( Ksif) Ml MM* f) + 

+ jf M,<jrtfj£ (SjW^'U j£_6-zCK, (',?'*] ^ > 

A J/' ti r.K #i j.\ ^ 


and 


^ “c 

$ , E 


j (p = -_L f m 5 MJli? s, CX»F'jf Mw M]^ s 
4 4 ij L X 1 

+ X f MfOr 6.W. rtp 4 V*L 6, 

4 -I _ i* ! X 


Mpj^ jM 


f KslfOf 2.X 

\ M 


(5-16) 


(5-17) 


(5-18) 


(5-19) 
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Equations (5-16), (5-17), (5-18) and (5-19) are four coupled integral 
equations with four unknown current components Ky(^), K g (^), My(p ) 
and M g (p ). They are solved by dividing the contour into N sampled 
values of each unknown current components and solving the resulting system 
of 4N linear equations for the unknown sampled values of Ky(jj), K s ("p), 
My((p, and M s ( p ) . The solution of four coupled integral equations is 
similar to the solution of two coupled integral equations, which is 
discussed in section (5-3). 

If we assume that ky = 0, then equations (5-16) and (5-18) uncouple 
from equations (5-17) and (5-19) and can be solved separately for each 
axial component of the total field. However, we will confine our study 
to E -polarized incident fields, in which case it is sufficient to consider 

y 

equations (5-16) and (5-18) alone. It should be noted, however, that the 

numerical solution of equations (5-17) and (5-19) is the same as the 

numerical solution of equations (5-16) and (5-18). In fact, the calcu- 
lation of the scattering matrix would not be much more difficult even 

if ky were not assumed to be zero. The only differences would be that 
a larger storage would be required for the matrix and its inversion 
would be more time consuming. 

Assuming that ky ** 0 and confining our study to Ey-polarized 
incident fields, we find that the integral equations reduce to 



(5-20) 
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and 


A 

■S » 


ky f JfO S ^ )£ Jt, W? Ik, I p "C ^ n) + ^ ^ C ^ i p '-l p ' l)J ^ " 


if M s (.pO H^ } Ck,Jp- p-,)-*. _k, ^UfplpM)U J / 

U, 2_ J 


+ if NlsCpOo.C^-^j.) HJUljfl^)) + H^(.klp-p'i)7JLs 

4 c < i ^\ff zjf-y i - 


(5-21) 


Assuming that the incident field is the transmitted field of a Ey- 
polarized plane wave incident normal to a conductive half-space, then 


£$ If) * ^ 2 . H< 


.ik.s" 


(5-22) 


l + Vz a 


and a 'S"‘ / , a . . -ik,**" (5-23) 

- gjggj Ho jl . 

I + z '/z 0 

Once K^( ^ ) and M s ("jJ ) have been estimated from equations (5-20) 
and (5-21), the total or scattered magnetic field intensity is calculated 
from equation (5-3) or (5-8) by setting = 0. Depending upon which 
region the fields are desired, it is found that 

= tig f Ky (p\) _ Y l)( sUv m - cospz)<ks' 

A \ 


- K 


42, J c, 


f Ms Cf>') Cai(fi-oL,) Ht > Lk,l^-p ,, !)(3VHp / 5c - COS p Jl$ l 

J/l 


42 | Jc ( 


C»&(Ap_-0&)db f ? 


(5-24) 
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exterior to the cylinder, and 

14^M - - ik ' Ky ip") M 

« Jc, 


-f k f M s (pOcos(^OH* J (|Jp-f'0( 

42^ Jc ( 


+ W f Mg(pO H; (jjgj P~ f 7 ^ ^ ***) g )^ S J 

■ c ! 3T lf-f ' * (5-25) 

interior to the cylinder. 

It is interesting to note that the general integral equation (5-20) 
can be reduced to the integral equation for cylinders with low surface 
impedance (equation (4-7)) by a different derivation than discussed in 
sections (2-5) and (4-1) , If it is assumed that becomes very large 
with respect to and that the curvature of the cylinder is small, then 
only that part of the contour in the immediate vicinity of the point 
at which the boundary condition is being enforced contributes to the con- 
tour integral of the interior field representations. Under these con- 
ditions, the contour integral can be replaced by an integral over the 
interval (-oo, oo) since the fields will be independent of the contour 
chosen outside of this region of interest. Consequently, equation (5-20) 
can be rewritten as 



(5-26) 
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Noting that 


Cos = S 

1 


(5-27) 


for a point of observation just inside the cylinder, then the last two 
integrals of equation (5-26) can be evaluated analytically from equations 
(4-28) and (4-29). Thus, equation (5-26) reduces to 

E‘( ? n= ^ i F ") h i r > - «il£2 


- i_ 

A 


1 Jc, 


(5-28) 


On replacing M (p) by “Z 2y Ky as discussed in section (2-5), we find that 


+ ik f kfi'CMf* f '0 ^ s< ‘ 

4 Jr 


(5-29) 


Equation (5-29) agrees with equation (4-7) . 

As in Chapters 3 and 4, it is necessary for integration accuracy 
to make a small argument approximation to the Hankel functions whenever 
IMlpfc'l 1S less than 0,3 and inte S rate th - e resulting expressions 
analytically. Fortunately, the necessary integrals have been studied 
already in Chapters 3 and 4. From equations (3-8) and (3-j.O), we can 
write that the small argument contribution of the jth interval to 
equation (5-20) is 
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+ ikl + kkLV tj jfV')> 

4 ^ (5-30) 

where E (k ; p') is given by equations (3-58), (3-59) and 

yj V r 

(3-60) with current component Ky(p), 
and B _(k,p # , ^') is given by equations (3-62b) and (3-63) 
with current component Mg ( p ) . 

Similarly, the contribution of the jth interval to the scattered 
magnetic field intensity (equation (5-24)) can be written from equation 
(4-10) as 


(?) ' ~ ik 

J 1 4 


( A os oc^ +- 8 & ) ^ + (_rt slvvotj - "B S Cc,s ° £ j 


i T S(c h ^ t> s ^ v ' oi j) ~ — (( ^> 1 ^~ p) + 2^6 Cos.otj)J-4. 

1 A 


+ £>(cs\fi\ oL- - i 2)C&Sc<' • ) 4. J_ I IS^t- F)Cc>$°{\ - iSG- z f ) 

L V J > (5-31) 

where A and B are given by equations (3-62) and (3-63) with current 

component K (p), 

y r 

C, D, E, F, and G are given by equations (4-12), (4-13) 

(4-14) with current component M (jj ). 


and 
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The small argument approximation to equation (5-21) is obtained by 
taking the dot product of s and equation (5-31) with opposite sign. 
Thus, from equations (5-10) and (5-31), the small argument contribution 
of the jth interval to equation (5-21) is found to be 


\ 

s 




where 


oC- ) 4* B S C.&S (cl - 
J t 



(5-32) 



Sun 


— a-j) + 3) Cd i — <=<j ) J 


-J. ^ (S X £- F) Co s 



(5-33) 


where A and B are given by equations (3-62) and (3-63) with 
current component K^(|T)j 

and C, D, ^E, F, and G are given by equations (4-12), (4-13) > 

and (4-14) with current component M (p). 

s r 


5-2 Integration Through the Point of Singularity 


Integration through the point of singularity in equation (5-20) 
presents no problem if it is noted that the limit taken when applying 
the boundary conditions is from both sides of the interface. Thps, we 


find that 
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(5-34) 


for a point of observation just outside the cylinder, and 



for a point of observation just inside the cylinder. Since 


(5-35) 


,Yf; 

tu / Jl£ i - [Ur> 

V 4 J Wi ■ X 


(5-36) 


the contribution of each term in the second integral of equation (5-20) 
adds to zero. The contribution of the first integral of equation (5-20) 
m the singular interval is found by applying equations (3-31), (3-32) 
and (3-33) to each of the Hankel functions and summing the results. 

It should be remembered, however^that equations (3-31), (3-32) and 
(3-33) are accurate only if the small argument approximation to the Hankel 
function is valid. When examining scattering from highly conducting 
cylinders, it is quite possible that although JkJ £ .3, jk^J .3 

where w is the half width of the singular interval. This problem can 
be overcome by going back to equation (3-28) and rewriting HK00 and HK02 
as 



Ill 


- /' w 

\C (UO As' > 


5^0 J J 


(5-37a) 


(5-37b) 


where 

c - .3 (5~37c) 

Tu ' 

Thus, the integrals are approximated analytically over the interval 
(-£,£.) from equations (3-30a) and (3-30c) and numerically over the 
intervals (-W;,£.) and (g^w; ). These results are then inserted 
into equation (3-28) to obtain the contribution to the singular interval. 

In the singular interval, od.= <*, so that equation (5-21) can be 
rewritten as 



(5-38) 
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It is evident that the contribution of the first integral of 

equation (5-38) is zero by comparing it with the second integral of 

equation (5-20). The contribution of the second integral of equation 

(5-38) also is zero and this can be shown by examining the term involving 

,(«) 


H v " (k^ ) j> A ^1) . On expanding M ( p') the contribution of the terms 




\ 
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(5-39) 


(5 -40a) 


(5-40b) 


and 


' S~6 1 i w . ' 


(5-40c) 


Equations (5-40) have been evaluated in Appendix D assuming that 

Ak /££?) is given by equation (3-26) . It is found that on taking 
o i 


the limit 


L . r £ A JU l£j£l A ] 

J S^oL 2 J. Wl . J 


(5-41 a) 
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* O > 


(5-41b) 


and 


I 7 - 0- 


(5-41c) 


Rewriting (5-41c) as 
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(5-42) 


then it is apparent that on taking the limit as 
jCj- - o . 


(5-43) 


Since a similar analysis involving H (k t-'X a 1 1 ) would yield an analogous 

o 2 r 

result/ the contribution of the second integral of equation (5-38) is 
zero . 


footing that 


3. ( & » 2. c<n ( /2> — • I ' 2 £ — I 

TsV 7 ) 


(5-44) 


‘•then the contribution of each term in the third integral of equation (5-38) 


can be written as 
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where 1^, I,^ 1^, 1^, I^q, and are given by equations (3-36) 
and (4- 14) . 
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These integrals have been evaluated in Appendix D assuming that 
H^(k/*% s'*’ ) and lP\k /FTT 7 ^* ) are given by equations (3-26) and 
(3-27). Thus, if the integration is performed over two intervals as in 
equations (5-37), then it is found that 





and 

fa (af r -tJ . jL/a^H^CK) - « Ko °) (5- 46c) 

4 

where HKOO is given by equation (5-37a) 
and £ is given by equation (5-37c) . 

The contribution of the third integral in the singular interval is 
thus found by applying equations (5-45) and (5-46) to each of the Hankel 
functions in the third integral of equation (5-38) and summing the results. 

5-3 Numerical Solution of Coupled Integral Equations 

The coupled integral equations (5-20) and (5-21) can be solved in 
essentially the same manner as when solving the Fredholm integral equations 
of Chapters 3 and 4. The integrals in each equation are approximated by a 
weighted sum of N sampled values of K^(^') and M g (p'). When this is done, 
the integral equations reduce to a system of 2N linear equations of the form 
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(ejl . 


(5 -47a) 

(#0 = 

(Ca.JCKj) + (C B )( Hd , 

(5~47b) 


where (C^), (C^), (C 21 ) and ^ 22 } re P resent square matrices of order 
N and (Ep, (H^,), (K^) and (M g ) represent column matrices of 
length N. 

Being cognizant of the matrix nature of equations (5-47), (K ) and 


y 

(M ) can be obtained by treating (5-47) as two linear equations with two 
s 

unknowns. On solving these equations, it is found that the best form of 
the solution is 


(Mh- J[(c y )(c„rj(c, l ) - (cjj'j[(c ll Xc„)‘](^)- «)} , (5 . 43a) 

and 

Up* (c»r'(Ej) - (c„) ' (C, 2 )(M s ) ; (5 -48b) 

where the superscript -1 denotes the matrix inverse. 

Equations (5*48) are preferred to other forms since the solution is 
well behaved even when studying reflection from a flat half-space. In 
this case, C ^ ~ ^l = ® ancl equations (5-47) uncouple, but equations 
(5-48) still give the correct solution. 


5-4 Numerical Examples 


Equations (5-20), (5-21) and (5-24) were programmed assuming that 


the incident field is the transmitted field of an E^-polarized plane 
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wave incident normal to a conductive half-space (equations (5-22) and 
(5-23)) . As m Chapter 4, their validity is demonstrated by comparing 
the numerical results with the analytical results obtained for the case 
of scattering from finitely , conducting circular cylinders (equations 
(4-15) , (4-16) and (4-17)). 

In Chapters 3 and 4, it was seen that an n of 2 in Simpson's rule 
was sufficient in general for integration accuracy across each interval. 
Although this will still be true in equations (5-20), (5-21) and (5-24) 
for those integrals involving the exterior Green's function, it may not 
be sufficient for those integrals involving the interior Green's function. 
The reason for this is that with large values for the conductivity of 
the scatterer, the modulation of the integrands of the interior integral 
representations by the Hankel functions becomes sufficient to affect 
integration accuracy. 

To overcome this problem, equations (5-20), (5-21), and (5-24) 
were rewritten in the form of equations (5-30), (5-32), and (5-31) 
respectively. In this form the order n in Simpson's rule is specified 
separately for each region and all independent integrals in equations 
(5-30), (5-31) and (5-32) are performed numerically. 

Figs. 52 through 55 demonstrate that the general integral repre- 
sentations are valid, although the accuracy and convergence obtained is 
not as good as that realized m Chapters 3 and 4 for more highly conducting 
scatterers. In this example, we have assumed normal incidence of an 
E -polarized plane wave and that: the depth z to the top of the cylinder 

y L 

is 20 m, the cylinder radius is 100 m, the incident field frequency is 

-3 

1000 hz, the conductivity of the whole-space is 10 mhos/m and the con- 
ductivity of the cylinder is 10"^ mhos/m. 
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It Is seen in Figs. 52 and 53 that except for imaginary (M g ) , an 
accuracy to about 1% has been obtained with 40 sampled values of the 

*S 

equivalent current densities. However, an accuracy to only about 5% 
has been obtained for imaginary (M s ) and, as a result, the fields cal- 
culated from these currents have an accuracy to about 3% with 40 sampled 
values . 

Fig. 56 indicates that much better convergence is obtained if the 

-2 -1 

conductivity of the cylinder is increased from 10 mhos/m to 10 mhos/m. 
For comparison, this example has been calculated also using the integral 
representations derived for conductors with a low surface impedance and 
a | small curvature. It is evident that in this example the general integral 

representations are required to predict imaginary (1^) accurately. 

-2 

As the conductivity of the cylinder is decreased below 10 mhos/m 
(assuming that all other parameters of Figs. 52 through 55 remain fixed), 
the accuracy and convergence of the solution become increasingly unsatis- 
factory. In Fig. 57, the horizontal component of the magnetic field 
intensity scattered by a cylinder with a conductivity of 0 has been 
plotted for 20, 30 and 40 sampled values of the equivalent surface current 
densities. Although an accuracy to less than 3% is attained m estimating 
the equivalent surface current densities, a surprisingly large error of 
14% to 30% results for real (H ) . 

X. 

Consequently, the scattered magnetic field intensity was calculated 
using analytical values of the equivalent current densities on the surface 
of a circular cylinder with a conductivity of 0, Since an accuracy to 
better than 3% was attained, this indicated that no programming errors 


existed which had not been detected in earlier tests. 
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However , the possibility existed also that these large errors might 
be due to the fact that the numerical cylinder was a polygon rather than 
a circular cylinder. To test this likelihood, a twenty-sided polygon 

i 

with a conductivity of 0 was sampled 42 times. It is evident in Fig. 57 
that the results closely duplicate the data for a forty-sided polygon, 
suggesting that the errors are not dependent upon curvature. 

Nonetheless, Neureuther (1969) believqs that this test still does 
not exclude the possibility that the solution is curvature dependent. 

The reason for this is that by sampling a twenty-sided polygon 42 times, 
the sampling points have been shifted closer to each corner of the polygon 
where' the influence of the polygon curvature will be more important. As 
a result, it will be necessary to solve the integral equations assuming 
that a circular cylinder has been sampled 20 times rather than assuming 
that a twenty-sided polygon has been sampled 20 times before the impor- 
tance of curvature can be 'resolved. 

However, it would appear that the accuracy could be limited also 
by requiring a- very accurate knowledge of the equivalent surface 
current densities whenever the reflection coefficient is small. The 
reason for this can be understood best by examining Figs. 58 through 60 
in which the problem of a cylinder having the same electrical parameters 
as the surrounding whole-space has been considered. Since no field is 
scattered, the tangential surface current densities are equal to the 
tangential components of the incident field, Ky = H^, M g = -E^, and 
equation (5-24) reduces to 
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(5-49) 


This result, which always must be true for any point of observation 
exterior to the cylinder, is deduced mathematically in section (6-1). 

As a result, it is evident- that the accuracy obtained in predicting 
the scattered fields is directly dependent upon an accurate knowledge of 
the scattered fields around the boundary of the inhomogeneity, since in 
general 
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Thus, it is expected that the accuracy and convergence of the general 
integral representations will be poor for those problems m which the 
reflection coefficient is small. 

-2 

In Figs. 58 through 60, a maximum error of 8 x 10 Araps/m (z. 2.%) 

-2 

in imaginary (M ) yields a maximum error in real (H ) of 3.5 x 10 Amps/m. 
s x 

_2 

Similarly, for N = 30 in Fig. 57, a maximum error of 7 x 10 Amps/m 

-2 ' 

(1.87 0 ) in imaginary (M ) yields a maximum error of 3.0 x 10 Amps/m 

s 

(197) in real (H^) . Since the scattered field is small with respect to 
the incident field at the contour (|E st ’|/ 1 E L | ••w .043 and )R 5 C |/]h 1 | .12 

for & = 0 and, |e SC | /{e 1 ) ^ .089 and ^ .16 for & =1Y), a 

significant percentage error is obtained in calculating the scattered 
fields even though the equivalent surface current densities are known to 
an accuracy of 27 

This conclusion is substantiated by the results of Fig. 61. In 

this example we have assumed normal incidence of an E^-polanzed plane 

wave and that: the depth z, to the top of the cylinder is 20 m, the 

cylinder radius is 100 m, the incident field frequency is 1000 hz, the 

conductivity of the whole space is 0, and the conductivity of the cylinder 
-3 

is 10 mhos/m. It is seen that an accuracy to about 2.57 in iH^l 

is achieved since the scattered magnetic field (which is more important 

than Ey in this example) now is much larger than the incident field 
sc i sc 1 

()E 1 / | E y |*» .24 and Ir^ l/jH s |/v 17.5 f or Q = 0.). It should be noted^ 

also, that the general integral representations are required to compute 
this example since the curvature approximation would have been violated 
had the solutions of Chapter 4 been used (compare Figs. 45 and 61). 

To illustrate the application of the method to problems which cannot 
be handled analytically, the field scattered by the vertical slab of 
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Fig. 24 has been plotted in Figs. 6'2 and 63. It is assumed that an Ey- 
pol-arized plane wave is incident normal to the slab and that: the depth 

to the top of the slab is 20 m, the incident field frequency is 1000 hz, 

-3 

the conductivity of the whole-space is 10 mhos/m, and the conductivity 
of the slab is 10 mhos/m. 

Since only a 1% change is obtained in the scattered magnetic 
field intensities as the number of sampled values of the equivalent 
current densities is increased from 30 to 62, we have assumed that a 
convergent solution has been found. This fact is suggested also by Fig. 

56 in which a cylinder having the same electrical parameters and about 
the same contour length as the slab of Figs. 62 and 63 has an accuracy to 
about 1% for 30 sampled values . 

It is interesting to note that at this frequency and with these elec- 
trical parameters , the inflexions of imaginary (H) are indicative of 
Cylinder width. It is evident from Fig. 62 that the slab brings about 
a significant decrease in imaginary (1^) between x = ± 50 m whereas 
the circular cylinder of Fig. 56 has a flat response between x = ± 50 m. 
(Actually, there is a very slight decrease in imaginary (H^.) over the 
circular cylinder.) The reason for this behaviour is that we are 
observing each corner of the slab .(where the radius of curvature is 
small) respond to the incident field. This will be discussed more 
thoroughly in section (7-3) where it will be seen that the corners 
mainly affect imaginary (K^) . As a result, the effects of corners are 
apparent in imaginary (H) . 
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CHAPTER 6 

TWO-DIMENSIONAL AND THREE-DIMENSIONAL INCIDENT FIELDS 

In Chapters 3 , 4 and 5, an integral equation solution to scattering 
from cylinders in a conductive whole-space was examined. In the derivation 
of these integral equations, however, it was assumed implicitly that the 
incident field was that of a plane wave. 

We will show now that this assumption is unnecessarily restrictive 
and that the integral representations are valid for any two-dimensional 
source configuration. In fact, any three-dimensional source configuration 
can be considered by expanding the primary current distribution and the 
field it radiates into a Fourier integral over a continuous mode distrib- 
ution and solving this two-dimensional problem. 

6-1 Derivation of the Integral Representations Assuming Any Two-Dimensional 

Source Configuration 

If the incident field is other than that of a plane wave, it is 
invalid to write that the axial components of the total electric and 
magnetic field intensities satisfy the homogeneous Helmholtz equation in 
the exterior region (assuming that this is the region which contains the 
sources) . However, since the scattered field intensities always satisfy 
the homogeneous Helmholtz equation, it is valid to write that in the exterior 
region 


(v\kM E^(j5) = 


) 


o 


(6-1) 
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and 

i^t 2 ) Uj'cpi = o. 

(6-2) 

The exterior Green's function remains unchanged and must satisfy an 
inhomogeneous Helmholtz equation given hy 


. -Hf-fl. (# . 3) 

Thus, it follows from section (2-1) that if equations (6-1) and (6-3) 
and equations (6-2) and (6-3) are introduced separately into equation (2-2), 
then 



since the contribution from the outer contour C tends to zero as C 

2 2 

approaches infinity. 

However, the incident field does satisfy the homogeneous Helmholtz 
equation in the interior region since by assumption this domain is source 
free. Thus, 


'f) Ej (jrU o 

and 

c^+k x ) C j) = 0 


(6-6) 


(6-7) 


in the interior region. In addition, the exterior" Green 's function is 
defined in the interior region by equation (6-3), 
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As above, it follows from section (2-1) that if equations (6-6) and 
(6-3) and equations (6-7) and (6-3) are introduced separately into equation 
(2-2) for the interior region, then 


0 * 


f G (j^') Or') - £y ( .j ') ^ 

>J(2 J 


and 


0 


* ( f if) {?') - Ujtf) 

J<1, ^n.' dh J 


( 6 - 8 ) 


(6-9) 


since the Green's function possesses no singularities m the volume 
integral inasmuch as ^ refers to any point in the exterior region. 

Adding equations (6-8) to (6-4) and (6-9) to (6-5), it is found that 


and 


X t L V ** J 

= f ( c * c FiP ^ L P - 2 ^ 

Jc l 1 M J 


(6-10) 


( 6 - 11 ) 


Upon adding the incident field intensities to equation (6-10) and 
(6-11)’, the desired integral expressions are obtained: 


£y(f) 


e j L f) + f l (p ~ Ey C pM- iff' ) 

Jc L da' dn.* 


(6-12) 


and 


Hj if) - ) *- 




1 if, ft if) - if ijif'Hh'- ( 6 - 13 ) 

v <w J 



125 


It is evident on comparing equations (6-12) and (6-10) with equations 

(2-9) and (2-10), and equations (6-13) and (6-11) with equations (2-12) 

t t 

and (2-13), that the integral expressions for E y (^) and H y (^ ) are unchanged 
Since the integral representations of the fields in the interior region 
remain unaltered for two-dimensional sources, it follows that the general 
field representations given in sections (2-3) and (2-4) are valid for any 
two-dimensional source distribution. 

6-2 Numerical Examples: Electric Line Source 

The electric field intensity at a point of observation p radiated 
by an electric line source located at is given by equation (3-6) as 

- - At# (■& p-fi. ; ) > (6_14) 

A 

where I is the electric current flowing in the source. 

For simplicity, we will confine our examples to scattering from 
perfectly conducting cylinders, in which case the desired integral equation 
is found to be (from equation (3-8)) 

. - f ■ <6 - 15) 

Jd, 

Once equation (6-15) has been solved for K y (|? / ), th e scattered transverse 
magnetic field intensity is given by equation (3-12) as 

H fe (f)* A f tf/ U*i ds - ( 6 - 16 ) 

A J <y 

To demonstrate that the integral representations are valid for two- 
dimensional source configurations, we will compare the numerical results 
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with the analytical results for the case of scattering from circular 
cylinders. It can be shown by following a development similar to that of 
Harrington (1961, p. 236) that the field scattered by a perfectly conducting 
circular cylinder in the presence of an electric line source, which is 
parallel to the cylinder is given by 


eft *#). < 6 - 17 > 

A 


A 




and ' 


4> 


H** C\$) * 1 1 2 ^ m ; 

A 1 — 


Co* 


(6-18a) 


(6-18b) 


where 


/H 


- 4- G, R) , 




and 


r ‘* 4 


a= e 
ft 2> * 

7 


(6-19a) 
(6 -19b) 


The coordinate system is defined in Fig. 64. 

It is interesting to note that the "reflection coefficient" a n is 
the same as that obtained for the case of an Ey-polarized plane wave 
(equation (3-48)) scattered by a perfectly conducting circular cylinder. 
Harrington (1961, p. 237) notes that in general the "reflection coefficient" 
is independent of the incident field. Thus, we can write immediately from 
equations (4-15) and (4-16) that the field scattered by a finitely conducting 
circular cylinder in the- presence of an electric line source which is 
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parallel to the cylinder axis is given by equations (6-17) and (6-18) with 


and 


s 


W -Cn.'C(l?,R) 


(6- 20a) 


C* - 


e, -U (1^*0 


(6-20b) 


where the derivatives in (6-20) are with respect to the argument 
(kR). 


This solution is similar to that given by Wait (1952). 

Figs. 65 through 67 demonstrate that the integral representations are 
valid for two-dimensional source configurations and that the solution 
converges as the number of sampled values of the current density is increased. 
We have assumed in this example that an E -polarized electric line source 
is located 20 m above the top of a perfectly conducting cylinder and that: 
the cylinder radius is 100 m, the incident field frequency is 1000 hz, and 

_ O 

the conductivity of the whole-space is 10 mhos/m. 

It is seen that a larger number of sampled values of the current 
density is required to obtain a 1% error than when the incident field 
is a plane wave (compare Figs. 65 through 67 with Figs. 11 through 16). The 
reason for this is that the field of a line source falls off much more 
rapidly than that of a plane wave. As a result, a small sampling interval 
is required for line source scattering problems to describe accurately the 
rapid current variation around the cylinder. 

However, it is evident from Fig. 65 that a small sampling interval is 
required only on the upper portion of the cylinder since this is where the 
largest current densities are located. At an angular distance of from 
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the top of the cylinder, the surface current density is 1/20 of the maximum 
value whereas at an angular distance of { 2, the surface current density 
is 1/100 of the maximum value. Thus, a small sampling interval is required 
only between = ± 7f [ 4. 

With this fact in mind, the cylinder was sampled 26 times using the 
sampling distribution of Fig. 68. Intervals i through 8 and their mirror 
image about the z-axis have the same sampling distribution as used in 
Figs. 65 through 67 for N equal' to 50. Then, the cylinder was sampled 3 
times to a & of/T/2 and the lower half of the cylinder was sampled only 
5 times . 


A maximum change of 0.1% from the results of Figs. 66 and 67 was 
obtained m the horizontal interval (-100,100). The changes were largest 
at x = ±100, since the point of observation is closer to the less accurate 
sampling here than at x = 0. However, the fields close to x = 0 are of most 
interest to the field geophysicist since this is the region of maximum 
field- intensities . Here, the changes are less than 0.01% from the results 
of Figs. 66 and 67. Thus, by an appropriate choice of interval distribution, 
it seems possible to obtain about the same numerical accuracy independent 
of the incident field using the same number of sampled values of the 
surface current density. 

To illustrate the application of the method to problems which cannot 
be handled analytically, the field scattered by a perfectly conducting 
vertical slab (Fig. 24) has been plotted in Figs. 69 and 70. In this example. 
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The slab was sampled 30 and 42 times using the general sampling 
distribution of Fig. 24. Than, choosing a sampling distribution with a 
line source in mind, the slab was sampled twice as often as for N = 42 
to a depth of 125 m on the slab (the depth at which the current was 1% 
of the peak value), and only 9 times on the rest of the contour. Since 
a change of less than 17, m the predicted field occurred when the contour 
was sampled twice as frequently (effectively), we can assume that a 
convergent solution has been found. In addition, F igs. 66 and 67 suggest 
that m a similar problem, an accuracy of greater than 1% can be obtained 
if a smooth contour Is sampled 500 times per wavelength. Since this sampling 
rate is true for N = 50 on the slab contour above a depth of 125 m, this 
result suggests also that a convergent result has been obtained for N = 50 . 


6-3 Integral Representations Assuming Any Three-Dimensional Source 
Configuration 

Harrington (1961, p. 292) notes that ”a three-dimensional problem 
having cylindrical boundaries can be reduced to a two-dimensional problem 
by applying a Fourier transformation with respect to y (the cylinder 
axis)." For three-dimensional source configurations, we seek a solution to 
the three dimensional Helmholtz equation 


(jL +1. + 

\ 3. 1 V Si 1 ' j 


Thus 


•iL 


)-OtS 

will be a solution to the two-dimensional Helmholt 2 equation 

[L * r ) 


(6-21) 


( 6 - 22 ) 


(6-23) 
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where X 


We can write immediately the solution to equation (6-23) from equation (6-10) 
as 





h 


/ 

> 


(6-24) 


where ^ ^ ,}') ’ M,/ C^l p_p' l) , 

or, on applying a Fourier transformation to the incident field with 
respect to y 








(6-26) 


Equation (6-26) asserts that the y dependence of (x,y, z) (and 

all other transformed quantities) is of the form e^ * v y ^ . Since the inte- 
gral representations derived in sections (2-3) and (2-4) are for just such 
an axial dependence, equation (6-25), can be solved once k is specified. 

y 

Thus, it is evident that a three-dimensional problem having cylin- 
drical boundaries can be solved by expanding the three-dimensional incident 
fields into a Fourier integral over a continuous mode distribution and 
solving equation (6-25) for a number of discrete values of ky. Conse- 
quently, the three-dimensional solution is obtained by fitting a. curve 
to a weighted sum of sampled values of c (^ky) and performing (6-26) 
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numerically. 

If the incident field is a plane wave in free space incident obliquely 
to the cylinder axis, then we can write that 

^ +• Zcoi 6 tos.^fa ) 

> ( 6 - 27 ) 

where is defined in Fig. 2 , 

and 6 is the angle in an axial plane between the direction 

of propagation and the plane normal to the axis of the cylinder. 
The transformed incident field is obtained from equation (6-22) and is 
found to be 

. , . 1^0 Usih <&* + ZCoi<Pi«) 

3 37T.^> \\ oJ i - £((^5 £h ) (6-28) 


where 


to = k«> * / \C~ . 


Equation (6-28) indicates that whenever the incident field varies har- 
monically as along the cylinder axis, equation (6-26) need be solved 

i 

for one value of only. In the case of a plane wave in free space 
incident obliquely to the cylinder axis, ky is given by 


fey a Ko Swt ^ • 


(6-29) 


It should be noted, however, that this value of k is to be used solely 
m the exterior integral representations. For interior integral represen- 
tations, k is determined m each region from Snell's Law. 

y 

If the incident field is a plane wave in a conductive whole-space 
incident obliquely to the cylinder axis, then equation (6-28) is invalid 
and instead, equation (6-25) must be solved for a number of discrete values 
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of k . Unfortunately, it is seen from the physics of this situation that 
the problem does not possess a solution since the plane wave amplifies 
xn the negative y direction. 

Nevertheless, xt xs not surprisxng that a solution does not exist since 
the problem is unrealistic in that it requxres a source with an infxnite 
power supply at infinity. This complxcatxon has been ignored in earlier 
chapters by assuming that a plane wave incident normal to a conductive 
half-space exists at the earth-air interface and attenuates from this 
point as it propagates into the earth. Essentxally, it was sufficient 
to assume that a plane wave generator existed at the earth-air interface 
since this closely approximated the actual situation. 

If oblique incidence is required, it will be necessary to assume 
that a plane wave generator exists at (x, o, z) and that the wave attenuates 
in the positive and negative y directions from this "source”. However, 
obliquely incident plane waves propagating in a conductive whole-space 
are not important in most problems of geophysical interest. The reason 
for this is that the transmitted field at the earth-air interface propa- 
gates normally away from the boundary even for grazing angles of incidence. 
Thus, when discussing plane wave scattering by cylinders m a conductive 
whole-space, it is adequate to consider normal incidence. 

These problems do not arise when a plane wave is incident obliquely 
to a cylinder in a conductive half-space. In this case, the plane wave 
is harmonic m the exterior region (air) and the problem need be solved for 
one value of 1c, only. 



CHAPTER 7 


SCATTERING FROM CYLINDERS IN A CONDUCTIVE HALF -SPACE 


In Chapters 3, 4 and 5, integral equation solutions to scattering 
'from cylinders and from topography were examined. It was seen that any 
problem could be formulated by choosing an appropriate integral repre- 
sentation in each region present and a solution obtained by solving the 
resulting integral equations. 

Having established the validity of the integral representations and 
the subroutines written to effect their solution, it is possible now to 
combine the cylinder and topographic results to study scattering from 
cylinders of arbitrary cross section in a conductive half-space with an 
arbitrary earth-air profile. Although this chapter will be concerned 
with normal incidence of an Ey-polarized plane wave, the results of 
Chapter 6 can be used to extend the analysis to sources of arbitrary 
configuration and polarization. 


7-1 Derivation of the Integral Equations 


If it is assumed that the primary field is a normally incident Ey- 
polarized wave (ky =0), then it follows from equations (2-56) and (2-66) 
that the integral representations for E and H in the region exterior to 


the half-space are 

= fcy'tf)- jj o co f Lf') H 

+ ike f Cel(ji -ot). > 

i J, ' ‘ 


(7-1) 


133 



134 


and 


H* if) = (j) +i L f W?, CL'rfW 

1 Jc. 


4ZoJ 


where 




(7-3a) 


and 


(fe>^pO - Kp > CU^-^'()(sh» ) 


(7-3b) 


** ( j> C 5w ( 3^jS - •*) y)c — Ceri (3^> -«*) * 

~ i-if-f'l 

The parameters of the exterior region have been referred to by the subscript 
0 and the equivalent surface current densities on contour 1 have been 
superscripted 1. The nomenclature for the contours and regions present 
is shown in Tig. 71. 

As in Chapter k, it is assumed that the half-space contour is flat 
outside an interval bounded by (-a, a) and that outside (-a, a) the equi- 
valent surface current densities arise from fields reflected by the half- 
space alone. Thus, equations (7-1) can be rewritten from equations (4-22) 
and (4-37) as 



, d- 4 






- <*- S' 
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and equation (7-2) can be rewritten from equation (4-45) as 


*n t 


'V * 0 = [ f } + j 


i «-s; 


- 4 - 5 / 




(7-5) 


where is the x coordinate of the point of observation, 
reflec 0 

Ey ( jj ) is the electric field intensity reflected 

by a flat conductive half-space and is to be used when 
the point of observation is above or on the flat half- 
space, 

-E^( p " ) is the incident electric field intensity and is to 

be used when the point of observation is below the flat 

half-space, 

reflec 

H fcr ( |s '*■) is the magnetic field intensity reflected 

by a flat conductive half-space and is to be used when the 
point of observation is above or on the flat half -space, 

l 

-H tr ( ^ ~ ) is the incident magnetic field intensity and is to 
be used when the point of observation is below the flat 
half-space, 

F f( K, p t p') and F^( Iq, ( p r') are given by equations (7-3) 
j 1 A f-'i u«h 


r e * /tw f * i<j Lf') 0 

4 J . 




- it Lf) cm t 


t 


“-'I 


(7-6) 
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and 
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-*-* c 
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P*^*h 


W S C^ic^rt (j6-x) 14* Cklj-p'l) Lz* i j!,& -c«j3« )is 

(7-7) 
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tj*') H, (klf-f'f) L$i*~Cl.p - «*V£ - Cro £^yi- *.)£ } , 

Jll/y'l 


-«.- St 


It should be noted that >as in Chapter 4, the contour integrals involving 

(Ky i are along a flat half-space while the contour integrals involving 

11 

(Ky s ) are along the topographic profile. 

The integral representations for E and H in region 1 do not follow 
directly from Chapter 4 since it is not evident what result the infinite 
integrals involving (Ky , ) will yield. To establish the value of these 

integrals, we will examine the electric field intensity transmitted into 
a homogeneous earth assuming normal incidence of an E^-polarized wave 

(ky f 0) * 7n this, case, we can write from equation (2-72) that 

\ 

4 J-no 4 • 

** (7-8) 

As in the exterior region, it is assumed that the half-space contour is 
flat outside an interval bounded by (-a, a) and that outside (-a, a) the 
equivalent surface current densities arise from fields reflected by the 
half-space alone. Thus, equation (7-8) can be rewritten as 


E 3i^'*/4^2 


-a Si OO . . y 

" | K* If') H? tfa.ljr-p'l M. (t.lpj'l) 4 «'j 


CL- S,- 
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- ck, | Cp') ora(j3--i) rf* C^i jji-f'OdL' + (p')aa(p -*) J 

^ - J-«5 4t-s c 

x S' ^ (7-9) 

4_ /x, u> (pQ _ X f B$(j'J c« (^3 -°0 ^(.fq 1^-^' •) «^ s "• 

*4 J . -4 


By rewriting the infinite integrals of (7-9) as 
o& oo 

^ K^(p') - ik w s(f') °* 

4 - —•*. ^ - _ «Q 


(7-10) 


" A ^ J C^r') H a tkj ji-j'i)^ -f ik, Mj Cj>') Usi H| G,|p^'l) JU } 


-*-S f 


it follows from (7-8) that for ^ below the half-space, the integrals 
over ( _ oo co ) must equal the total electric field intensity transmitted 
into a homogeneous earth by the flat half -space. Consequently, equation 
(7-9) reduces to 

S. ( p> • % Cf 4 ) - r e 

4 \-o-- K 

If j> is above the half-space, the integrals over (_<o ( £0 ) must 
equal zero. The reason for this is that the integral representations pre- 
dict a null field for points of observation outside the volume of interest 
since the boundary conditions have been satisfied by the equivalent surface 
current densities (see Harrington, 1961, p. 106). 

Thus, in general the fields inside a conductive half -space can be 
written as 


* - Ms (p') onCp- <*-) tfKklf-f '!)&$'. 


(7-11) 
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+ /is£ 

4 J 


3--S l 


-*"S c 


a-s: 






(7-12) 


where 0( p + ) designates that this term is to be set equal to 

zero when the point of observation is above the flat half- 

space, 

trans 

and Ey ( jj ) is the electric field intensity transmitted 

by a flat half-space and is to be used when the point of obser- 
vation is below or on the flat half-space. 

It can be shown analytically that the infinite integrals of equation 
(7-10) yield the transmitted electric field intensity for p below the 
half-space and zero for J above the half -space if it is assumed that the 
incident field is an E^-polarized plane wave impinging normal to a flat 
half-space. In this case, is constant and is given from equation (4-21) 
as 


*3 * 


.5 2, 




^ + 2* 


(7-13a) 


Similarly, is constant and is found from (2-39b) and (4-25) to be 


K 


- (e^+eJ JU c ) = 


- H, 


(7-13b) 


+ z i 


Introducing equations (7-13) into (7-8) and assuming that p is below 
the half-space, then equation (7-8) becomes 


E*. (f~3 a f Mf + i&( H?* Cfc, y 

5(2 0 +2 f )LI 0 o -L /sVi' 4 


(7-14) 
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Using the results of equations (4-28) and (4-29) , then (7-14) 
becomes 






(7-15) 


On comparing equation (7-15) with (5-22), it is evident that the infinite 
integrals of equation (7-10) yield tke transmitted electric field inten- 
sity when p is below the half-space. If jj is above the half-space, 
(7-14) becomes 






ii f ill" (1/5*77^ 1 • 
L i/sV*'* < 7 - 16 > 


Introducing the results of equations (4-28) and (4-29) into (7-16)^ it 
is found that E^( p + ) is zero. 

A similar analysis can be followed for ) in region 1 begin- 

mng with equation (2-~74), or by computing H from (7-12) using Maxwell's 
equations. It would be found that 



where 0 (J + ) designates that this term is to be set equal to zero 

I 

when the point of observation is above the half-space, 
trans 6 . 

H tr ( p » ) is the transverse magnetic field intensity trans- 

mitted by a flat half-space and is to be used when' the point 
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of observation is below or on the half-space, 

) and ^ .f 7 > are siven by e <l uations (7-3), 

and ) is given by equation (7-7). 

When a cylinder is present in the ground, equations (7-12) and 
(7-17) must include a contour integral around the cylinder boundary. 
Thus, in the present problem, (7-12) and (7-17) can be rewritten as 




The parameters of the half-space have been referred to by the subscript 
1 and the equivalent surface current densities on contour 2 have been 
superscripted 2 (see Fig. 71). 

It should be noted that the integral representations for an interior 
region, equations (2-72) through (2-75), were derived assuming that the 
normals on the contour are into the interior region. Since the normal 
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on C 2 has been chosen as into the inhoraogeneity in Fig. 71, the opposite 
sign to that of equations (2-72) and (2-74) has been taken for all inte- 
grals in region 1 around C 2 « 

Finally, the integral representations for E and H in region 2 
follow directly from equations (2-72) and (2-74) . Assuming normal inci- 
dence of an Ey-polarized field, it is found that 


- 




4 4 






and 


H i -iL- f tk,f + k_(’ M s (7-21) 

where F^ ( ' ) and F^ ( ^ ^ f are .given by equation (7-3), 

and the parameters of the inhomogeneity have been referred to by 
the subscript 2. 

I The unknown current components Ky( jT ) and M g ( p ) on contours 1 and 

« «* r 

2 are obtained by enforcing the boundary conditions on tangential E and 
H and solving the resulting integral equations. We have from equations 
(2-53) that the tangential unit vector s is given by 




d ' /X. + Coi o£? £ 

L c > 


(7-22) 


where d.- L is the normal at the boundary point. 

Thus, the tangential component of in each region is found by 
taking the dot product of (7-10) with equations (7-5), (7-17) and (7-19) . 


It is found that 
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and 
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If we let the point of observation approach a boundary point on C.. 


given by the position vector jj - “ and equate tangential components of 
E and H, it is found that 




e / i n 


°<j P"‘) { - 


JU* f K^CjroFw?^lf' -p' 1 ) +/*i 


(7-28) 
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Similarly ^ if we let the point* of observation approach a boundary 
point on C£ given by the position vector p "* and equate tangential 
components of E and H_, it is found that 



and 


0 = 


0(p"'*) 


£ 
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A -M i W U"' v ) 
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Equations (7-28) through (7-31) are the desired coupled integral 



are solved as in Chapter 5 by dividing (-a, a) and into a total of N 
sampled values of each unknown current component. "When this is done, the 
integral equations reduce to a system of linear equations of the form 


t e )' = (<,)( u s,; x §) + ( x Hj) + c ) 

(m' ■ ( a ;, x k; ) f * (niXMi)+ 

Co)' - + CBJWS*) + W,\X Hi)* C^JCHJ), 

(°V = ( ) + Cb^Xh) ) . 


The superscript on A^j , } refers to the boundary upon which the boun- 

dary condition is being enforced when each matrix is determined. The 

primes on (E^), (H^) and (0) indicate that these column matrices are 

-i 1 

modified due to the contour integrals involving (K , M ) on contour C . 

s s \ 

Equations (7-32) can be rearranged into a system of 2N linear 
equations of the form of equations (5-47). It is found that 



(7~33a) 


(7-33b) 


Equations (7-33) are solved in exactly the same manner as equations (5-47) . 

Their solution is given by equations (5-48) by treating each partitioned 

matrix (A., . B, 1 ) as C 

' ij ij 

Once K^( ^ ) and ( S' )* have been estimated from equations (7-28) 
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through (7-31) , the total or scattered magnetic field intensity is calcu- 
lated from equation (7-5), (7-17) or (7-19), depending upon which region 
the field is desired. Exterior to the half-space, the scattered 


magnetic field is given by 
-*"■ refUc 


- j v;v- + > j - visf.fO 

- H h ~) J 
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where 1^ ( | , f ' ) is given by equation (7-7). 

It should be noted that as in Chapters 3, 4 and 5, it is necessary 
for integration accuracy to make a small argument approximation to the 
Hankel functions whenever lkilp-f't<,5 and integrate the resulting ex- 
pressions analytically. However, since the integral representations of 
this chapter are the same as those of Chapter 5, this problem has been 
discussed previously in section (5-1) (equations (5-30) through (5-33)). 
Similarly, integration through the point of singularity has been examined 
in section (5-2) . 

It is interesting to note that equations (7-30) and (7-31) reduce to 
(5-20) and (5-21), the integral equations for scattering from cylinders in 
a conductive whole-space, if we assume that the earth-air interface is 
horizontal and ignore coupling* between the cylinder and the interface. 


* The term coupling will be used to describe the electromagnetic inter- 
action between the cylinder and half-space contours which modifies the 
equivalent surface current densities on the cylinder surface from their 
whole-space value. 
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In this case, K'(^) and M*( ) arise from the incident field alone and, 

as a result, (-a, a) can be set equal to (-0., 0.). Thus, equations 
(7-30) and (7-31) reduce to 






H 0 l) 



(7-35) 



and 



Assuming that the cylinder is below the half-space, it is evident on 
comparing equations (7-35) and (7-36) with (5-20) and (5-21) chat these 
are the same integral equations as derived to examine scattering from 
cylinders in' a conductive whole space. In' (7-35) and (7-36), the inci- 
dent field is the field transmitted into a flat, homogeneous half -space. 


7-2 Numerical Examples 


(7-31) 

Equations (7-28), (7-29), (7-30)^and (7-34) were programmed assuming 
that the incident field is an E -polarized plane wave incident normal to 

y 



148 


i i 

a conductive half-space. In this case, and M g are constant and given 
by (7-13). Furthermore, the integrals 1^ and I reduce to expressions 
which are similar to equations (4-39) and (4-48) . Thus, we can rewrite 
(7-4) and (7-5) as 




i! -4-5 : 


+S* J*s 


f i ik£ Hrcfc/s?) ^; (7 _ 37) 
i i AS j ' 1 


J- H (.t,*,*') = Ky ( c> T, ^ 2 ) 4. kHs ( & £5- - k I, - . 

A AB L > 

+ & / X*, - a T l0 \ e 1 

1 k J J 


(7-38) 


where 1^ , , Ig , 1^ and I are' given by equations (3-63a), 

(3-63b), (4-13a) , (4-13b), (4- 14a), (4- 14b) and (4-l4c) res- 
pectively, with a = -a-s^ and b = a-s^. 

I When the point of observation ( jp ) is on the half-space, a limiting 
process similar to integrating through the singular interval (section 5-2) 
must be taken before (7-37) and (7-38) can be evaluated. Noting that this 
limit is taken from opposite sides of the interface in each region, it is 
found that for points o£ observation in air, equations (7-37) and (7-38) 
reduce to 


a jO = JU 0 U> Kj { HKoO ? + Ms] 

4 ^ J 2 
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and 



Kj / \ .£ + t ko ( X 2 1 

[*)- Mi iu f^V r,,)^ 

(7-39b) 

W 4 U->* - 
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Similarly, for 

points of observation 

in the ground, 



//., u) |<^ JU { HKoo 

}-? * 

(7-40a) 

and 
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If follows from section (5-2) and Appendix D that 

§~*.o v. J S-*“0 J 

J -*-K 

= A z (f*-^ 3 ) + i]l 
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where 

Equation (7-41) states that HKOO is approximated analytically over (-e^£ a ) 
by taking the limit as S »*■ O of the small argument solution given 
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m Appendix D, and numerically over the intervals (-a-s^ , - ) and 

(g , a-s.). The subscripts on £ imply that if £ p > I -a-s.'j or 
8 > (a-s ), - En or is to be replaced by -a-s. or a-s. respectively, 
and the first or second integral of (7-41) ignored. 

It is shown in Appendix D that can be evaluated analytically, and 
it is found that 
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^ — *>- o 
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finally, it follows from section (5-2) and Appendix D that 
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As in equation (7-41), or ^ is to be replaced by -a-s^ or a-s^ 
respectively, if £ > | -a-sjor <f > (a-s^) and the first or second 
integral of (7-43) ignored. 
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Assuming that the incident field is an E^-polarized plane wave im- 
pinging normal to a conductive half -space, then 


Ej‘ If") 




^ If") = -sCvtod- h o ui 


_ iU 


The reflected field intensities are given by equations (4-24) and 


(4-25) as 


_ I'e'ff-** 
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Finally, the transmitted field intensities are given from equations 
(5-22) and (5-23) *ls 
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It should be remembered, thatl equations (4-40) through (4-44) 
must be considered at the edges of the interval (-a, a) to ensure con- 
tinuity of the unknown current distributions with K 1 and M L . 

y ' s 

Figs. 72 through 79 demonstrate that the integral representations 
derived for scattering from cylinders in a conductive half-space are valid 
and that the solution converges rapidly for these electrical parameters as 
the sampling rate is increased. In Figs. 72 through 78, we have assumed 
that the conductivity of the air is equal to the conductivity of the 

_3 

half-space (10 mhos/m), so that the solution must reduce to the solution 

i 

for scattering from circular cylinders in a conductive whole-space. It 
should be noted, however, that since the entire integral equation must be 
solved (i.e., no expressions in equations (7-28) through (7-31) reduce to 
zero when this assumption is made), such a test forms an accurate check on 
the formulation of the problem for scattering from cylinders in a conductive 
half-space. 

In addition, we have assumed normal incidence of an E -polarized 

' y 

plane wave in Figs. 72 through 78 and that: the height Z Q of the point 

of observation above the half-space is 1 m, the depth Z^, from the half- 
space to the top of the cylinder is 19 m, the cylinder radius is 100 m, 

the incident field frequency is 1000 hz, the conductivity of air is 

—3 —3 

10 mhos/m, the conductivity V; of the half-space is 10 mhos/m, and 

-1 

the conductivity of the cylinder is 10 rahos/m. 

Fig. 72 indicates that with these electrical parameters a contour 
section bounded by (-lOOOj-f- 1000) is sufficient to describe the unknown 
surface current densities. Fig. 73 shows that an error of less than TL 
in the peak value of real (H x ) is achieved if the half -space profile is 
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sampled 20 times. It is seen also that a slight instability arises in 
the field, although the amplitude of the oscillations decreases as the 
sampling width on the topographic profile decreases. 

As in the examples in Chapter 3, these oscillations arise from the 
inaccuracy of the numerical integration across intervals of the contour 
m close proximity to the point of observation. It was noted m section 
(4-3) that the integration inaccuracy is a result of the manner in which 
the program has been written since the small argument approximation is 
used only when |kt'| < Re for all points on an interval. As a result, 
it was found that the numerical integration was inaccurate across any 
interval in which part of the contour is greater than Re while most of the 
contour is much less than Rc. 

It was seen from Table 1 that for such an interval, solution accuracy 
could be increased by raising the number of sampling points within the 
interval (-a, a) or by increasing Re. Note, however, that Re cannot 
exceed the validity of the small argument expansion for the Hankel func- 
tions. Figs. 73 and 74 indicate that these oscillations are a result of 
integration inaccuracies since their magnitude is quite dependent upon 
the sampling rate (Fig. 73) and the value of Re (Fig. 74). 

To overcome this problem, the program should be rewritten so that 
both the small argument approximation and Simpson's rule can be used in 
integrating whenever Re falls within an interval. Since these oscillations 
amount to less that 17 a of the total field with an Re of .5 and careful 
sampling, the present program has not been changed. However, before a 
general analysis of scattering problems is undertaken, this change should 
be made to avoid an unrecognized error. 
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Fig. 75 shows that by sampling the topographic contour 29 times and 

the cylinder contour 30 times, an overall error less than 2% of the peak 

value of H has been achieved for the horizontal magnetic field intensity, 
x 

Similarly, Fig. 76 shows that with the exception of several discrete 

points, an error of less than 2% of the peak value of H z has been achieved 

for the vertical magnetic field intensity as well. 

The discontinuities in H are observed only when the point of obser- 

z 

vation is near the contour and directly over the edge of two intervals. 

In this case, the contribution to the fields comes mainly from that part 
of the contour in the immediate vicinity of the point of observation. 

As a result, the discontinuity in the parabolic fit to the current density 
(see Fig. 7) predominates -and a poor estimate of the scattered field is 
obtained. 

Further observations which point to the break in the current density 
as the origin of these discontinuities in H z are: 

1) The magnitude of the error in H decreases (Fig. 76) as 

z 

the parabolic fit to the surface current density improves 
(i.e., as the number of sampled values on the topographic 
contour increases) . 

2) The magnitude of the error is larger for points of obser- 
vations near the cylinder (i.e., the parabolic fit is 
least satisfactory in those regions of the contour where 
the field intensities vary most rapidly) . 

Thus, the discontinuities in H are predictable and can be excluded from 

z 

the smooth profile of H whenever they occur. 

z 

Figs, 77 and 78 show that when the topographic profile has been 
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sampled 29 times. and the cylinder contour 30 times, an error of less than 
1% has been obtained in estimating while an error of less than 6% has 
been obtained in estimating M s . However, since the peak value of M g is 
about 1/3 the peak value of K^, this decreased accuracy in estimating M g 
is not surprising. 

Figs. 72 through 78 have demonstrated that the formulation of the 
problem of scattering from cylinders in a conductive half-space reduces 
to that for scattering from cylinders in a conductive whole-space if the 
conductivity of the air is set equal to the conductivity of the ground. 

A second check on the formulation of the problem and the program is to 
set the conductivity of the cylinder equal to the conductivity of the 
ground. In this case, the solution must reduce to the solution for plane 
wave reflection from a homogeneous half-space. 

Fig. 79 shows that the solution does re'duce to that -for reflection 
from a homogeneous half-space. In this example, we have assumed normal 
incidence of an E^-j^olarized plane wave and that: the height Z Q of the 

point of observation above the half-space is 1 m, the depth from the 
half-space to the top of the cylinder is 20 m, the cylinder radius is 
100 m, the incident field frequency is 1000 hz, the conductivity <p o of 

air is 0., the conductivity < 3“, of the half-space is 10 mhos/m and the 

-3 

conductivity of the cylinder is 10 mhos/rn. 

It is seen that an error of about 3 % has been obtained with 20 
sampled values on the half-space contour and 40 sampled values on the 
cylinder contour. This accuracy is not surprising since it was observed 
m section (5-3) that solution accuracy was poorest when the reflection 
coefficient was small. (In this case, the reflection coefficient is 
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zero.) However,.. it does indicate that when the reflection coefficient on 
the cylinder is small, the numerical results will have an error of less 
than 3% of the half-space value. 

Figs. 80 through 84 give the convergent numerical solution for scat- 
tering from a circular cylinder in a conductive half-space. We have 
assumed normal incidence of an E^-polarized plane wave and that: the 

height Z Q of the point of observation above the half-space is 1 m, the 
depth Z^ from the half-space to the cylinder is 20 m, the cylinder radius 

is 100 m, the incident field frequency is 1000 hz, the conductivity 

-3 

of air is 0., the conductivity of the half -space is 10 rnhos/nv, and 

"1 

the conductivity of the cylinder is 10 mhos/m. The topographic 
profile has been sampled 41 times between (-2100 2100) and the cylinder 
contour has been sampled 30 times. It is estimated that with these 
parameters an error of less than 2% of the peak value of H has been 
obtained . 

Fig, 80 shows that the peak value of real (H ) for a circular cylinder 

X 

in a, conductive half-space is about 10% greater than that obtained for 

scattering from a circular cylinder in a conductive whole-space when the 

\ ' 

field transmitted into the half-space is. taken as the field incident upon 

the cylinder. It should be remembered that this whole-space solution 

ignores coupling between the cylinder and the earth-air interface, and 

the boundary condition which must be met by the scattered field at the 

earth-air interface. The large response of real (H ) observed on the 
is due mainly x 

flanks of the anomaly^ to the magnetic field intensity reflected from a 
homogeneous half -space. 

Fig. 81 shows that the peak value of real (H ) is about 107 o less 
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than that obtained for spattering from a circular cylinder in a conductive 

whole -space when the field transmitted into the half-space is taken as 

\ 4 

the field incident upon the cylinder. As above, the whole-space solution 
ignores coupling between the cylinder and the earth-air interface, and 
the boundary condition which must be met by the scattered field at the 
earth-air interface. The reason that the peak value of (H^) is greater 
than the whole-space solution while the peak value of (H ) is less is 

Ct 

that the horizontal magnetic field intensity reflected from the half- 
space has not been added to the whole-space scattering solution. This 
has been done in Fig. 82 and it is seen that now the peak value of (H ) 
is about 20% less than that obtained for a cylinder in a conductive 
whole-space. 

Figs. 83 and 84 show that the effect of coupling between the cylinder 
and the earth-air interface on the induced surface current densities is 
to reduce the electrical surface current densities by about 10% and the 
magnetic surface current density by 20% in this example. Note from Fig. 

78, however, that the same accuracy cannot be attached to as to K^. 

It is easy to argue that the effect of coupling between the cylinder 
and the earth-air interface is to decrease the equivalent surface current 
densities. This is seen by examining the. polarization of the fields as 
they are modified by each boundary: 

1) When the incident field is transmitted into the conductive 
ground, the polarization of the wave remains unchanged 
at (+E , +H ). 

y f X 

The field scattered by the more highly conducting cylinder 
will possess an E -polarization 180° out of phase with the 


2) 
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incident electr-ic polarization. Directly above the top 
of the cylinder, the polarization of the scattered field 
will be (-E , +H ) . 

3) The field reflected by a dielectric medium back into a 

conductive medium will have the same electric polarization 
as the original wave, but the magnetic field intensity will 
undergo a 180° phase shift. As a result, the polarization 

i 

of the field reflected back into the ground directly above 
the cylinder will be (-E , -H ) . 

y x 

The "effective" incident field that is scattered by the cylinder is the 
sum of fields 1) and 3). Consequently, it is apparent that coupling 
between the cylinder and the earth-air interface will decrease the equi- 
valent surface current densities. 

The surprisingly small coupling effect can be reconciled by noting 
that the half-space contour in the immediate vicinity of the cylinder is 
in the static region ( of the fields scattered by the cylinder. 

In this region, it appears that the reflection at the half-space interface 
is small. By analogy, we can examine the response of a sphere in a con- 
ductive half-space to a static uniform electric field. The solution to 
the sphere problem has been investigated by Grant and West (1965, p. 425), 
and they note that interaction between the sphere and the earth-air inter- 
face is less than 10% if (Z^ + R) > 1.3 R, where is the depth from the 
earth-air interface to the top of a sphere of radius R. In Figs. 84 and 
85, <Z + R) = 1.2 R and it is seen that the half-space solution differs 

from the whole-space solution by 10% for K and about 20% for M . 

y s 

To illustrate the convergence of the method for VLF scattering 
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problems, the field scattered by a circular cylinder has been examined in 
Fig. 85 assuming ,that the incident field frequency is 30,000 hz. Further, 
we have assumed normal incidence of an E^-polarized plane wave and that: 
the height Z q of the point of observation above the half-space is 1 m, 
the depth from the half-space to the top of the cylinder is 19 m, the 
cylinder radius is 100 m, the incident field frequency is 1000 hz, the 
conductivity of air is 10 mhos/m, the conductivity Tj of the half- 
space is 10 mhos/m, and the conductivity of the cylinder is 10"' L 
mhos/m. 

Since the conductivity of the air and ground are taken to be the 
same, the numerical solution must reduce to the analytical solution for 
scattering from a circular cylinder in a conductive whole-space. It is 
seen that an error of less than 2 % of the peak value of real (H ) has 

X 

been achieved with 39 sampled values of the topographic contour between 
(-285, +285) and 30 sampled values of the cylinder contour. 

The field scattered when the vertical slab of Fig. 24 is placed in 
a conductive half-space has been plotted in Figs. 86 through 90. In 
addition, the vertical field scattered by this slab in a conductive 

whole-space (Fig. 62) has been plotted in Fig. 87 for comparison of the 

✓ 

half-space and the whole-space solutions. We have assumed normal inci- 
dence of an E^-polarized plane wave in this example and that: the height 

Z q of the point of observation above the half-space is 1 m, the depth 
from the half-space to the top of the cylinder is 20 m, the incident 

field frequency is 1000 hz, the conductivity^ of air is 0., the con- 

-3 

ductivity of the half-space is 10 mhos/m, and the conductivity 
of the cylinder is 10*"^ mhos/m. 
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It is evident that the main effect- of the_half-space is to reduce 
the amplitude and phase of the scattered field, as was the case for the 
circular cylinder of Figs. 80. In particular, it is seen that the in- 
flextions of imaginary (H) which were observed- for scattering from 
finitely conducting slabs in a conductive whole-space are still present. 

It was noted in section (5-4) that these inflextions m imaginary 
(H) were caused by each corner of the slab' responding to the incident 
field. This is seen most clearly in Fig. 89 where the equivalent elec- 
tric surface current density on the right half of the slab has been 
plotted. It is evident that the upper right hand corner of the slab, even 
though it has been smoothed, brings about a negative peak in imaginary 
(Ky) at a distance of 25 m' around the contour from the top of the slab. 

The upper left hand corner of the slab will bring about a similar res- 
ponse. 

Thus, the observed inflexions in imaginary (H) can be explained by 
treating the peak m imaginary (K^) at each corner as a line source. 
Whenever the radius of curvature is small, the current density will 
remain finite but exhibit an increase in magnitude. It should be noted, 
however, that the work of Mei and Van Bladel (1963b) in studying scat- 
tering from perfectly conducting rectangular cylinders indicates that 
this behaviour does not occur for H^-polarized incident fields (or at 
least it will not be as marked). 

The field scattered by the topographic profile of Fig. 47 has been, 
considered in Figs. 91 through 95. In this example, we, have assumed 
normal incidence of an E^-polarized plane wave and that: the height 

Z q of the point of observation above the half-space is 150 m, the 
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incident field frequency is 1000 hz , the conductivity of air is 0., 

-3 

and the conductivity try of the half-space is 10 mhos/m. 

In Figs. 91 through 93, imaginary (H x ) has been examined to establish 
solution convergence. Since imaginary (H ) is less than 10% of real (H ), 

X 1 X 

a solution convergence of 1% in Figs. 91 through 93 establishes conver- 
gence accurately for other components also. It is evident from Figs. 91 
and 92 that an error of about 3% (at the peak value of imaginary (H ) 
can be obtained if the contour is sampled 33 times within the contour 
section (-1100, +1100). 

The decrease in accuracy obtained when the contour is sampled 43 
times is a result of integration inaccuracies as was discussed for 
Figs. 72 through 74. This is seen even more clearly in Fig. 92 where 
the solution has been obtained for an Re. of .3 and .55 when the contour 
has been sampled 57 times between (-1100, +1100). Note that the largest 
departure from the- convergent solution occurs in both Fig. 91 and 92 
at x = -900 m, since the integration inaccuracies occur at the edges of 
the interval (-1100, +1100) where large interval widths have been taken. 

As in earlier examples, it can be demonstrated also that this is an 
integration accuracy problem by increasing the sampling density. Conse- 
quently, the contour has been sampled twice as frequently between 
300 4 I xl < 1100 as for N = 57, and an Re. of .3 was retained. It 
is evident from Fig. 92 that an accurate solution has been obtained in 
this manner. 

Fig. 93 indicates that, a contour section bounded by (-5000, +5000) 
is sufficient to describe accurately the unknown surface current density. 
On the basis of Figs. 91 through 93 it is seen that an error of about 3 % 
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(at the peak value of imaginary (H ), or about .2 % at the peak value of 

X 

real (H ) )' can be achieved if the contour is sampled 43 times within 
x 

the contour section bounded by (-5000, +5000). 

! 

In Figs. 93 and 94, the magnetic field intensity scattered by the 
hill has been plotted assuming a constant flight level of 150 m above 
the half-space and -a contour flight level of 150 m above the topographic 
profile. It is seen that as in Fig. 51, the peak electromagnetic response 
of H is reduced by a factor of two by contour flying. However, even by 
contour flying, there has been an approximate increase in the peak value 
of real (H x ) of about 6% as a result of topography. It should be noted, 
however, that this value is approximately one half the response predicted 
in Chapter 4 when the low surface impedance solution was applied to this 
problem. (It is interesting to note that in a similar problem. Ward 
(1967b, p. 271) deduced a topographic response of 10% on the basis of 
physical reasoning alone.) 

It is evident from Fig. 95 that by contour flying the peak response 

of H is reduced by about 10% and the position of the peak value (H ) 
z z 

shifts away from the cross over. The shift in peak H away from the 

z 

cross over is in accordance with the fact that the point of observation 
is further from the current sources than when the fields are observed 
on a constant flight level. 

To illustrate the application of the method to the most general 
problem of scattering from cylinders in a conductive half-space, the 
field scattered by a vertical slab within a hill has been plotted in 
Figs. 97 through 102. The geometry of the problem is shown in Fig. 96. 

In Figs. 97 and 98, we have assumed normal incidence of an E^-polarized 
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plane wave and that: the height Z Q of the point of observation above 

the half-space is 150 in, the height from the half-space to the top 

1 

of the cylinder is 50 m, the incident field frequency is 1000 hz, the 
conductivity <3^ of air is 0,, the conductivity cr, of the half-space 
is 10 mhos/m, and the conductivity of the cylinder is 10~^ mhos/m. 

The field scattered by the hill alone has been plotted for comparison, 
and it is evident that with these electrical parameters the fields scat- 
tered by the slab predominate over those scattered by the hill. It is 
not evident from a study of imaginary (H ), however, that the cylinder is 
a slab. The reason for this is that it could be argued that the inflexion 
of the peak value of imaginary (H^.) arises from the contribution of the 
fields scattered by the hill. Nonetheless, it would be possible to 
identify the scatterer as a slab. since the hill does not give rise to 

the observed inflexion in scattered H . 

z 

In Figs. 99 and 100, the conductivity of the slab has been increased 

from 10 mhos/m to 10. mhos/m. All other parameters of Figs. 97 and 98 

remain unchanged. It is seen that the amplitude of the scattered fields 

has increased and the inflexions in imaginary (H) are absent. It is not 

surprising, however, that these inflexions will disappear as the slab 

conductivity is increased since we observed in Figs. 25 through 28 that 

they are not present when the slab is assumed to be perfectly conducting. 

In Figs. 101 and 102, the conductivity of the half-space has been 
_ 3 / 

decreased from 10 mhos/m to 10~ mhos/ra, but all other parameters of 
Figs. 97 and 98 remain the same. It is evident again that the amplitude 
of the scattered fields has increased and the inflexions in imaginary 
(H) are absent. In addition, it is seen that the fields scattered by 
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the hill are only about 2% different from the half-space value. Note, 
however, that topography still is important since if modifies the field 
transmitted into the half -space from a plane wave to a wave propagating 
normally away from the hill profile. 

The increased amplitude observed when the ground conductivity has 
been decreased is expected from Tig. 36 since a A , where A is the 
largest dimension of the body, has increased by a factor of 3. The dis- 
appearance of the inflexions in imaginary (H), however, was not antici- 
pated. A study of imaginary (K ) around the slab contour indicates that 
there is a slight increase in the magnitude of imaginary (KA) at the 
corners, but not as marked as in Figs, 97 and 98 and not sufficient to 
be observed on a plane of observation 150 m above the half-space. 

Mei and Van Bladel (1963b) have shown that the current density is 
singular at the corners of a perfectly conducting rectangular slab when 
the incident field is an E^-polarized plane wave. Figs. 25 and 26 indi- 
cate that the current density on a perfectly conducting slab is well 
behaved if the corners are rounded as much as in Figs. 24. 

Figs. 99 through 102 show that for a large reflection coefficient, 
no inflexions in imaginary (H) are observed. However, Figs. 97 and 98 
indicate that inflexions are observed with a smaller reflection coefficient. 
Thus, we conclude that the inflexions observed in imaginary (H) over wide 
slabs are dependent upon the reflection coefficient and the radius of 
curvature at the slab corners. 



CHAPTER 8 


CONCLUSIONS 


8-1 Summary 


The purpose pf this dissertation has been to consider the theory of 
integral representations as applied to the solution of two-dimensional 
geophysical scattering problems. The examples given have demonstrated 
that for two-dimensional source problems, equations (2-56), (2-57), 

(2-66), and (2-67) are the most general integral representations of the 
fields m the exterior (source) region. Similarly, equations (2-72), 
(2-73), (2-74), and (2-75) are the most general integral representations 
of the fields in each homogeneous interior (source free) region. In 
addition, it was shown in section (6-3) that these two-dimensional inte- 
gral representations can be used to solve for scattering from cylinders 
assuming three-dimensional source configurations by expanding the primary 
current distribution and the field it radiates into a Fourier integral 
over a continuous mode distribution. 

Using these integral representations, it was shown that the solution 
to any two-dimensional scattering problem could be obtained by choosing 
an appropriate integral representation in each homogeneous region present 
and solving the resulting integral equation. It was found that for plane 

wave scattering problems, an error of less than 5% can be obtained without 

-2 


difficulty if <T^| > 10 mhos/m or /q~ 10 for frequencies up 

to 30,000 hz. Since these cylinders constitute a large proportion of the 
scatterers encountered in electromagnetic prospecting, this solution 
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accuracy is considered satisfactory. In those cases where < 10 mhos/m 

or /<T i < 10, the reflection coefficient of the cylinder is 

’dr • 

small and solution accuracy decreases. 

Solution convergence generally was established in several steps: 

1) The numerical problem was examined first for the case of 
scattering from a circular cylinder assuming that the 
cylinder was about the same size and had the same electrical 
parameters as the general problem. The conductivity of the 
air was set equal to the conductivity of the ground and the 
numerical solution was compared with the analytical solu- 
tion for scattering from circular cylinders in a conduc- 
tive whole-space. In this way, a general impression of 
solution accuracy and convergence was obtained. 

2) The width (-a, a) of the contour .section required to des- 
cribe the unknown surface current densities on the half- 
space was determined for the general problem. 

3) The integration accuracy on both the topographic and cylinder 
contours was decided. 

4) The distribution of sampling points necessary on each con- 
tour was found by increasing N. Although general rules 
cannot be given to suggest the sampling densities that 
•might be required, several observations can be made: 

a) The magnitude of the equivalent surface current 


densities will be large and vary rapidly on those 
parts of the contour having a small radius of 
curvature (see Pigs. 25, 26, 88 and 89). Conse- 
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quently, a small sampling interval is required in 
these regions to represent accurately the surface 
current densities. 

b) The equivalent current densities are equal to the 
tangential components of the total electromagnetic 
field quantities on the boundary. Thus, in those 
problems where attenuation is significant or the 
incident field strength falls off rapidly, the scat- 
tered field is determined primarily by that portion 
of the scatterer contour closest to the source and 
point of observation. Consequently, small sampling 
intervals are required in these regions only. (For 
example, see the fields scattered by a circular 
cylinder in the presence of an electric line source. 

Figs. 65 through 68.) 

5) Although not always investigated, an estimate of the size 
of numerical round-off with the above numerical parameters 
can be obtained by setting all conductivities equal to that 
of the half-space (see Fig., 60). 

It is important to note that in those examples where integration accuracy 
was not maintained (see F,igs. 13 through 18, 72, 73, and 88), the result 
was apparent as an oscillatory behaviour in scattered H when none was 
anticipated. 


It would be misleading to state that the program accompanying this 
thesis solves the problem of scattering from a vertical slab buried under 
a hill (Figs. 95-97) in 2 minutes on a CDC 6600, or equivalently, about 
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10 minutes on a CDC 6400, The reason for this is that the convergence 
tests listed above require a considerable amount of computer time before 
the desired equivalent surface current densities can be computed. However; 
it should be remembered that once an accurate estimate of K and M has been 
obtained^ the fields can be studied in detail using a small amount of 
computer time since it is only necessary to solve the integral equations 
once for each scattering problem. 

It would appear that solution accuracy can be improved and/or solution 
time reduced in several ways: 

1) The work of Zaki (1969) can be followed to account for contour 
curvature . 

2) The work of Green (1965) in studying finite difference problems 
indicates that it might be possible to extrapolate the results 
obtained at several sampling densities to a more accurate 
solution. 

8-2 ■ Conclusions 

As a result of applying integral representations to investigate 
geophysical scattering problems (assuming normal incidence of an E^- 
polarized plane wave in most cases); the following important observations 
and conclusions have been reached. 

1) The phase of H is dependent upon the position of the observer 
in space; even for perfectly conducting scatterers (see 
Figs. 27 and 28). 

2) The phase of H normally is not zero; even for perfectly con- 
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ducting scatterers because of the importance of the 
axial electric field intensity > (see Fig. 37). 

3) Coupling reduces the induced surface current densities 
approximately 10% to 207 o assuming a standard ground con- 

-3 

ductivity of 10 mhos/m, a cylinder conductivity of 
10 mhos/m, a frequency of 1000 hz, and a separation 
of 20 m between the top of the conductor and the half- 
space (see Figs. 83 and 84). 

4) Inflexions in imaginary (H) can occur in the field scat- 

/ 

tered by a single conductor if the radius of curvature 
on the upper portion of a finitely conducting scatterer 
is small at several sites (see Figs. 86 “through 90). 

These inflexions are dependent upon the radius of cur- 
vature and the -reflection coefficient on the slab contour 
(see Figs. 97 through 102). Note, however, that these 
results probably are applicable to E y -polarized incident 
fields only. 

5) A smooth hill 600 m wide and 100 m high with a maximum slope 

of 31° increases real (E^.) approximately 11% over that field 

reflected by a flat half-space assuming a standard ground 
-3 

conductivity of 10 mhos/m and a constant flight level 
of 150 m. The result of contour flying over this hill at 
150 m is to reduce the peak horizontal response of the 

hill by a factor of two and yield a constant field intensity 

\ 

over the central position of the hill. The peak response 
of H z is reduced about 10% and shifted away from the cross 
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over (see Figs. 94 and 95). When the conductivity of the 
ground is reduced to 10 rahos/m, the fields scattered by 
the hill are less than 27» different from the half-space 
values (see Figs . 101 and 102) . 

8-3 Extensions and Applications 

The above observations and conclusions were reached in demonstrating 
the validity and utility of the integral representations. It is clear 
that many important questions will be answered when a general numerical 
analysis of scattering from cylinders in a conductive half-space is under- 
taken. Furthermore, the work can be extended to study some of the most 
basic scattering problems encountered in geophysical exploration. Some 
of the general problems which now can be investigated include the 
following: 

1) The effects of overburden on the scattered field can be 
studied like the half-space solution by choosing a fourth 
integral representation in the homogeneous overburden layer. 
Outside (-a, a) we assume that the layer is flat and the 
equivalent surface current densities arise from the primary 
field incident upon a flat layer overlying a conductive 
half -space. Alternatively, a finite overburden layer can 
be considered by treating the overburden as having a 
catenary-like cross section. 

The fields scattered by cylinders in a conductive half- 
space assuming H^-polarized incident fields can be examined 


2) 
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by programming the Hy-polarization solutions given in this dis- 
sertation. 

3) Scattering from multiple conductors can be considered by choosing 
an integral representation in each additional scatterer present. 

4) The magnitude, phase, and geometry of the fields scattered by 

cylinders in a conductive half-space assuming finite sources can 
be studied by following the methods outlined m Chapter 6. When 
the finite dimensions of the dipole become important and/or 
coupling between the ground and source is important, the tecta-- 
niques of Chapter 6 may not be satisfactory. Instead, the 
problem should be considered as an antenna scattering problem 
and the current distribution on the antenna included as an 
unknown. Work in this direction for dipoles over a lossy earth 
has been reported in an abstract by Arens and Embry (1968) . 

5) An analysis of scattering from finitely conducting bodies 
of revolution can be undertaken by extending the work of 
Andreasen (1965 a) to include interior integral representations 
and applying the results to geophysical scattering problems. 

Note, however, that this solution requires that the incident 
field be expandable into a set of orthogonal TE and TM modes 
propagating along the symmetry axis of the body considered. 

6) the- time- response' of 'fields scattered from cylinders of arbi- 
trary cross section can be studied by transforming frequency 
domain results into the time domain. Although this may not be 
the most economical manner in which the problem can be solved, 
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it is straight-forward in principle. 

7) Since elastic waves obey the same Helmholtz equation as electro- 
magnetic waves, the analysis of this dissertation can be repeated 
to investigate seismic scattering problems. Acoustical scatter- 
ing problems have been investigated in some detail in the liter- 
ature (for example, see Mitzner (1967b), Copley (1968), and 
Schenck (1968)), but the results were not applied to seismic 
exploration probl ems . 

In each of the above cases, an analysis should be under- 
taken to determine where approximate solutions that have been 
developed are valid and more expedient. 'In those situations where 
approximate solutions are not available, it is important to deter- 
mine when approximations can be made to the numerical scattering 
solution. For example: 

a) The much simpler solution discussed in Chapter 4 can be 
used to describe scattering from the cylinder (but not 
the half-space) whenever the -cylinder can be treated 
accurately as having a low surface impedance. 

b) If coupling between the cylinder and the half-space can 
be ignored, the scattering problem can be considered m 
several parts and a solution can be obtained' more ac- 
curately and quickly than if the entire problem were sol- 
ved once. 

The application of the results of this thesis to plane wave scatter- 
ing by cylinders encountered in geophysical exploration leads to the 
following important conclusions: 
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1) In AFMAG surveys, the ratio of real (H) to imaginary (H) is a 

function of traverse position x and ground conductivity , as 
well as the cylinder conductivity In addition, the magni- 

tude of the scattered field is highly dependent upon the con- 
tribution of the electric field intensity E, so that it cannot 
be ignored. As a result, investigations such as that by Ward 
and Fraser (1966) in which the contribution of the electric 
field intensity and host rock conductivity have been ignored 

are likely to have a limited application to AFMAG interpretation 
for cylinders . 

2) Topography can give rise to a tilt angle of about 5° at an oper- 

-3 

ating frequency of 1000 hz-., ground conductivity of 10 mhos/m. 
and normal incidence of an “Ey-polarized plane wave according 
to the analysis presented in Chapter 7. However, before these 
results can be applied directly to AFMAG interpretation, an 
investigation must be. undertaken in which the incident plane 
wave can assume an arbitrary polarization and angle of incidence. 
The reason for this is that- in air, the incident AFMAG fields 
are propagating at grazing angles to the earth-air interface 
with a vertical electric and- horizontal magnetic polarization. 
This direction of propagation, will not affect the conclusions 
reached in 1-) since the field transmitted into- the half-space 
will be propagating normally away from the earth-air interface 
even for grazing angles of incidence. Thus, both Ey and Hy 
polarizations are possible. 
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3) In no case has a zero phase been observed. Even perfectly 
conducting scatters buried in a conductive half-space will 
produce an out~of -phase response for plane wave incident 
fields. This result is particularly significant since early 
work such as that by Ward and Fraser (1966) had indicated that 
AFMAG fields scattered by a perfectly conducting cylinder would 
yield no quadrature response. 

4) Conclusions 1) through 3) are equally important for V.L.F.- 
E.M. studies. The conductivity of the bedrock and the elec- 
tric field vector must be included in any analysis, and the 
incident field should possess an arbitrary polarization and 
angle of incidence. 

The above extensions and applications indicate that integral 
equation formulations represent a powerful technique for solving 
many of the complicated ’'electromagnetic scattering problems encounter- 
ed in geophysical exploration. By verifying approximate solutions 
which have been developed and continuing the analysis of this disserta- 
tion where none exist, electromagnetic interpretation in mining geo- 
physics can be based on sound theoretical investigations. 
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APPENDIX A 


DERIVATION OF THE TWO-DIMENSIONAL GREEN'S FUNCTICN 


The two dimensional Green's function in an infinite region is given 

! 

in many texts (see Noble, 1962, p, 239), but the proof is generally left 
to the reader. For this reason and for the sake of completeness, a general 
derivation of the two-dimensional Green's function is given. It follows 
a development similar to Fuller (1968) in obtaining the- three-dimensional 
Green's function. 

We wish to obtain a ^particular solution to . the inhomogeneOus scalar 
Helmholtz equation, " • 




(A-l) 


Assume that both G(x,y,z; x',y',z') and $ (x - x') £ (y - y')cT (z - z') 
have a Fourier Transform given by 
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-Substituting equations (A- 2) into (A-l) and carrying out the 
operations, we find that a sufficient condition for a solution is 
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Transforming to polar coordinates such that the position vector A 

A A * 

is xi + yj + 2 k, then 

i/k - m' ) (.2-2') - j? ( . ( a - A. ) = 1?, 1 k-J.' I tos© ^ 

4 tj 1 ^ ■= H, X , (A-5) 

&K* - Jk* Slm & J.& d.^) . 


Equation (A- 4) becomes 

,^o ^ f zr ivA-vS'lcos© ^ 

<S Ca, a') - L_ 

Liuf 



o Jo 


iv-n 


£wv© Jlk { Afedity 


_J 

(Att) 1 



iUa-,*' / A 2. 

.^'P. £'n\ 


(C-t 1 ) 



- 1 

nr* lx -*' I 




- -- L 

r r w> n it, ( A-^'l 

1 f A* it, * 

r° C* U'Vl A 

«C ^ oU 

<^7 T^Ia-A 1 

L J> J 

U x -P) J 


- c 

■4 Tf 1 I a- a'J 



(A- 6) 


The integrand of equation (A-6) has poles at ± k where k is a complex 
number given by . In the complex plane, = Ik-J =|fc,l^os ©j + 1 sm©,). 

For 0^ S, - W } sin& ( is positive, so that 


C^Ja-a'I } I a - a'I IM c<>s ©, - iX-x'l Ik, I sew©, 
JL ~ *- 


(A- 7) 
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Thus, the integrand of equation (A-6) approaches zero as k^— «. oo , 
and we can add the dashed contour shown without changing the value of 



Thus, by Cauchy's Residue Theorem, equation (A-6) becomes 

i kU-;;'! 

a Jt. (A-8) 

-t TT Ia-a' I 

Equation (A-8) gives the three-dimensional Green's function m an 
infinite region. The .two-dimensional Green's function is obtained by 
integrating out the axial dependence of equation (A-8)'. Thus, in two 
d imens ions 


= * 


,<*> it /t*-*’)*' + t 1 i (2 ~ Z ' ) 


At 


(At 9) 


/ (/<- a') Z + i 4- (B - Z ') 4 


where 


t - 




Ward (1967, p. 131) shows that 


1 


00 ikfc /— — \ 

JL_ J* = _i_K 0 C-ik /* +J( ) 


%T\ 


4T J-c« ^ 

where R -JtT + y + z“ and is a modified Bessel function 


(A- 10) 


2 . ~2~ 2 -1 


of the second kind of order zero.. 



i/a 


We have the relationship from Watson (1966, p. 78) that 



± 

X 


uzrr 

ti JL 



^ - n < i r/z 


Thus, 



l (kj {*-*'>* +1*-*^) * ± 

7 < 


(A- 11) 


(A- 12) 
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APPENDIX B 

ALTERNATE DERIVATION OF THE FIELD SCATTERED 
FROM PERFECTLY CONDUCTING CYLINDERS 


We wish, to solve for the field scattered from perfectly conducting 
cylinders in a conductive whole-space, assuming an incident electric field 
polarization which is parallel to the axial direction of the cylinders. 

To support this scattered field, we postulate the existence of electric 
currents which are induced such that the boundary condions on E and H 
are satis'fied. 

We have assumed that E 1 possesses an axial component only. Further- 
more, -we will assume that the incident field is constant in this direction. 
As a result, the induced currents will possess an axial component only 
and, from this, the scattered electric field must be axial also. 

This scattered electric field must satisfy an inhomogeneous Helmholtz 
equation given by 


{ It" ) Ej 5 ‘ C (a ) a - Kj Ca) 

To solve equation (B-l), assume that both E SC (x,y,z) 
Fourier Transform given by 


(B-l) 


and K^x,y, z) have 


( a ( -J 


OQ 







te 




(B-2) 


« 

When equations (B-2) are introduced into equations (B-l) and the 
operations carried out, we find that a sufficient condition for a solution 
is 
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~ fe e ). ^ 


(B-3) 


where 


5 W) = — . 


-t' 

Then, by the convolution theorem, we can write the solution in (x,y, z) 
space as 


£j C ^M> 2 i = Ky^*>3> z ) 


(B-4) 


where G(x,y,z) xs the three-dimensional Green's function and 
is the inverse Fourier Transform of gO 5 ^ ky> k z ) . 


Thus, 


= — - 


l(^V 




ut) j jjj ~ J (B ‘ 5) 


We have evaluated equation (B-5) in Append xx A so that we may write 

(?£<■, j,?) ~ -£- . 

4 TTVi (B-6) 

When equation (B-6) is introduced into equation (B-4), we find that 
the scattered electric field intensity is given by 







2 .) - 



U-A 'I 


-/ 








(B-7) 


where 
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However, since l^x'^y'^z 1 ) is constant along the, axial direction, 
equation’ (B-7) reduces to 


ty 2) = iM-U> 


re r° iL / iy-yJV cs'^o 2 n 

tyuo ji - Jty 1 d*.' A?' 

-nr JJ J Tcz^PT^ 

—cO -00 J J 


(B-8) 


This integral has been evaluated also in Appendix A, from which 
we may write 

£y Cfcjg) = -><.0 jj K^C^e') Cie + C£-z'-> . (B-9) 

Since the cylinders are assumed to be perfectly conducting, K(x,'z') 

3 

is a surface current density. Thus, equation (B-9) reduces to a contour 
integral representation for the scattered electric field intensity. 


by JlU 

A 


tyf') K ( 'l) ds > 


- (B-10) 


where 


I j-p'l - 
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APPENDIX C‘ 


SCATTERING FROM A FINITELY CONDUCTING CIRCULAR CYLINDER 


The solution for plane wave scattering from a finitely conducting 
circular cylinder in a conductive whole-space is very important for 
testing the validity of the integral equation scattering programs which 
we have developed. A solution which can be compared with that given 
by Wait (1959) is given below. 

If we assume that the incident plane wave is E -polarized, then we 

y 

can write, referring to Fig. 103 




(C-l) 


Using the cylindrical wave transformation given by Harrington (1961, 
p. 231), equation (C-l) can be rewritten as 

~ ;- n J /L ^o cnq/ s f J 

(C-2) 

where 


'> rftpoo > o I l ‘ 

K - f 

v 2- , (\ > l • 


The total electric field intensity at any point in space is represented 

i 

by the sum of an incident electric field E and a scattered electric field 

y 

sc J 


= + V * 


(C-3) 
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To represent outward-travelling waves external to the cylinder , the scat- 
tered field must be of -the form 




(C-4) 


Inside the cylinder, we expect standing waves which are finite at the 
origin. Thus, we seek solutions of the form 

ef - § • (c 

J sh.~0 1 

At the cylinder boundary, the condition that .tangential E is con- 
tinuous must be met. Since we equate' each coefficient of cos#ify , 
this implies that 


( k, R) + '* (G-6) 

To obtain a second equation between a n ' and b n , we will enforce 
continuity of tangential H, that is . We have from Maxwell's first 
equation that 



Jl to 




Thus, from equations (C-2), 


( _[ kg ^ 0 

' j 5 "bty 

(C-4) and (C-5), we have that 


= ~L JL 

A* 




CO 




/h>0 





JA,U> 


(C-7) 


(C-8a) 


(C-8b) 
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In equations (C-8), differentiation is with respect to . Continuity 


of tangential H implies that 


^[-L' (fe f R) + <vO k « R) I = z x CkR) , 


(C-9) 


where we have equated each coefficient of cos and now differentiation 
is with respect to the argument (kR) of equation (C-9) . 

Solving equations (C-6) and (C-9) for we find that 


= * 


4. U, R ) - c n 

-<a 


(C-10) 


where C, * , 

2, -4 (M> • 

m_ 4, Ck|i^) 
i.R 




dJM- , 

4 R 

ifCM'- £(W 

■ n 


If we set 


££ * z, Ho «, 


- (C-ll) 
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then the scattered electric field intensity is given from equation (C-4) 
as 



,?)* Ho 



o 





(C-12) 


where a is given by equation (C'-IO). 
n 

From equations (C-7) and (012), the transverse magnetic field inten- 
sity is given as 


and 


i,A 


ikk 


Z- L 

l 


~ K a /* H* C^) s “' /n $ 


(C-13) 


(.«)- id./' 11 2- 


H jC ft $ 


<3-„ 


SK-O 



where a is given by equation (C-10). 
n 

The equivalent electric and magnetic current densities on the 
surface of the cylinder., assuming Ey-polarization, are given from equations 

(2-38a) and -(2- 3 9b) as 

• r 

Ky * U s * If Q > 

Ms ~ ~ • 

Thus, from equations (C-8a) and (08b), the equivalent electric current 
density is 

-J ' ■ / h<0 fr> 


(C-15) 

(C-16) 


and from equations- (C-2) and (012), the equivalent magnetic current density 
is 
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. 1 r <20 . 

K s (^) - ^ l' + 4 


<*»«> £>‘ 




Cos /H 


4“, 


(C-18) 


where a is given by equation (C-10). 
n 

The derivatives in equation (C-17) are with respect to the argument 
(kjR), and are given in equation (C-10). 

A useful approximation and check on our solution is to study the 
asymptotic results as •» oo . In this case 


IU - 



S <=© . 


Since 

-vJ^Ckzf?) - f A Co * ( k ( ft - If ~ 'H ir j } 

A— /irk t R i z 


(Watson, 1966, p. 195), 


then 

£i~ 

4,—'“ 



and 

Ckl £) 

^ (M 


Cos ( k, ft - V} , ^ Ve ) 

cos (fe, ft - V? - £'*»«■») Vi ) 


(C-19) 


(C-20) 


(C-21) 


(C-22) 


Thus, 

Cos ( fe, R - - * Vi) o 

^ ^4,. Cos (li, R - V$ -N+i) Vi ) 




(C-23) 
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and hence the reflection coefficient reduces to 


^ ^ ( k > p) 


K Ck,e) 

Except for a different time dependency, this is the same reflection 
coefficient as obtained by Harrington (1961, p.233) for scattering from 
perfectly conducting cylinders. 

Introducing equation (C-24) into equations (C-12), (C-13), and 
(C-14), the scattered field quantities become 


(C-24) 




t ...» 


f'H^O 


(C-25) 


• -2ih 2- rffCk/ 1 ) i 


and 


w„% 4 ) . -i ^ ,k %c\ 4M k - CA^-4 

(Af i 


The equivalent electric and magnetic current densities reduce to 


(C-26) 


(C-27) 


;W 


and 


KqU)-* 2 (^))(OS / H^’ > (C-28) 

' ^ W) 7 


Ms (#) = O . 


(C-29) 


If we introduce the Wronskian relationship (modified from Watson, 


1966, p. 77) 



(C-30) 


4» (fc *) C, W- -L 0, » * IM . 

Tt,k 


then 


4,'tw IU*0 - J«mC(k,V - -!jkx)n~%x) 


= t_AL 

rk,R 


(C-3 1) 


Thus, equation (C-28) reduces to 


Kj ($> - *tf. 

rfc,R 


l^k 

JL 


CO 

z 

y 9 * *0 



£,©S<«(-P 

iFTPJ 


(C-32) 
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APPENDIX D 

SOLUTION FOR POINTS OF OBSERVATION CLOSE TO THE SCATTERER CONTOUR 


When investigating the scattered fields close to the scatterer 
contour (lkf||r-^'i < < \ ), we found that a parabolic approximation to the 
integrand of each integral was inaccurate across those intervals near 
the point of observation. This was due to the pseudo-singular behaviour 
of the Hankel function, and we stated that this problem could be overcome 
by making a small argument approximation to the Hankel function in this 
region of the contour integration and integrating the resulting expression 
analytically. 

To integrate analytically, it is convenient to translate the 
(x - x 1 , z - z ') coordinate system to one which has its axis parallel 

' A 

to the contour interval and its £ axis normal to the interval and through 
the point of observation. Thus, it is evident from Fig. 104 that 
this yields 


' / 

5 



- (js-f'V 


A/ 

s 

. £ 


(D-l) 


Since 

& ST CoS £ oC -f -f Svn j£ — -Sivi ot A. + Cos at. Z 

8 4 + Cos Coif s — Coi ca -x — 2 


> 


(D-2) 


we find that 




CbS. ai 


(D-3a) 


< b — — j {/K- /v ! ) L° s <k * oij 


(D-3b) 
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Similarly 


(✓*--*') zr s's^Kot - § doS ^ 

(g-e') = — 4 S / dbSoi! +5 iixa;} 


(D-4a) 

(D-4b) 


D-l Perfectly Conducting Scatterers: E y -Polarization 


The problem of the pseudo-singular behavior of the Hankel function 
first arises when estimating the surface current density on a perfectly 
conducting scatterer assuming E^-polarization. Thus, we must evaluate 

.b 


ejlf) *£k( K Ijlfiuflklf-ft)*'' , 


4 

(D-5) 

for « 1 . The surface current density Ky(^') over this 

interval is given from equations (3-19) and (3-21) as 


kj Cs') , | ■ U-A,Kdc : i . U'- «:;)• ] K 

L H % -I J 


b- + 


LijCl+AjJ tj Cl+Ap. 


k« 

Jj+i 


k. , 

Aj(»+Xj) ( i+y J 


(D-6) 


where 


A j ’ W j * *S« > 

A i * a t v , 


/V 


w j +vv j-l 


- (Wj +v# ifl y« , (j+w l . I ) , 

<£-, - ) £Vh oi- — (2-2- ) 6»Soi. 

J J J a J 

When equation (D-6) is introduced into equation (D-5), we find the 
following solution to our problem: 
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1 (o'O S yxw f Hkdo - (t-Ap (hKq! * gj ttfcQfr) - (. H«o^- - -3d. V c~Hk*o )l k fj 

T U A i % J 


♦ A 


i±t 


kAj 


iifce-l - Cj tf£t>0 ^ HK6?. - HKbJ + tj HKOO 


4 j 


\ Vr 

A- (Hfco‘-c’s - HK62- 2c,/f/roi 

l 

HAj 

A- T; J . 

^ J J 

j (D-7) 

where 

HK»« = 3 

J CL 

(D-8a) 


HKot * j 4 t* }£ > tt/}> l 4& ,l )ic' , 

(D-8b) 


**» . f b ** c 

(D-8c) 


X 



Equations (D-8) can be evaluated analytically after H 0 (k/5 ) 

has been approximated by equation, (3-26). However, the following inte- 
grals will be required to carry out this integration: 


fjU d ( S ) + c , (d-9.) 

J 2 1 4 


fortes’. ifV'Xdt U'M'Au*"-- S l X(s'*u ! >) 

2 2 


i |VjU U'*sY\ I /_ s V.ftrV j 

2 3j[ 1 ^ s 5^ 


+ C > (D-9b) 


C (D-9c) 
+ ) 


/turf -i/ihsVfStiu'is^lLc/D-sd) 

y 11 2 2^ -J 
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(2ms' 1 ) JU'* ±Js*JL]ik(g / \S £$'* + l 4 s‘ _ S^V)L C , (D-9e) 

2 si * ^ 3 S J 


s'Sirf (i\s x )%' * i [s'jLYi UV) V ll(/_ -Sl(s't^))] + C > (D-9f) 

2 6 [ 2 2 ^2 2 J 


where we have used the facts that 


J j>(5 ') U U (s V' 2 )^' . AXkJjLs')JU f Jju'Un 

j A j(s)ls' = s'JkjCsO -y s'l'Ul l*' > 

3 ^ 

ha') AjtsOds' = fG')JL j(s') - J 7 rtspyCsO A ' , 

3 


(D-IO) 


(D-ll) 


(D-12) 


with 


FCs‘) 


and g(s') and f(s') are assumed to be rational functions. 

Introducirg equation (3-26) into equation (D-8a) and carrying out 
the integration, we find that 

H Koo ~ / |_ hV)( U*) - A 1 (k 3 -* 3 ) f A. | +Jl U*' 4 ^ | 

\ 4 ll 2F L 3 J 


' i c 1 ) ^ 


at f b(i_ £S‘ _ »(i- Ef- 


+ £ . (i- ££*)(!.— | 



Similarly, introducing equation (3-26) into equation (D-8b), we find that 



Finally, introducing equation (3-26) into equation (D-8c), we find that 

HKOI ^ I /, . i^y b‘-^) - It 1 lb S - + i f (iwu e 

3 1 A A 2irl ^ s 

+ 1L { i> 3 (l- £**- (tV') Vl - a*/ 

3ir 1 A i i 4 zo > 2 


♦- (i - il f)( .-U1)))|. (d is) 

Having estimated the current density on a perfectly conducting 
scatterer, we now wish to calculate the magnetic field intensity close 

to the scatterer. Thus, 

^ * 

N Co) - it- f i 

4 (D-16) 

must be evaluated for ) kll^t-p 7 ! J 

via equations (D-4), then 


. Translating our coordinate system 
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and 

Col A — /K~ ■/*■' - s $ 1/vt c<- — O CoS o< * 

Ip.jT'l (SV‘) W (D-17b) 

Thus, after approximating the current density by equation (D-6), the 
following form of equation (D-16)is obtained: 


_ - 1 k (ft Cc.5c< + ^ +- (ftSb*oA — SBcfeSoC^^j ^ 

^ J 


where 



(D-18) 


(D-19a) 



(D-19b) 

(D-20a) 

(D-20b) 



195 


4 ■'* »? (k/sW r )A\ 


'A ' 


WTPW 


(D-20c) 


Introducing equation (3-27) into equation (D-20a) and carrying out 
the integration, we find that 

k f (i_ - i - 1 - j; fi_ - t, j i) 

f? /I •* A^« I ^ i ^ 


(D-20d) 


r i 




f 


1A 


J_ h , S£(b 3 -«?)| 

1 ik % J 


3r 


$ 2A ' 2 V X x 




(D-21) 


Equation (D-20b) can be integrated analytically for all values of the 
argument, and is given by 

A 


j » . - 


i 


r jL f «? ( ust ). 


Equations (D-20c) and (D-20d) can be rewritten in terms of equations 
(D-8a) and (D-8b) through integration by parts. Thus, 


r* . . _L f l A , 

i L iA -l 


- J 


a Cik/^P) - i uftt/Fs) J i’a/A? 


fl (D-23) 



196 


and 



The remaining integrals m equations (D-23) and (D-24) have been estimated 
previously by equations (D-13) and (D-14). 


D-2 Perfectly Conducting Scatterers: H^-Polarization 


In estimating the surface current density on a perfectly conducting 
scatterer assuming H^-polarization, the expression 


* b 

^ if") = - K,(f) ~ d( P-25) 

3. A J/t 


mus 


2 A 'q 

t be evaluated for } kl ) jt"~ f>' I I » Since 

ot) - - S 


A 


•f S 


(D-26) 


and K s (p ) is given by equation (D-6), equation (D-25) can be rewritten 


as 


^ c*r> - - + iU b 

^ A 

(D-27) 

where B is given by (D-19b), 

The magnetic field intensity scattered by perfectly conducting cylinders 
assuming H^-polarization is given by an expression similar to equation (D-25) 
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W*Lj). i_J f Mf'J -*) (D-28) 

*7 - c 

As a result, it is evident that the integration across the jth interval 
for points of observation close to- the contour can be written as 

H,“ c ? ) • - i £4 B » 

J (D-29) 


D-3 Scatters with Low Surface Impedance: E -Polarization 

y 


The integrals that are necessary to evaluate the electric -field 
intensity for points of observation near the surface of a highly con- 
ducting scatterer assuming E^-polarization have been estimated also. 
The expression 



(D-30) 

must be evaluated for Ikllp-A 7 1 I . -The first integral has been 
estimated by equation (D-7) and the second integral is just Z ( p') 
times the integral of equation (D-28) across the jth interval. 

However, a rather complex expression results when we wish to estimate 
the scattered magnetic field intensity close to the contour of highly 
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conducting scatterers . The expression 



must be evaluated for tkll p-p'l^ I . Since the first integral has been 
estimated by equation (D-18), we will consider only the second integral 
of equation (D-31). Once the trigonometric functions are expressed in 
the (s, 1 &) coordinate system as 


CoS ^ ~ 

Co-i 


_ Cs Co s<=t + 

~ S CoS <X. 

CS%s'‘> Vt - 


Co ■* jj Ces oC 4- 5wi 


■SUWot — 


- 5 


C&N-s'*)'* 


77 > 


Son(^-^)s It'S ei - CottjS 


tl?S oL — Cert R oi — ~ S 




(D-17a) 


(D-17b) 


(D-26) 


(D-32a) 


si^tafr— 0 = vcoSjli<^C/a-««)= t (D-32b) 

a%s-) 

£oS<.2 ( &~oO r Co-S^ <^a(p-=<) - Svn ^ a (£*-*' )ce>S ~ , (h-32c) 

( s' J ) 
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the second integral of equation (D-31) can be written in the form 


X. 




k TS&Fi 


A. 

A. 


£>/ s^sCviot - ^ CoSot - Is'Ssw^j ^ (kffi 





(D-33) 


Thus, introducing equation (D-6) into equation (D-33), our solution is 


*k ^(j3 j) |TS( Cc*s<x * $ ^ •iUs.at. +iSGcosa)|^ 


+ j^S ( C%<.*ct - SD costi) +]^ € s.ia» ot)| 2 ^ (D-34) 


where 


e - [ It - ihhi (4 '<j j 7 (5 

A' t J 

J J 

J 7 ~ c i J * + -r 2c i V £i 


Jj 


* Jjiit 


*+ A: 


J 


- J g -2^j J 7 +glj , 

l4Aj L 4 - 1 

* - [% - (/-A;) 




+ F ^Sjg + j? ~ ZCj 3 ^ gj Ig ] 
^ L a 4 ^ J 




1 4 A L 


(D-35a) 


(D-35b) 



zoo 



(D-35c) 


(D-35d) 


(D-35e) 

(D-36a) 

<D-36b) 

(D-36c) 

■ (D-36d) 

(D-37a) 

(D-37b) 

(D-37c) 

(D-37d) 
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h Tan^) **' ■ 


(D-37e) 


Before equations (D-36) and (D-37) can be evaluated, the following 
integrals will be required: 


' JL X -T L^) %' , i(jLVk\^ l i_ / J ,_L ( JgwCs%^) i^ 4. 

^ z S ' 2 J (D-38) 

' * a 4 3 jusy* )(jtyk . 

£\s ' 1 a i 2 4 J (D- 39 ) 


where equation (D-39) has been obtained through integration by parts. 

We can show that the remaining integral in equation (D-38) does not 

possess a closed form solution by transforming variables according to 
i i- , 

s' = £x. Thus, 


\> ws k/5 y* 

P A 5. f = jill (aaT^x I + _L 




J 50+A 1 ^ 

d/T 


a /<T 


a 


t-s « l 

«}S 


(D-40) 


The last term of equation (D-40) is given by Grifoer and Hofreiter 
(1961, p. 112) as 
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where 


i a (*) . - ( X (m) <6s 
X * 



ft ~ 3r 




j^or j ^ | ^ 1 



(D-42) 


Since the evaluatioh of (D-41) would invdlve the numerical calculation 
of four integrals , the quickest numerical solution would be to evaluate 
equation (D-38) as it is written. However, by rewriting the integral 
of (D-38) as 


■/V3o 

Ms'S 

-As. 



f h. 

hy 

'30 





1 O.S)/ S , 

? Vo > J S/3o S -Ss‘ Ik 3oS J 


where we have assumed that a<S < b, 

the numerical integration is performed over the smallest width necessary. 

Introducing equation (3-26) into equation (D-36a) and carrying 

I ^ ' 

out the integration, we find that 


i s ^ Wft) - 

S S % t 4 T 

+ $ 1 (4»~' ^ oj) + 1 I Ujt) 

1 $ 2 & 5 2- 


_£f 

^ L 1 


a 


a $ £ ! 


l (D-43) 


where 


IU.M 




(D-44) 
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Similarly, introducing equation (3-26) into equation (D~36b), we find that 


~ J_ JL *>*+* * - |-*0 + l 

1 8 If y f 


JDtV.SL kjj 

1- & 




+ i 

4 


- e. l A.]A U l .5*)''* - i kVi! ).] 

1 2- «>+S Z J 


(D-45) 


By rewriting equations (D-36c) and (D-36d) in terms of proper 
fractions, we obtain integrals that have been estimated already. Thus, 


I, a 


/V (k]z\s tl )£ 
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where the first integral has been estimated by equation (D-13) and the 
second integral has been estimated by equation (D-43) . 

Similarly, equation (D-36d) becomes 
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where the integrals have been estimated by equations '(D-14) and (D-45) 
respectively. 

Introducing equation (3-27) into equation (D-37a) and carrying out 
the integration, we find that 
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where I(a, b) is given by equation (D-44) . 

Equation (D-37b) can be rewritten in terms of integrals already 
estimated if it is integrated by parts. Thus 
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where the remaining integral has been estimated by equation (D-45) . 

The remaining three integrals, equations (0-3 7c), (D-37d) and 
(D-37e) can be written in terms of estimated integrals by expressing them 
as proper fractions. Thus, equation (D-37c) becomes 
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where the integrals have been estimated by equations (D-21) and (D-48) 
respectively. 

Similarly, equation (D-37d) becomes 
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where the integrals have been estimated by equations (D-22) and (D-49) 
respectively. 

Finally, equation (D-‘37e) becomes 
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where 'the integrals have 'been estimated by equations (D-23) arid (D-50) 

\ « 

respectively. 
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APPENDIX E 


EVALUATION OF AN INTEGRAL OF THE SONINE-GEGENBAUER TYPE 


To obtain the desired integral, we will evaluate 




first. (E-l) can be rewritten in terms of the modified Bessel function 
since we have from Watson (1966, P.-78) that 
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Consequently, (E-l) becomes 
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However, since (E-4) can be considered to be a Hankel transform, we find 
from formula (35) of Erdelyi et al (1954, p. 72, Vol. 2) that 
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provided b> 0, Ira(k)> 0, Re(<5)>0, and Re^tt )>•-!, 
On rewriting (E-5) in terms of Hankel functions, we find that 

/ 7 °° . , jx , f— t-r—, 
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provided b> 0, Im(k)> 0, Re(S )> 0, and Re(a)>-1. 

Although it is not obvious, it would seem that equation (E-6) is 
valid in the limit as Im(k) approaches zero and Re(k) is greater than 
zero. This is easy to demonstrate for the special case when all the para- 
meters in (E-6) are real. 

LO 

Watson (1966, p. 179) gives an integral representation for H^(z) as 
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Setting m equal to zero in (E-7), then (E-l) can be rewritten as 
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On interchanging the order of integration and setting 6t= __*• . 

(E-8) becomes 
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(E-9). 


The inner integral of (E-9) can be evaluated from equation (4) of Watson 
(1966, p. 394), so that (E-9) becomes 
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provided Re (^ )> -1, Re(k)>0. 


Setting x = / - b^~ 
5 a 


and assuming for convenience that k>b, then 


(E-10) reduces to an integral representation of the form (E-7). Thus, we 


find that 
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provided k>bi>~0, S> 0 and Re(^)>-1. 

By assuming that k>b and that all parameters are real, we have avoided 
the problem of considering the validity of deformed contours for the 

integral representation of H (z) which otherwise would have arisen. 

0 

Consequently, it would seem that (E-6) is valid at least in the 
special’ case Re(k)>b>0, Im(k) = 0, and Re(S)>0. 

If we take the~limit of equation (E-6) as b-*-0, we obtain the formula 
of interest to our scattering problem. Noting that 
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then equation (E-6) reduces to 
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provided Re(k)>0, Im(k)> 0, Re(£>):>0, and Re(^)> -1. 

It may seem difficult to justify using (E-12) for all s' since 
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at some point during the integration the product bs ' will violate this 
small argument assumption. However , it can be shown by using a stationary 
phase argument (see Papoulis, 1962, p,. 139) that if b is small enough so 
that (E-12) is accurate up to that value of s 1 where the large argument 
approximation 
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is valid, then that part of the integral for which (E-12) is inappropriate 
contributes a negligible amount to the integral. 

It should be noted that (E-13) could have been established by rewriting 
(E-13) m terms of the modified Bessel function as was done -for equation 
(E-l).- Then, considering (-ik) to be a complex constant and setting 
x = (s ' 2 + £^)^^, the resulting integral can be treated as a K-transform 
and evaluated using formula (13) of Erde’lyi et al (1954, p. 129, Vol. 2). 
However, since it would have been difficult to confirm the validity of 
(E-13) for Im(k) equal to zero using this approach, ($-13) was established 
through (E-l). 



